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PREFACE 


On May 15-16, 1975, a symposium was held at the 
Goddard Space Flight Center on the Study of the Sun and 
Interplanetary Medium in niree Dimensions. The symposium 
brought together more than 200 scientists from the U. S. 
and Europe to discuss the importance of e}q>loring the inter- 
planetary medium, and viewing the Sun over a wide range of 
heliographic latitudes. Among the topics discussed were the 
missions that NASA and ESA are currently considering for 
possible flight out of the ecliptic plane, and the likely 
scientific returns of these missions in the areas of solar, 
interplanetary and cosmic-ray physics. The symposium was 
sponsored by NASA and ESA, and was a topical meeting of the 
American Astronomical Society and the American Geophysical 
Union. 

In this proceedings an attempt has been made to pro- 
vide a complete summary of the content of the symposium. 

J. A. Simpson has kindly summarized the various options for 
out-of-the-ecliptic missions. The speakers who addressed 
the likely scientific returns from these missions have in 
most cases provided detailed summaries of their remarks. 

In cases where a full summary was not available, an abstract 
has been provided. In addition, K. C. Hsieh has kindly 
provided an article on direct measurements of neutral gas 
out of the ecliptic. Time was not available at the symposium 
for presentation of these interesting ideas. 

iii 



The organizing committee of the symposium is very 
grateful to the authors, all of whom have busy schedules, 
for providing detailed summaries of their talks. We would 
also like to express our sincere appreciation to Mrs. Sandy 
Schraeder and Mrs. Martha Harding for their invaluable help 
in organizing the symposium. 

A list of all the speakers at the symposium, and 
their topics is given in Appendix A to this proceedings. 

A list of all the attendees at the syiq>osium is provided 
in Appendix B. 

L. A. Fisk 
W. I. Axford 
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Some Introductory Remarks 


Over the years studies of the Sun and interplanetary 
medium have occupied the attention of a sizeable fraction of 
the space community, and have constituted a considerable 
portion of NASA and ESA's efforts in space science. And 
rightly so. We attempt in this subject to understand our 
local environment in space. 

With our continuous and in some cases in situ obser- 
vations, we also provide in this subject detailed testimony 
as to v;hat physical processes are possible on a star and in 
an astrophysical plasma. Our studies of the solar environ- 
ment thus influence our thinking as to what is possible in 
other astrophysical settings. In some cases the influence 
is quite direct. For example, the solar wind should exhibit 
many of the properties of stellar winds. Shock waves in 
interplanetary space may resemble those in interstellar 
space. The behavior traits of cosmic rays propagating in 
the interplanetary magnetic field should be similar to those 
of cosmic rays in the galaxy. And so on. 

All of the spacecraft that NASA and ESA have flown to 
study the Sun and interplanetary medium, however, have been 
limited in one major respect. None of these spacecraft have 
penetrated off the equatorial plane of the Sun by more than 
about + 10° in heliographic latitude. We have thus sampled 
particle emission from the Sun, or interstellar matter 
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impinging on the solar cavity (e.g. galactic cosmic rays), 
only over a narrow range of latitudes. The look-angle for 
studying photon emission from the Sun is similarly limited. 

This latitude limitation would not be serious, of 
course, if we had any expectation that conditions in the 
solar environment were invariant with latitude. However, 
even the most cursory of examinations of a white light 
photograph of a solar eclipse reveals a strong latitude 
dependence for conditions in the corona. Sun-spot and 
flare activity is known to be concentrated in the mid- 
latitude regions on the Sun. Interplanetary scintillation 
studies suggest that the flow of the solar wind is faster 
and more turbulent over the solar poles, than it is in the 
equatorial plane (See Coles and Rickett, in this proceedings). 
The expected Archemedes spiral pattern of the interplanetary 
magnetic field may cause the ratio of thermal to magnetic 
pressure of the solar wind to vary strongly with latitude, 
thus resulting in a plasma over the poles which has substan- 
tially different properties from the local solar wind. This 
same field pattern may let low energy galactic cosmic rays 
( '' 100 MeV ) penetrate unopposed into the region of the 
solar poles, whereas such particles are excluded from the 
inner solar system near earth. And the list goes on. 

We are thus forced to conclude that, to date, space- 
craft exploration of the interplanetary medium, and measure- 
ments of the Sun which are look-angle dependent, have studied 
a non-representative sample of our total solar environment. 
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We have studied in only two dimensions, what is a three- 
dimensional structure. We have not sampled the full variety 
of astrophysical conditions that are available to direct 
measurements in our local environment. Thus, we have not 
adequately broadened the base of knowledge that we can use 
for deciding what is possible in other astrophysical settings. 

An Out-of-the-Ecliptic (0/E) Mission, in which the 
interplanetary medium is e}q>lored, and the Sun is viewed 
over a wide range of heliographic latitudes, will provide 
the measurements that can help replace our current parochial 
view, with a more accurate assessment of our local environ- 
ment in space. 

The Current Status of the 0/E Mission 

At the time of the publishing of these proceedings, 
NASA and ESA are engaged in a Phase A study of an 0/E mission, 
which could be launched in the early 1980' s. Two possible 
missions, a primary and a backup, are currently under con- 
sideration. In the primary mission two spacecraft are to 
be launched simultaneously by the shuttle, with the Interim 
Upper Stage (lUS) - 4 stage booster. The spacecraft, which 
would be similar in design to Pioneers 10 and 11 , would then 
fly to Jupiter. There, one of the spacecraft would be 
targeted so that after encounter it passes up out of the 
ecliptic plane and over the north pole of the Sun. The other 
spacecraft would be targeted so that it passes over the south 
solar pole. Following their polar passes, each of the space- 
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craft would return to the ecliptic plane, and then would fly 
up over the opposite solar pole. In the backup mission, one 
spacecraft is launched from the shuttle /lUS - 2 stage. This 
spacecraft is targeted at vTupiter to pass over one of the 
solar poles, return to the ecliptic plane, and then fly up 
over the opposite solar pole. In both missions the space- 
craft can obtain high heliographic latitudes ( 80^). 

The major advantages to the primary mission are: 

(i) It will provide simultaneous measurements of conditions 
in the northern and southern hemispheres of the Sun. (ii) The 
total science payload (both spacecraft) carried in this mission 
is large ( 60 kg), (iii) The two spacecraft passing by Jupiter 
will provide an interesting opportunity to distinguish spatial 
from temporal effects in the Jovian magnetosphere. Indeed, 
this mission may provide the only opportunity in the forseeable 
future when two spacecraft are simultaneously near Jupiter. The 
advantage for the backup mission is, of course, in cost. It 
costs less than the primary mission, and yet it still provides 
a survey of conditions at high heliographic latitudes. 

In the primary mission, it is planned that ESA would 
build one of the spacecraft, and NASA, the other. In the 
backup mission, the spacecraft would be provided by ESA. 

In both missions, NASA would provide the launch vehicle and 
the RTG power supplies. 

The Phase A study will be completed in the spring of 
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1976. After this time, NASA and ESA will decide whether or 
not to proceed with a joint 0/E mission. A decision can be 
expected in the summer of 1976 as to whether a Phase B study 
of one of these mission options (or some other option) will 
be undertaken. 

Details on the primary and backup mission, as well as 
on other options for 0/E missions, can be found in the review 
paper that J. A. Simpson has kindly written for this proceed- 
ings. It should be noted also that a particularly interesting 
variation on the 0/E mission is discussed in the proceedings 
in the paper by G. Colombo, D. A. Lautman, and G. Pettingill. 
In this variation, one of the two spacecraft in the two- 
spacecraft, Jupiter-swingby mission is a solar probe, which 
after encounter with Jupiter is directed back at the Sun 
along a nearly rectilinear path. Accurate measurements of 
the quadrupole moment of the Sun can be made from the solar 
probe, as well as in situ measurements of the solar wind in 
the corona. 

0/E Science 

Some information on the flow speed of the solar wind 
at the higher heliographic latitudes can be obtained from 
observations of interplanetary scintillations, and also from 
observations of ion comet tails. In this proceedings we 
have included a paper on scintillation observations by 
W. A. Coles and 3. J. Rickott, and one on studies on comet 
tails by J. D. Brandt. It is interesting to note, however, 
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that the conclusions of these two papers on the behavior of 
the average wind speed are substantially different. These 
differences may be reconcilable, as Coles and Rickett suggest. 
However, such discrepancies in indirect observations emphasize 
the need for direct in situ measurements. 

The remaining papers and abstracts in this proceedings 
discuss the impact that an 0/E Mission will have on studies 
of the solar corona, solar x-ray and EUV emission, solar 
radio astronomy, the solar wind, the interplanetary magnetic 
field, solar and galactic cosmic rays, interplanetary dust 
and zodical light, and interstellar neutral gas. It is 
recognized immediately from this list that 0/E science is 
highly multi-disciplinary, involving many areas of solar and 
cosmic-ray physics, and all areas of interplanetary physics. 
There are also many different topics in each area which will 
be affected by 0/E measurements. It is not the purpose of 
this proceedings to give an exhaustive survey of all possible 
topics. Rather, we have simply highlighted here some of the 
more interesting problems in 0/E science. In this sane 
context, it should be recognized that an 0/E mission is 
exploratory, and as such is certain to uncover new phenomena. 
It is one of the axioms in the space program that when we fly 
where we have never been before, we uncover phenomena that we 
cannot anticipate in advance 

The papers and manuscripts in this proceedings have 
been reproduced directly from the typescripts provided by 
the authors. In rddition to considering scientific qiaentions , 
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the authors In some cases have discussed the hellogiaphic 
latitude that an 0/E mission must obtain for best results in 
their particular area. NASA and ESA requested that such 
information be provided in the talks at the symposium. 
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1 ^ 76-24120 

Experimenh Out of the Solar System 
Ecliptic Plone; An Introduction 
to the Execliptic Mtssion* 


J. A* Simpson 

Enrico Fermi institute and Department of Physics 
University of Qiicogo 
Chicago, Illinois 60637 


February 28, 1976 


U Introduction 

The dramatic step made possible by direct measurements in space with soteMites 
and probes during the past 18 years has totally altered our concepts of the Sun, the 
interplanetory medturrvand their influences upon Eorth. This has been achieved with 
observations confined solely to the vicinity of the equatorial plane of the solar system. From 
these two-dimensional investigations we have mode dubious attempts to extrapolate our 
knowledge to deduce what the Sun and space in the solar system is like in three-dimensions, 
However, the solar, interplanetary end galactic phenomena discovered in these years hove 
roised many urgent scientific questions which can only be answered by direct observations and 
experiments far out o\ the ecliptic plone and over the solar pole to achieve a "global" concept 
of the Sun, the interplanetary mediunvand their relationship to Earth and the boundory of the 
heliosphere with the Interstellar medium. We have been faced for many years with this age-old 
problem which occurs often in science, namely, the extrapolation of physical phenomena from 
two-dimensions to deduce phenomena in three-dimensions. The uniqueness and importance of 

♦ Based upon invited talk at "Symposium on the Study of the Sun and Interplonetary Medium 
in Three Dimensions" — May 15-16, 1975, Goddard Space Flight Centei, Grccnbelt, Maryland, 
and "Charged Particle Astronomy in Interplanetary Sparc ot High Solar Latitudes", J. A. Simpson, 
Chap. 5,3; NASA collected papers for ESRO Document MS(74) 34 (1974), 
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o scientific mission which can directly achieve this global study and the recognillon of its 
potential for discovery has been clear for many years. In describing such an exploratory 
mission, it is not unfair to moke a comparison between the importance of Man's exploration of 
the spherical surface of the Earth, and on execlipHc mtssiory which is qualitarively 

similar in its conceptual and practical consequences for space science to the impact of the 
full exploration of Earth on Man's intellectual advancements. 

If a mission out of the ecliptic is so vital to the advancement of science then 
why has it not become a reality by now since the technoiog'' for its accomplishment has been 
with us for several years, and Pioneers 10 and 11 have demonstrated that a Jovian gravity- 
assist to drive a probe out of the ecliptic is sofe? Among die reasons appears to be the broad. 
Interdisciplinary character of the most important investigations on the mission which, on the one hand, 
represents the greatest strength of the mission, but, on the other hand, becomes a source of 
weakness for marshalling the sources and the support of the scientific community, or even 
leadership within the federal agencies where the missions must successfully compete with other 
important types of space missions. Since the late 1950*s an out of the ecliptic mission has been 
under discussion. ^ Now, hopefully, it is a mission whose time has come since it has recently 
elicited support from a wide segment of the scientific community on the basis of its uniqueness 
and importance for science ond the applications of science to our understanding of the Sun and 
its influence upon Earth. Furthermore, as has been the case for the space program to date in 
the equatorial plane, out-of-the-ecliptic observations are almost certain to yield important, 
unanticipated discoveries. We can best describe the quality of the mission objectives as 
exploratory and interdisciplinary and, therefore, the investigations must be designed tc encompass 
the unexpected. 

The purpose of this note is to summarize the most likely alterna rives for carrying out 
a mission to achieve these broad scientific goals and to illustrote wirfi specific examples d*awn 
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from charged particle astronomy. Interplanetary and solar physics some of the experiments and 
observations which may be carried oi*r. It Is not within tlie scope of this note to describe In 
detail the many exciting scientific challenges opened by the mission, but the reader easily will 
perceive this wide range of posiible investigations, many 6f which are discussed in the 

2 3 . . . 

Proceedings of this Symposium or outlined by Page. As a basis for discussion of the regions 
of the Sun and interplanetar}' space which may be explored byexccliptic missions I have 
prepared in Figure 1 a schematic represcniotion of the main solar latitude zones of highest 
interest and their possible interfaces with the interplanetary medium. In the fol!ov;!ng discussion 
we win summorize some of the alternative missions which reach into these regions of space and 
the constroints they place upon experiments. As a secondary objective we sholl ol;c> describe 
the opportunities which some of these missions provide for unique studies of the magnetosphere 

of Jupiter since, for the most likely mission choices, Jupiter becomes 

the *'gofeway'' to space out of the ecliptic, 
li* Ten Ways to Tet Ttiere 

Tuble ] is a summary of the most outstanding tnisslon alternarfves. There is 
a h-ng Iilsloiy o" -jccessivc proposols for these mlssiunr bored primarily on avuMc4ble launcli 
vcincle techno' ogy, hence some of the a! iei natives picsented in Table 1 nov; are ordy of 
historicoi int. r ’sU Basically, llie mission options are dependent on laundi veliicle capobi I it ics. 
Direct balli.-ilc injection ai ~ 1 a. u. (option 7 and C) even under optimum conditions will tcJo3 
c spacecraft only to ~ 37*^ hcHographic latitude. Ihc ocidition o! solar electric propulsion fSt P) 
makes it possible to achieve a spacecraft trajectory oscillating iti latitude ot ) a. u. , and 
synchronous with Earth thus scanning over a north^south latitude ror>gc that reaches a limit of 
t < 60^ within 2, 7 years, under maximum launch capobi litics (option 9, for dclotls sec 
rcfcirnce 4 and 5), Tiu: trajectory plotted in a plonc at 1 o. u, is sliown in figuic 2. This ty|.)o 
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bbmgrinadevqKEqrln'Ibstorbs^ 


E iir^e^ten^^ V|be irate teebnieol duvelopmerit of k>lor eteefrtc pn:^uisSoa> 

? jri^ the early !900'|. 

^ from Figure I see thbt olf of ibe riirect infection mluion^ discussed above do not 

■fij^bW ebseiyqHons iii si tu beypftd the zones of solar activity and leave (^en all qoesHotis vs^lich 

:^(»u!d only be answered by attaining higher solar lotitudes. As will beshown below, the po lor 

:, !|^*to“polar cop passes of Jupiter swingby missions will achieve the latitude scons of the StP 

pylons, aldiough for a shorter intervol of the solar cycle. 

The alternative to direct infection is a gravity assist by Jupiter, i. e. by a 

Jupiter swingby (JSB) — a technique first proven by Pioneer- 10 to achieve a soioi system 

6 

escape trajectory. This technique makes it possible to achieve solar polar cop posies with 
thdximum latitudes depending on the launch constraints (opiions 1 6). Among those options 

that achieve polar cop posses of > + 80^, options 2, 3 and 6, offer the greatesi potential 
for achieving our stated prime objectives of fixp^lorotioo and discovery. Option 6 enables a single 
Spacecraft to reach 79^ solor latitude v/ilfi an Atlas/Ccntaur vehicle by requiring on Earth 
swingby o$ shown in Figure 3. However, in add! lion to the increased time to reach solar maximum 
latitude (1. 7 years longer), this mission suffers further from reduced reliability because it 
requires on additional spacecraft propulsion subsystem to undertake tv/o additional and critical 
spacecraft maneuvers as described in Figure 3. We focus on the most fruitful of all the JSB 
options, hopefully the most likely to be adopted, namely a dual spacecraft launch to Jupiter 
(options 1 or 2) which will result in spacecrafts over both solar polar caps simultaneously with 
tre Rectories as shown in Figure 4 passing fiom pole-to-polc in opposition, Wc refer to this type 
of misiton os the tandem Jupiter swingby (TJSB), 
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alt at a fiifiatlaa of tlma* AiHratqih ittt« Mtampla U Ibt option 2 In Tobta the tta|aafOnai-^ 
flpHon 1 wtli ba tlm^ from tb< point of view of tclantiflc InvastisoHons fti different * , 
;^O^I{inat^ Oig. IntarplonetoryplasnM ond imigneHe flald, tbe mognetoailiara 

^-Jupllar, eosinie ro;irs« ate., tbo miulon Is topicaily to ba dlvldad lnio five pheiBs {derttlfled 
by portiom of tbe spaceeroft tro{ectorles. Tbesaore: (1) from Earth to Jupiter, (iD lhrough 
;lha Jovlort meQnetosphare, (III) out-of-the>eciiptic slmul^cmeousty In 4te northern and southern 
hemi^eres over the radial range ~ 1.5 • 5 au., (IV) over the solar pole, and pole-to* 
pole trandts of the ivro spacecrafts, and ^/) post solar pole tra{eelories. In each phase 
diere ore some prime mission goals for one or more of the scientific investigation^ In the 
following discussion we take the reader on a "guided tour" through these five phoses of the 
TJSB missions using illustrative scientific investigations which will lead to discovery or the 
answer to old questions. Although it is not possible to discuss in this note all the important 
scientific (^jeetives of each phase, these trajectories provide a rich source of new investigations 
with each phase of the mission possessing its own set of unique scientific objectives. 

Phose (I): Earth to Jupiter 

Ihe two spacecraft travel near the ecliptic plane with a radial-spatial separation of 
order 10^ kilometers and with simultaneous transmission of data. This separation makes it possible 
for us to undertake a new family of studies, since never before have spacecraft been so separated 
for a long period of time free from the influence of a nearby planet and never before have we hod 
the opportunity to do correlative studies between closely spaced observatinn points over a large 


radial range. For example we may undertake: 
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a) the study of cliargod particle-magnetic field interactions, especially 
for very low energy nuclear particles in die range 0. 1 to 1 MeV. 

This spacecroft separation distonce becomes comporable to the correlation 
length of the interplanetary magnetic field and to the scattering scale 

g 

size of the particles. 

b) the study of the modes of propagation and interaction with magnetic 
fieicbin interplanetary space of the electrons which have been recently 
found to be escaping from Jupiter. 

12 

c) meosurements of the interplanetary occelerotion of protons and electrons 

13 

in the regions surrounding blast waves from die sun. It will also 

become possible to investigate in detail the forward-backward moving shocks 
which are now observed to b*» associated with so-called ** interplanetary 

14 

active regions*' . These active regions and shocks are also associated with 
enhanced fluxes of 1 MeV protons. 

This phase of the mission corresponds to on interplanetary version of the smaller 
scale Mother- Dough ter satellite combination devoted to mognetospheric studies in the period 
1977-1980. No other interplanetary mission studies of the above type hove been mode, or are 
contemplated in the foreseeable future. 

Phase (II): Jovian mognetospheric studies 

Since for operational reasons the two spacecrafts in opposite hemispheres will 
have timesof closest approach 2 to 3 days apart, we obtain a unique and valuable separation 
of the two spacecrafts in the Jovlan magnetosphere capable of attacking problems that could 
not be investigated by the Pioneer 10-11 spacecraft, the Mariner- Jupiter-Saturn spacecraft, or 
even a single Jupiter orbiter spacecraft. In Figure 5 we disploy o. meridional plane projection 
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of iIk* traj.^ctories of fhe tv/o spacecrafts A and B. Figure 6 is a projection of the spacecraft 
trajectories on the ecliptic plane. For comparison the Pioneer-11 encounter trajectory 
to Saturn is shown. At the time of closest approach for spacecraft A (position 1), spacecraft B 
is at a distance of 50 Rj, and when spacecraft B is at closest approach (position 10) spacecraft 
A hcis moved to 50 P^. Thus it should be possible to separate large scale spatial from 
temporal effects in the Jovian magnetosphere. It will also be possible to obtain measurements 
at four magnetic latitudes for each radial distance. Some of the key problems to be attacked 
are: (u) investigations of the varluiion of the radial position of the bow shock witi; tii* e, and 
(b) the distortions of the mognetospherlc boundar)^ in response to fluctuations in the strength and 
direction of the solar wind and the rotation of the magnetospiiere. An Important feature of the 
Jovian mognetosphoric observations possible with the TJSB mission Is the simultaneous 
measuremant of tho solar wind outside the magnetosphere by one spacecraft while meosuroments 
wiihit'. the mogneto:.nheie ore under v/ay with the second spacecraft, (c) The nature of the 
"global time depend .-nt 10 hour variations of electron intensiiy and spectrum within the 
magnetosphere^' ' how is this effect related to the rotation effects of the equatorial plasma 
sheet end the rpot'aMy dependent lO hour variation? This in turn is reloted to the problem or 
the mechanism foi the release of electrons from Jupiter into interplanetary space, (d) Jovian 
satellite interactions with the trapped radiation; «peclal opportunities exist whereby it is 
possible to cross the flux tubes associated with the satellite lo and thus to investigate the nature 
of the control excr led by lo over decametrlc rodio bursts. ^ 


In sum, the dual spacecraft out-of-the-cdiptic mission will answer questions 
whicli would otherwise remoin a puzzfe for studies made with a single spacecraft, 

Obsoivations made with Pioneer 10 and 11 have eslabliihcd the importonce of 
transient phenomena ' ^ for tho magnetosphe re of Jupiter, such as largo scale disforllo-if. 
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Observations v/lth a single spacecraft cannot unambiguously scporate the temporal and spatial 

dependences of such transient effects, so that a dual spacecraft mission offers our best hope for 

gaining o further understanding of the physics of the Jovian magnetosphere. 

Phase (III): Out-of">the-"ec l tptic at lorge radial dlstonces 

Figure 4 illustrates the trajectory characteristics for spacecrafts A andB as 

a function of time after a Jovian swingby. In a period of ^ 2 years, the tv/o spacecrafts 

traveling in the opposite hemispheres of the interplanetary medium cover a rodial distance of 

4 a.u. while slowly traversing a solar latitude range of up to nearly 90^. 

This phase of the mission offers the opportunity to take snapshots of solar 

active regions on the sun, (EUV, UV, x-rays, radio, etc.) for comparison with ldf=mtlcal 

observations from Forth to form stereoscopic pictures of solar phenomena. It also becomes 

2 

possible to Investigote new aspects of the Gegenschein. 

This phase of the mission also offers the opportunity to study the behaviof 
of magnetic sector structure at large radial distances from the Sun in the solar activity zone 
(Figure 1) to answer such questions as: (u) V/hot role does th<^ region cf solar activity (10 to 
35 degrees) ploy in determining the sector structiiro of magnetic Ticids extending From 1 to 
5 a. u. ? (b) To what extent does tfie magnetic sector structure persist at high latitudes and at 
great distances from the Sun? (c) How does the Sun‘sdifferential rotation, which produces a rotation 
period 9 days longer at the pole than at the equator, change the structure of interplanetary megne tic fields ? 
These and questions concerning the"global shape" of blast waves in the two hemispheres constitute 
major magnetic field and plasma studies for this phase of the mission. The characteristics of this 
region for charged particle propagation for both solar and galactic particles are entirely 
unknown, and their determination would be the pilmc goal of clicirgod particle studies in this 
region. 
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Phase (IV); The polqr observations of Sun from ~ K 2 ** 2 o» 

A. Solar observaHons from i!ie polar viewpoint 

The evolution of coronal feoturcs above solar active regions, coronal streamers 
and related transient, large scale phenomena now observoble only by solar limb studies 
from Earth -- may be undertaken from the A B spacecraft by time-lapse observations 
obtoir>ed simultaneously at all solar longitudes. Coronagraphic studies on a spinning spacecraft 
are difficult, but the potential is great for understonding the origin and dynamical structure 
of the inner and ot’er corona from simultaneous, polar and equatorial observations (J« A. 

Simpson to G* Nev'^ k, private communication, 1968). 

B. Solt interplanetary studies* 

Somewhere in the region tentatively identified as the transition region above 

solar latitude 60^, characteristics of the interplanetary medium increasingly become 

determined by the properties of the sun and corona in the polar regions. For exomple, it is 

currently believed that the polar region may be represented by a coronal hole where a 

continuous emission of the solar wind at high velocities (> 700 km/sec) is expected. It is in 

this region that the rotational effects of the Sun cease to play a major role in the large-scale 

structure of the magnetic field carried into interplanetary space by the solar wind. The properties 

of this region are unknown and expected to be totally different from those so far studied in the 

vicinity of the equatorial plane. For studies in this region it is vital that measurements in the 

north and south polar regions be made simultaneously since it is well-established thot both the 

temporol and spatial distributions of observable solar phenomena in the polar regions are frequently 

18 

different at the two poles. In oddition to the particle, magnetic field and plasma interactions 
which will be studied for the first time under these new physical conditions, we point out that 
low energy particles from the galaxy moy find a relatively easy entry to this 

. U.i: UV or' THO 
i 'GO IS 
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r« 9 )on, as discussed below. 

As shown in Figure- 4, the two spacecrafts pass from pole to pole in periods the 
order of 260 days (or on on overage of ~ 0. 7^ per day). This corresponds to an elapsed time of 
~ 10 solar rotations when the two spacecrafts are at radial distances 1.2 to 1. 5 a. u. Cleorly 
radially dependent effects are likely to be small compared with latitudinaily dependent 
phenomena being studied simultaneously on the two spacecraft. Hie observations to be made 
and the scientific objectives of Phase IV of the dual spacecraft are essentially the same as those 
for the SEP mission (options 9 and 10). Both missions provide a scan of solar latitude at a 
rate of ~ 0. 4 - 0. 7 degrees per day at approximately constant radial distance from the Son, 
although the dual spacecraft mission provides coverage of the polar region of the Sun, while 
the SEP mission does not. Among the latitude dependent phenomena to be investigated are 

(a) the effect of differential rotation on the magnetic field structure in interplanetary space; 

(b) the nature of the transition region from the polar coronal holes to the band of solor activity 
(Figure 1); (c) the nature of transient phemomena such as shocks and high velocity streams 

at high solar latitudes. Measurements relating to these questions made in Phase IV are distinguished 
from similar measurements in Phase III by the fact that the radial position of the spacecrafts is 
not an important parameter during Phase IV, thus providing a clean separation of latitudinal 
and radial effects. 

It is likely that the combined direct solar observations, magnetic field, plasma 
and high energy particle studies will introduce a qualitative change in our understanding of 
the differential rotation of the Sun, of the 22 year magnetic cycle and, thereby. In our 
understanding of the internal dynamics of the Sun. 

C. Galactic composition of cosmic radiation. 

In part B above, measurements during the pole-to-pole excursion of the 
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two spacecrafts were concentrated on the electrodynamics of the interplanetary medium 
and the rote of the solar features in determining the dynamics of the medium. If conditions 
over the solar poles are anywhere neor those predicted !t would appear that cosmic ray particles 
of low energy from the galaxy which cannot otherv^ise propagote into rhe inner part of the solar 
system near the equatorial plane becaus; of solar modulation may be oble to penetrate 

by way of solar poiar magnetic fields to within ^ i - 2 a. u. If so, we may be able to 

obtain for the first time samples of th ; l energy composition of galactic cosmic rays; 
that is, the relotive abundances of tite elements in the nuclear component of the 
cosmic rays and the relative isotopic obundances of hydrogen to nickel. 

Through such studies it may become possible to identify the low energy component of cosmic rays 
accelerated In our local region of the galactic arm. These studies will be of vital Importance 
for deciding among models of nucleosynthesis of the elements in the sources of cosmic rays. 
Furthermore, under such circumstances, it v/ould become possible to obtain the energy densities 
in interstellar space for these very low energy particles ( a problem concerned with the heating 
of interstellar clouds). 

Finally, all of those investigations over the solar pole v/hen takei together with 

observations in the equatorial pL’ne will yield a **global'‘ model for solar modulation which 

takes account of the propagation of nuclear particles and electrons extending dov/nward in energy to 

19 20 

energies where at present their modulation is not understood. Recent observations ' shov/ 

that evenathigher energies, revisions of our ideas about modulation may be required, possibly 

21 

involving processes taking place off the ecliptic plane , or Introducing iiverstcllor neutral 

22 

particles in the heliosphere. The IMP satellites and Pioneer 10/11 deep space probe observations 
have raised a number of interesting questions regarding whether or not low energy particles 
could have occcss to the solar equatorial zone. 
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D. Energcifc portlclcj* of solor origin. 

Iho role of coronal Ironiipoi i in the pio|joyotion of solar poiiidc*-. from More 

23*^26 

sites to the interplanetary medium has been much studied but is not v/ell understood. 
Out“of-the-cdipl5c missions will play a major role in deciding on tf»e transport mccharMsm, 
on the storage lime ond distribution of porticles at ihc sun, ond, in turn, v/hile using solar 
particles as probes of the intervening magnetic fields, will obtain information on fheneai-sun 
magnetic field structure including the distribution of irregularities in the magnetic field. It 
may be possible that the effects of differential rotalion can be onalyzrd best by studying t!ic 
emission of solor flare particles at high solar latitudes. 

E. Models of the heliosphere. 

At the present time we cannot choose conclusively among models of the 
heliosphere with boundaries for particle modulation which, for example, could be a) pancake- 
like in character, extending the order of soy 20 to 50 a. u. in the equotorial region, but only 

a few ostronomical units thick over the poles, or b) with distant boundariics over the solur 

27 

poles and therefore much more spherical in character, os sketched in Figure 7. Although 
these two examples represent extremes, it would appeorthat data obtained out-of-lhc- 
eciiptic plane could assist in deciding between thern by using galactic cosmic rays as probes 
of the outer magnetic fields of the interplanetary medium. Predictions for th harged porticlc 
gradients, particle anisotropies ond energy spectral changes as a function of solor la ‘^ude can 
be made; therefore a wide range of models can be tested by the dual spacecraft mission 
(e.g. reference 28)* 

Phase (V): Pest solor polar observotions 

The above four phases of the dual mission iMustrate the wide range of physical 

processes which can be studied during the mission. After leaving the sun Hie dual spacecrafts 

again travel outward from the Sun where they are able to repeat some of the observations winch 

were obtained between Jupiier and Iho Sun in Phase III us mucli os a half s olnr cycle earli er. 
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Ittese loiter meosuremeni$ would Irdccd be very intcresiing since they con provide further 
evidence on the long-term chatiges in the heliosphcrt^ especially in tite latitudinal structure 
of the interplonetory mcgnetic fields. 

IV. Summory Remoris 

We conclude from the mission alternatives that the latitude scon missions 
such os the SEP options 9 arKl 10 (Figure 2) are dedicated to more detailed observations 
ond exploration of the soiar active zones (Figure 1) with o strong emphasis on extending 
the stereoscopic viewing of the solar phenomena now under observation with Earth orbiting 
satellites. On the other hand, the solar polor cap passes further extend the exploration 
of new regions of the Sun and interplonetaiy medium, the evolution of solar coronal features 
seen simultaneously otoli solar longitudes, and the possible access of low energy particles 
from interstellor space via the polar magnetic fields. Such particles can not be detected 
with deep space probe missions near the equatorial plane in the foreseed^le future. Tlius, 
in many aspects both types of missions are inqjortont for science, but are qualitatively differen 
in their goals. 

With regard to strategies, the TJSB missions offer the unique advantage of 
providing well-defined interfaces for international collaborations since one spacecraft could, 
for example, be the responsibility of the European Spoce Agency while the second spacecraft 
could be the responsibility cf the NASA. 

Some instruments sliould cover the same measurements simultaneously on both 
^acecrofts, e.g. magnetic fields, plasmas, charged porticics ond x-rays. However, in 
addition one spacecraft could carry a set of complex instruments to complement the other 
spacecraft, viz, a coronagraph on one spacecraft, and a complex super-thermal particle 
spectrometer and solar radio emission detec.or c the other spacecraft. 

n .,- rJCIBlLITY OF THE 
0UiOLwVI> PAGE IS POOR 
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AllhoMsb Tdbie 1 suggests o wide range in costs bosed upon the Icnmch 
vehicle, the ultimate difference between o single ^cecroft JSB mission and TJSb mission 
is reoll^ less than indicated beccAise the same imqpsitude of commitment otkI resources is 
required for scientific instrument preparation and integration, for data acquisition throughout 
the yeors of the mission and, most iipportont, for the level of commitment of those in the scientific 
community motivated to undertake such o lortg term enterprise. 

Ihe need for simultaneous measurements at Earth during an execliptic mission 
rnmst not be overlooked in order to separate spatiol from temporal changes in solor inteplanetary 
phenomena, and to relote these observations to the present day scientific knowledge derived 
from equatorial meosurements. 

Ingenious experiments and observations have been reported for many years 
to explore the high solar latitudes near the Sun and interplanetary space. These include 
die use of radio waves from distant stars to study the magnetic irregularities and electron densities 
near the Sur^ the observation of comet tails at high latitudes ond the scintillation effects 
of galactic cosmic rays to deduce properties of the solar wirid, ond the large-scale probing 
of the interplanetary medium by high energy cosmic rays to estimate tlie scale size of the 
heliosphere. However, they cannot substitute for direct observations in the regions of the 
solar system to be penetrated by an execliptic mission. 

The author apologizes for not adequately covering in this note the many alternate 
mission options with their unique scientific objectives. 

He wishes to thank Or. Bruce McKibben and Mr. John Niehoff for 
assistance with the preparation of materials for this manuscript, and Mr. Don Herman 
of NASA and Mr. H. F. Matthews of the NASA/Ames Research Center for providing essential 
background material. This work was supported in port by NASA grant NGL 14-001-006. 
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Figure Cop<ion$ 


Figure 1 


Figure 2 


Figure 3 


Idealized mcriodional plane view of the interplanetary regions associated 
with principal features on the sun. Ihe shaded region represents a region 
-t- 7^ in solar latitude within which all measurements to date have been made. 
The region of principal solar activity over o solar cycle extends from 10 to 
35-40 degrees north and south latitudes and is highly variable. A region 
from ~ 40 to 70 degrees is a transition region betv/een the region of solar 
activity and the polar region where the rotational effects on the magnetic 
field carried out by the solar wind begin to subside. It is believed that the 
polar region is mainly occupied by a coronal hole-like structure and therefore 
that the solar wind has a high velocity in this region. 

The solar electric propulsion mission is one in which the spacecraft remains in 
a 1 A. U. orbit from the sun and is therefore synchronous with earth. The orbit 
inclination is increased by thrusting about the nodes. In a period of ~ 3 1/2 to 
4 years a full excursion of the spacecraft is expected to be between 50 and 60 
degrees. The spacecraft is normally thrusting except for ~ 100 days per year 
at the anti-nodes. However, the spacecraft propulsion can be turned off for 
a day or so during the normal operating periods lo obtain scientific data. 

(See references 4 and 5. ) 

The AVEGA flight mode (AV-Eartli Gravity Assist). In this mode the 
transfer event points are: 

Point 1: Earth laurich 4/18/82, = 27 Km^-sec 

Point 2: Perihelion modification maneuver, AV~ 900 m -sec \ 
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Figure 4 


Figure 5 


Figure 6 


Figure 7 


Point 3: Earth powered swingby, AV -■ 1000 m-sec . 

Point 4: Jupiter encounter 3.2 years after launch. 

The trajectory is not to scale. 

The out-of-the-ecliptic trajectory for the dual mission after Jupiter swingby# 
showing the radial distance of the spacecrafts from the sun and the 
heliographic latitude of each spacecraft as a function of time. (Adapted 
from reference 7. ) 

The dual spacecrafts A and B enter the Jovian magnetosphere approximately 2 
or 3 days apart. The figure is a meridional projection of the spacecraft 
trajectory with Jupiter at the center of the coordinate system. Tlie fiducial 
marks on trajectory A represent 6 hour Intervals v^hich correspond in numbers 
to the 6 hour intervals along trajectory 6. Thus it Is seen that one spacecraft 
Is near closest approach when the other spacecraft is at 50 Jovian radii. 

The trajectory of Pioneer 11 is shown for comparison. 

Projection on the ecliptic plane of the trajectories of spacecraft A and 
spacecraft B. For comparison the trajectory of Pioneer 11 is sliown. 

Two alternate models for the shape of the heliosphere. (See reference 27. ) 
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Idealized Meridional Diagram of Solar Regions 
Connecting with the Interplanetary Medium 


Figure 1 
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AN ALTERNATIVE OPTION TO THE OUAL-PROBE OUT-OF-EaiPTIC MISSION 

VIA JUPITER SWINGBY 


6. Colombo 
and 

0. A. Lautman 
Center for Astrophysics 

Harvard College Observatory and Smithsonian Astrophysical Observatory 
Cambridge, Massachusetts 02138 

and 

G. Pettengill 

Massachusetts Institute of Technology 
Department of Earth and Planetary Sciences 
Cambridge, Massachusetts 02140 


We have recently conducted a preliminary study on the possibility 
of combining the out-of-ecllptic (OOE) mission with a solar-probe 
mission. In particular, we have been looking at the possibility of 
having a high-incllnation OOE probe complemented by a second probe 
going from Jupiter to the sun along a rectilinear path (at least for 
the segment from C.3 a.u. Inward to the sun). 

The scientific Interest In approaching close to the sun is obvious 
since it enhances observation of particles, fields, and gravitational 
harmonics. Our particular choice of path results from the associated 
simplicity of the spacecraft configuration needed to provide, for 
example, good thermal control, a drag-free system, and good communica- 
tions with the earth. 
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A prelialnanr error analysis conducted by J. D. Anderson of the 
Oet Propulsion Laboratory leads to very Interesting conclusions for 
an elliptical ortit with a perihelion distance of 16R^. Assuming ^at 
the nongravltatlonal forces are compensated by a drag-free system 
(with three degrees of freedom) and that the spacecraft Is tracked 
down to perihelion, the quadrupole moment of the sun can be determined 

7 

with an accuracy of 3 parts In 10 . Since the estimated value of J 2 
ranges from 3 x 10“® (applying Oicke's theory) to 1 x 10“^ (assigning 
rigid rotation of the Interior with the observed surface), the Inter- 
est In determining this moment with an accuracy of at least 1 part In 
10^ Is clear. He remember that ^2 gives & fundamental constraint to 
the moment of Inertia (or the ratio O/NR ) and, therefore, on the 
Internal density distribution of the sun. 

As mentioned above, the result obtained by Anderson implies a 

-8 2 

three-axis drag-free system with an accuracy of 10 cm/sec . A drag- 
free system having this accuracy has recently been flown in the 
TRIAD satellite. If. however, we choose a solar- impact 

trajectory, then by using a spinning spacecraft, a one-axis drag-free 
system can be implemented that requires much less complexity. In 
fact, a spacecraft spinning about an axis aligned with the rectilinear 
path would allow 1) gyrostablllzatlon (from 0.3 a.u. to the sun), 2) 
an easier design for thermal shielding, and 3) a one-degree-of-freedom 
drag-free system. In particular, a sphere with an electrostatic 
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suspension is free to move along the spin axis with no exchange 
of forces along the path. The displacement of the sphere along 
the spin axis will be sensed, causing the thruster (oriented along 
the spin axis) to compensate the nongravitational forces along the 
path to the desired accuracy. The spacecraft will be forced to 
follow the proof mass and, therefore, to follow a purely gravita- 
tional path. Transverse forces should be'^ orders of magnitude 
smaller and need not be compensated. The drag-free system can be 
calibrated when the probe is far from the sun (5 a.e.) in order to 
find the equilibrium position along the spin axis of the proof pass 
in the gravity field of the spacecraft. 

During the3.5days that the spacecraft will spend in going from 
0.25 a.u. to 0.01 a.u. and closer, the earth-sun- probe geometry will 
pe>”mit the earth to be In the beam of the 0.2-m-diameter antenna 
pointing parallel to the spacecraft spin axis. The Jupiter-swingby 
technique has enough flexibility to enable the mission to be timed 
so that the earth-spacecraft line remains within a few degrees of 
the direction of the spacecraft track. A 20-cm dish mounted on the 
spacecraft operating in the X band has a beamwidth of 8*^ and a gain 
of 27 db. By using a 61-m dish on the earth, this will allow a 
transmission data rata from 50 to 100 bits/sec, even viith the noise 
of the sun in the background. Doppler tracking using two frequencies 
in the X band (8 and 12 GHz) should yield a relative- velocity measure- 
mer.t accuracy of the order of 10 cm/sec (1 a ) with 60-sec integra- 
tion time. 
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From oupiter inward to O.S a.u. • tiie spacecraft will operate 
in the Pioneer mode* performing selected experiments related to the 
solar stereoscopic and 00c missions. The two-frequency radio- 
science experiment will allow the integrated electron content to be 
determined at each instant, and perhaps a weighted component of the 
integrated magnetic field (frc»n Faraday rotation). 

From 0.5 a.u. inward, however, the mission will become more sun- 
oriented. The spacecraft spin axis will be directed sunward at 
this time, and the drag-free servo system will be activated. Sw^ 
relevant probe parameters inside 0.5 a.u. are given in Table 1. 




Table 1. Probe parameters. 


Solar Distance 

Time to 4R 

Velocity 

Equilibrium 

(a.u.) 

(R ) 

(days) 

(km/sec) 

Temperature 

CK] 

0.5 

100 

10.5 

54 

215 

0.4 

80 

7.4 

60 

235 

0.3 

60 

4.8 

70 

275 

0.2 

40 

2.6 

85 

335 

0.1 

20 

0.9 

121 

470 

0.05 

10 

0.25 

171 

680 

0.02 

4 

0 

270 

1060 
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The ten^eratures shown In the table are surface temperatures 
related to a properly designed reflective heat shield that covers 
the "front** side of the spacecraft; the Internal temperatures of 
the spacecraft will not necessarily be so high. Furthermore, the 
last few solar radii are traversed in less than an hour, during 
which time, theroml equilibrium will not be established. It is 
entirely possible that the system will survive to 2 Rgi If it does, 

a straightforward calculation shows that tracking to a doppler 

.2 

accuracy of 10 cm/sec over the half-hour interval required for 
the spacecraft to fall from 3 to 2Rg would permit Jg to be deter- 
rained to an accuracy of 10 . Since a realistic estimate of the mag- 
nitude of J 2 is about 10"^, as is shown in the Appendix, an extremely 
valuable result would be guaranteed. 

Obviously, many details of engineering design and scientific 
applicability remain to be worked out for this mission. But the 
preliminary effort so far expended appears more than sufficient to 
warrant this further pursuit. 
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APPENDIX 


J 2 OF THE SUN 

The value of 0 ^ can be inferred for a rotating axially sym- 
metric body by means of the first-order formula (Jeffries, 1970) 



( 1 ) 


2 

where f is the flattening and m = “ Rg/9g is the ratio of centri- 
fugal force at the equator to gravity at the equator. The assump- 
tions made are that the gravitational potential is given by its 
first two terms only. 


V 


9 





(sin «) , 


( 2 ) 


and that the surface is rotating uniformly so that the centrifugal 
force can be derived from the potential 


V ^ 
c 2 


2 2 2 
( 1 ) r cos ^ 


(3) 


Then the actual surface will be a level surface of the potential 

V = V„ + V^. 

9 c 

if the surface is not rotating uniformly, we can modify equation 
(1) oy assuming that, at any latitude, the surface will be perpen- 
dicular to the resultant of the gravity force given by the gradient 
of equation (2) and the centrifugal force equal to u (<>)r cos<> and 
directed away from the axis of rotation. Assuming that the shape of 
tne surface is'given by r = (1-Y), we find, to first order. 
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Rc w^(^) 

(3Jp + -= )sin ^cos 


( 4 ) 


We assume that the angular velocity of the sun's surface can be 


approximated by w(^) ® - fc >2 sin^^ and obtain 


? Ja • f - i- ”o 


3 ^ 


(5) 


-5 


• (0 “ 
• 0 


where * <0^ R Fur the sun, we have m = 2.14 x 10 
14.4 Vday, and 4.5Vday» so the second term in equation (5) 
lies between 7.7 x 10”^ and 10.7 x 10”®, depending on whether the 
differential rotation is included or not. The best determination 
of the flattening of the sun (Hill, 1974) is f » (9.6 ± 6.5) x 10”®. 
It is clear that O 2 c®nnot be derived with any accuracy from equation 
(5), since it is the difference between two not very well-known quan- 
tities of nearly equal magnitude. 

If the sun is rotating uniformly and if the density distribution 
is known, J 2 can be directly calculated. Following an analysis by 
Sterne (1939a), we define the "apsidal motion coefficient," 


3 - n 


k = 


4 + Zt]' 


where is the value at the surface of the variable n, which is zero 
at r = 0 and which satisfies Radeau's equation: 


r ^ - n ” 6 + ^ (n <• 1) = 0 . (6) 

In equation (6), p is the density at r and is the mean density 
interior to r. Then, 
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Jg “ I km , (7) 

where m has been previously defined. The coefficient k depends 

solely on the distribution of mass within the star» ranging from 

zero for a completely concentrated star to 3/4 for a homogeneous 

star. Values of k have been calculated (Motz, 1952} for solar 

models by Schwarzschlld (1946) and by Epstein (1951). Motz obtained 

k » 0.00585 and 0.00599, which leads to Jg * 8.3 x 10“® and 8.5 x 10"®, 

respectively. Calculating J2 for three later solar models, we found 

J2 * 1.56 X 10”^ for a zero-age sun (Schwarzschlld, 1958)and J2 “ 

1.41 X 10"' and 1.20 x 10*' for two models of the present sun 

(Weymann, 1957, and Sears, 1964). Although we do not at present have 

detailed calculations of later solar models, we note that a recent one 

(Hoyle, IS 75), proposed to explain the low neutrino emission from the 

3 

sun, has the unusually low central density of 75 g/cm . The ratio of 
central to mean density is then 53.2, which is quite close to 54.2, the 
ratio of central to mean density of the "standard model," a polytrope 
of index 3. Russell (1928) found k = 0.0144 for a polytrope of index 
3, so we consider 2 x lO”^ to be a reasonable upper limit to the value 
of 02 for a uniformly rotating sun. 

It is of interest to consider a lower bound to k and, hence, to J2. 
The most concentrated star with a given central density is the general- 
ized Roche model, which consists of a homogeneous core, with a density 
equal to the central density containing all the star's mass, and an 
envelope with Infinitesimal 
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density. Radeau's equation can then be solved analytically (Sterne, 
1939b) to obtain 



With current estimates of the central density of the sun ranging from 
about 75 to about 150 g/cm^, we find the lower limit of J 2 to be 
between 1.4 x 10“® and 4.5 x 10“®. 
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Recently solar physicists have gained the ability to view the solar corona from 
a new perspective, by means of space observations at x-ray and extreme ultraviolet 
(EUV) wavelengths, which permits us to observe the forms of the hot, but very tenuous, 
corona against the relatively cool solar disk (Figure 1). Ground-based observaUuns 
of the corona, except for relatively low-resolution radio data, require natural 
(via eclipse) or artificial occultation of the bri^t disk, limiting coroiml observations 
to *'side views" throu^ the atmosphere exteixiing beyond the occulting disk. The direct 
face-on view provided by EUV and x-ray space observations, combined with detailed 
information on density and temperature of the emitting regions conveyed by the spectral 
character of the ultraviolet and x-ray data, has given us the first detailed information 
on the complex structures of the corona and their interrelations. 

The possibility of looking down on the sim from a spacecraft hi^ above the 
ecliptic plane opens vq> yet another perspective on the sun, one that is certain to help 
us understand the nature of coronal holes at high latitudes, and their relationship 
to the expansion of the three-dimensional solar wind. 

The large dark area on the x-ray image of Figure 1 is a coronal hole — a pheno- 
menon whose properties are just now becoming understood, largely through EUV and 
x-ray observations. Although the existence of large regions of low coronal density 
had earlier been inferred from coronajraph "side views", we are now obtaining detailed 
knowledge of their size, shape, density, temperature, magnetic field, and evolution 
from spacecraft observations at EUV and x-ray wavelengths. 
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Briefly, cormal holes are large regions of the corona v^iose density is some 

3-10 times less than that of the "averse” quiet corona. The temperature, at least 

at the level where mMsurable x-ray and EUV radiation is emitted, is also less, 

6 6 

having ^ value ai abmit 1 x 10 % instead of about 1. 8 x 10 *K, as in tiie aven^ 
corwa. Analysis of low-res<dution OSO data (Munro and Withbros, 1972) suggested 
that the temperature gradient b^ween the 10^ chromosphere and 10^*K conma is 
abont an order of magnitude less in coronal holes than in the average sun. As a result, 
the energy loss to die coromi by thermal conduction back to the chromosphere, which 
is a very important energy Sink for the average corona, may be much less significant 
for coromil holes. Due to the low density in holes the radiative losses are also less 
(which is of course why they appear dark in the x-ray image of Figure 1). Since 
energy losses by radiaticm and conduction are both drastically decreased in holes, 
we may conclude that either the heating of the corona is less in holes or some additional 
mechanism of energy loss is present in holes that is not found in other regions of the 
quiet sun. 

Coronal holes are assuming great importance today because of their apparent 
associatioa with hi^-speed streams of the solar wind. The association, first suggested 
by Krieger et al (1973) on the tetsis of rocket x-ray data, has been put on a much firmer 
footing from detailed correlative studies of OSO-7 and SIqrlab data during die period 
1973 to 1974. Nolte e'*' al (1976) found that every large near-equatorial hole observed 
during the Sl^lab mission was associated with a hi^-speed stream at 1 a.u. Furthermore, 
there was a clear positive coirrelation between the velocity of the observed solar wind 
stream and the area of the associated coronal hole as measured from Skylab x-ray 
photographs. Finally, a very high correlation was fouiio between the polarity of the 
interplanetary magnetic field associated with the high-speed streams and the magnetic 
field underlying the associated coronal holes. 

During the Skylab mission coordinated ground-based and space data revealed 
that the location of coronal holes can also be detected through ground-based observations 
of subtle properties of certain Franhofer absorption lines, notably He 1 ,1 10830 
(Harvey et al, 1975). This discovery has permitted the mapping of coronal hole 
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boundaries during the two years after the terminatioa of the Skylab mission. These 
tu*o years were a time of large and persistently recurring bis^-speed solar wind 
streams, and significantly, a pattern of l<mg-lived coronal holes was detected from 
the He 1 10330 observations. Harvey et al (1976) lutve shown that, as before, there 
is a good correlaticHi between central meridian passage of holes that cross the equator 
and hi^-speed streams during diat time. In addition, when the data also include 
high latitude holes that also e.'ctend down to ± 40* latitude, the correlation becomes 
even better. Thus there appears little question that the high-speed streams are 
related to coronal holes, and in additicm there is evidence that some of the streams 
observed in the ecliptic plane are associated with holes at latitudes at least as high 
as 40*. 


To establish that the holes are the origin of the high-speed solar wind streams, 
however, it is necessary to identify a physical mechanism in addition to finding a high 
correlation between the two phenomena. In a search for a physical mechanism, we 
first note that coronal holes occur over areas oS unipolar photospheric magnetic fields. 
Although magnetic fields are measurable only in the photosphere, they may be ma{^d 
upwards into the corona, using potential theory al(»ig with the assumption that above 
about 1. 6 to 2. 5 solar radii they are stretched out radially by the expanding solar 
wind (Newkirk, 1972). Such calculations show the magnetic fields underl 3 ring holes 
to reach the source surface and thus to open out into the interplanetary medium, while 
fields underl}ring other regions in the corona generally close back on themselves (Altshuler 
et al, 1976). Thus holes seem an easy pathway for the escape of coronal plasma 
into the solar wind. 

We have already noted that the density and temperature structure of holes in 
the low corona suggest that either the coronal heating rate is less in holes or excess 
energy may be available to accelerate the solar wind outward in holes. It is interesting 
that coronal holes are almost impossible to detect in the chromosphere or below, 
suggesting that at those levels the atmospheric structure does not depend on whether 
or not a hole e.xists in overlying corona. This suggests (although it does not prove) that the 
amount of mechanical heating that passes upward through the photosphere is Independent of 
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the existence of coronal holes. Under that assumption Pneuman (1973) 

and Noel (1973) showed that the solar wind, expanding outward in regions of open 

field lines, would carry off energy thieugh acceleration of the solar alnd and throu^ out\^’ard 

thermal conduction, sufficient to comnensate for the decreased energy losses from 

holes by radiation and inward thermsi conduction. The exact mechanism by which 

originally closed field lines break open to allow fiie expansion of tte solar wind and 

the creation of a hole is not yet clear, but arguments based on energy flow support 

the reality of the process. 

If radial outflow really occurs o\'er coronal holes, it should give rise to 
observable doppler shifts of XUV emission lines. Preliminary reports (Cushman and Rense, 1976) 
indicate the detection of outward velocities of the order of 16 km/sec, which may in 
fact be the beginnings of the solar wind expansion. However, the data are scanty 
and further verification is needed. 

What does all of this have to do with the out-of-ecliptic mission? The significance 
lies in the fact that coronal holes ha ve been found to occur very frequently at the solar 
poles. The polar holes appear to be very similar to equatorial holes in their physical 
properties, with a major difference that they are much larger. A plausible hypothesis 
is that they too give rise to high-speed solar wind streams. These streams may 
emanate from the poles but spread out to lower latitudes, even reaching the ecliptic 
plane if the polar hole extends to low enough latitudes. 

Unfortxmately observations of polar holes from the orbit of earth still suffer 
partly from the projection problem described at the beginning of this paper. Because 
polar holes always occur near the limb (as opposed to near-equatorial holes, which 
are carried past disk center by rotation), we always observe them from the side. 

This of course leads to loss of spatial resolution due to foreshortening. In addition, 
observations from the side are particularly troublesome for observations of coronal 
holes, which are by their nature only very weakly emitting, and therefore are very 
easily obscured by foreground and background emission from the neighboring "normal'' 
corona. 
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It is also extremely difficult to get accurate measurements of magnetic fields 
associated with polar holes. Firstly, projection effects lead to loss of resoluti<m. 
Secondly, the sensitivity of a magnetograph is proportional to the line-of-si^ 
component of the field direction; for vertical fields near the limb this becomes very small. 

Finally, of course, velocities radial outflow above a coronal hole at the 
polar limb would not give rise to a line-of-sight doppler shift when observed from 
the orbit of earth. 

What might one hope to observe from a spacecraft situated over the pole ? 

Figure 2 shows a reconstruction of the appearance of the south polar hole as it would 
have been observed from such a vantage point during nine months of the Slgrlab mission. 
The images of Figure 2 were rectified using Sl^lab data from the Naval Research 
Laboratory XU V monitor instrument (Sheeley, 1975, personal communication). The large 
size of the hole compared to typical equatorial holes (cf. Figure 1), and its extended 
lifetime, are immediately apparent. 

A spacecraft able to observe the sim from higher latitudes (say greater than 60') 
for several months at a time, and properly instrumented, should be able to accomplish 
many significant observations of coronal holes. Large and sophisticated instruments 
such as have been flown on Skylab are by no means necessary, and are probably out 
of the question for the foreseeable future. The following are examples of important 
observational objectives, that could be met by realistic instrumentation aboard an 
out-of-ecliptic mission. 

1) Continuous mapping of the location of polar holes, and study of their 
evolution. A simple imager at any of a number of XUV or x-ray wavelengths, chosen 
such that the emission within the band pass largely originates at temperature in 
excess of about 1.5 x 10 °K, would be adequate. Spatial resolution of about 30 arcsec 
would be sufficient. For a small instrument, count rates would be quite low, but time 
resolution need be only of the order of many hours, so long integration times are 
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possible. ImagBs could bo built up by scazmbig a point detector (perhaps using the 
rotation the spacecraft for scannins^. Data on the location of fiie hole and its 
boundaries would be correlated with measurements from tte same spacecraft of 
local plasma parameters (density, velocity, magnetic field, composition, temperature), 
and as the spacecraft traverses directly above different parts of the hole or its 
boundaries, some idea of the three-dimensional flow field could be obtained. In 
addition, from the absolute intensities recorded, some useful limits on the density 
and temperature of the emitting plasma inside die hole could be obtained. 

2) With imagery at two or more XUV or x-rav wavelengths, one can obtain 
much better information on the idivsieal conditions in the coronal hole itself. 

Approximate values of density and temperature can be determined independently, 

and combined with modeling techniques, the data can give information on the variation 
of the i»rameters with height. Since both the density and temperature distribution 
in the low corona stron^y determine the plasma flow properties at one a. u. , correlation 
with these properties measured at the spacecraft itself will be very important. 

3) XUV spectroscopy at high spectral resolution ( x 10 ) would be 

very useful to measure the outflow velocity of material in the polar hole, in the manner already 
reported by Cushman and Rense (1976) for equatorial holes. Unfortunately the weakness of 
XUV emission lines in coronal holes, combined with the requirement for high spectral 
resolution, implies either a rather large instrument or extremely long integration times. 

This experiment, while very important, may therefore not be a suitable candidate fora very 
early exploratory out-of-the-ecliptic mission. 

4) Measurement of the polar magnetic field from an out-of-ecliptic spacecraft 
appears to he a natural and important objective. As mentioned aboTO, there are 
considerable ach^antages in obser\ing polar magnetic fields from mote nearly above 
the poles. Spatial resolution of 30 arcsec would be adequate to detei mine the gross 
structure of the fields and to follow their evolution. A sjnall ria‘;neiograpb operating 
with a solid ctalon fabry-perot filter in visual wavoloagtlis might well \ e feasible 
for inclusion on an exploratory out-of-ecliptic mission. 
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All of the above objectives would be considerably furthered by simultaneous 
measures from the ecliptic plane, in order to obtain stereoscopic information. 

In the case of ma^etic fields, for instance, observations from the earth would record 
those fields not recorded from the out-of*«cliptic spacecraft, and \dce versa. 

Comparison of relative signal strengths from such paired observations could help 
determine the vector field in the photosphere, thus putting potential mapping of high- 
latitude magnetic fields on a more secure footing. Similarly observations from an 
earth-orbiting satellite of the XUV or x-ray structures in coordination with simultaneous 
out-of-ecliptic observations of the same structures would yield the 3-dimensional 
structure unambiguously. 

We should note that the success of a program to study polar holes from out-of- 
ecliptic may depend strongly on the phase of the cycle during which the program is 
carried out. It appears that polar holes may shrink and even disappear near simspot 
maximum, doubtless reflecting the shrinkage and disappearance of the unipolar magnetic 
field cap associated with reversal of the general dipole field of the sun, which occurs 
about that time. At about sunspot minimum, unipolar magnetic fields at the 
poles and polar coronal holes appear to reach their greatest extent. 

From the point of view of studying polar holes and their relation to the solar wind, then, 
it may be useful to time the passage of an out-of-ecliptic satellite over the solar poles, 
or at least over hi^ latitudes, to occur with a few years either side of sunspot minimum. 
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Figures 

Figure 1. Skylab soft x-ray image of the sun, June 1, 1973. Filter bandpass 
2-32A and 44-54A. See Vaiana (1976) for details. Photo courtesy 
American Science and Engineering, Inc. , and Harvard College Observatory. 

Figure 2. Diagram of the evolution of a coronal hole as it would have been seen from 
above the south pole. Data are rectified from ATM Skylab observations, 
May 1973 to February 1974, made with the Naval Research Laboratory 
XU V monitor instrument. Data courtesy N.R. Sheeley, Naval Research 
laboratory. 
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The state of the solar atmosphere, which is considered to Include 
the solar wind, Is largely determined by magnetic fields. Both the magnitude 


and the configuration of fields at photospheric level appear to determine 
the flux of non-radlatlve energy, mass, and momentum Into the base of 
the corona. Likewise the interaction between coronal field and pli’.sr.a 
modulates the flow and produces the density distribution we see In the 
corona as well as the state of the Interplanetary medium. Rapid annihila- 
tions of the field In the lower atmosphere associated with flares and 
slower readjustments of the field associated with eruptive prominences 
both lead to disruptions In the structure of the overlying corona and 
solar wind. 

Although a definitive test of the relation between coronal morphology 
and magnetic fields Is still lacking, several characteristic coronal 
structures are associated with distinct topologies of the field: 

. The tightly closed coronal loops and arches above active 
regions with similar configurations In the field. 

. Coronal holes at low latitude or over the poles with 
magnetic fields which open directly into the solar 
wind. Such regions appear to be the origin of fast 
streams In interplanetary space. 

. Coronal streamers with large arcades in the magneti'' field 
below 2 to 2.5Rg and a current sheet above that height. 

Such features are detected as the "sector" boundary between 
large scale, oppositely directed fields in interplanetary 
space. 

. Hot, dense knots of plasma in the very low corona and visible 
as "bright points" in X-ray images with minute bipolar regions 
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distributed more or less uniformly over the surface. In the 
upper corona and Interplanetary space such regions are hypothe> 
sized to give rise to minute current sheets which play a role 
In determlninq the electrical and thermal conductivity, the 
propagation of radio waves and energetic particles, and. 
possibly, coronal heating. 

The characteristics of the photospheric field In magnitude, spaclal 
extent, and lifetime suggest that the corona and Interplanetary medium 
can be divided Into 3 regions having relatively distinct properties as 
shown In the table. We note that although Interplanetary neasurements 
have sampled both the plasma originating in the latitude zone 5 10® and 
shocks originating at h gher latitudes, the influence of the 10® - 50® 
zone, which contains active regions, cn interplanetary medium is 
still uncertain. 

A real understanding of the structure and evolution of the corona 
and interplanetary medium can be claimed only after we have constructed 
a self consistent 3-n model of the entire region and have tested it with 
concomitant observations in the lower solar atmosphere, in the corona, 
and in interplanetary space. The Out-of-the-Ecliptic Mission will provide 
not only the critical tests of such models but will also afford insight 
into the fundamental mechanisms governing this entire region through 
the sampling of zones where different field topologies, magnitudes, and 
evolutionary timescales (and presumably different mechanisms of mass, 
energy and momentum transport) dominate. Several specific questions may 
be considered: 
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Are the polar regions Identical to the coronal holes at 
low latitude? 

. If this Is so, how do we account for the fact that polar 
regions with apparently Identical magnetic configuration 
display a vastly different appearance In the corona? 

. Does the coronal and Interplanetary microstructure change 
with latitude? What Influence does this have upon the 
heating, electrical and thermal conductivity, wave content, 
and energetic particle propagation of the medium? 

. How is the corora "mapped" into interplanetary space? 

. How does the solar-interplanetary field couple to the 
interstellar fielc? 

. What role does the continuous occurence of coronal transients, 
which originate at t< 60°, play in determining the state of 
the interplanetary medium and energetic particle propagation? 

. What mechanisms control coronal and interplanetary abundances? 

These are but a few of the questions concerning the rcle of solar magnetic 
fields in determining the structure of the corona and interplanetary medium 
which will be explored by an Out-of-Ecliptic Mission. 

Concerning the mission options two points should be kept in mind: 

1) Although the appearance of the corona suggests that a mission 
restricted to £<40° might be successful in reaching the polar 
zone, only a polar mission can guarantee that a truly new region 
of the interplanetary medium is to be explored. 
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2) The corona and interplanetary medium are continually evolving. 

To assure success the Out-of-Ecliptic Mission must enconq)ass a 
coordinated program of solar and interplanetary measurements 
so that a coherent attack can be made upon the important problems 
which beckon. 

This presentation was made possible through the generosity of colleagues 
at Kitt Peak National Observatory, High Altitude Observatory, Smithsonian 
Astrophysical Observatory, and American Science and Engineering who kindly 
allowed the use of illustrative material. 
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3*D Solar Radioes tronomy and the Structure of the Corona and the Solar Wind 
J.L. Steinberg and C. Caroubalos 
Obaervatoire de Paris ,Meudon, 92 190, France. 

Solar radio bursts are intense radio emissions from localized 
regions in the corona and interplanetary medium. Their brightness 
temperature is so much higher than the electron temperature of the 
ambient plasma that the mechanism which produces them is certainly 
non-thermal ; the necessary energy is brought into the source region 
by energetic electrons which can give rise to different types of 
transient radiation. One of them is the type III radioburst produced 
by energetic electrons travelling along open magnetic lines of force. 

Another quite usual form is the type I burst, typical of the meter wave 
lengths range. 

We do not fully understand yet how the non-thermal energy is 
converted into electromagnetic energy ; but we know that, to be efficient 
enough, the conversion must take place at frequencies close to the 
resonant frequencies of the medium. At these frequencies, the refractive 
index for radio waves takes extreme values : close to 0 for the plasma 
resonance, much larger than I for the gyro-resonances. As electromagnetic 
waves travel away from their source, those resonance conditions are no 
longer fulfilled because of the non-uniformity of the electron density 
or magnetic field ; the refractive index comes back quickly to unity ; 
such a variation of the refractive index is favourable to beaming effects. 
The radiation mechanisr itself which can involve some amplification may 
also produce a directive primary emission. In both cases, the beam will 
be oriented along or at definite angles to the principal directions of the 
medium : the electron density gradient or the magnetic field. 
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Using results from the STEBEO-I e^erimentt it is shown 
that stereoscopic observations in the deci to decameter X range 
can provide information on 

- the burst emission mechanisms 

- the local electron density gradient and magnetic field 
vector at the source 

- the macrostructure of the corona and the solar wind 

- the characteristics of small scale electron density 
inhomogene i t ie s . 

Radio bursts of type III can also be used to map solar 
magnetic field lines of force throughout the iuterplanetarv 
mediiim up to the earth orbit and beyond. 

Future experiments of these kinds should be carried out 
between an out-of-the-ecliptic probe and the earth or an earth 
satellite. 
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Beaming of the radiation of type 1 and type III bursts was 
predicted long ago» but not observed until recently. Rather than 
giving up the basic mechanisms which seemed capable of explaining 
many observed properties but implied beaming effects, radio astronomers 
suggested that random inhomogeneities of the refractive index scattered 
the radiation. (Roberts, 1959). This suggestion received strong support 
when scintillations from radio stars seen through the upper corona 
were discovered. Fokker (1965), Steinberg et al (1971), and Riddle (1972) 
carried out Monte Carlo numerical computations of the scattering of the 
radiation from a source embedded in the inhomogeneous medium. These 
authors used models which were extrapolated to low coronal altitudes 
from measurements of radio scintillations made for paths which did not 
cross the corona lower thar 5 solar radii or so. However, these studies 
were successful in accounting for several observations which could hardly 
be explained in any other way. They showed that the inhomogeneous 
medium produces a scattered image broader than tl- source and appreciably 
displaced from it ; at the same time, the random propagation tends to 
suppress any beaming of the radio waves and smoothes the radiation pattern 
of the source. 

Measurements of the angular distribution of the intensity of a 
source of radio bursts can therefore yield information on : 

- the orientation of the principal directions of the medium : 
grad Ng or the radiation mechanism and beaming processes. 

~ the characteristics of small scale inhomogencities which 
cannot be obtained in any other way. 
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To measure lUrectivity tobservations from the ground only are 
inadequate. For many years » authors tried to reach at least a 
statistical view of the directivity of radio bursts ; for instance, 
from their E.W. probability of occurence ; but, if the orientation 
of the radiation pattern of individual bursts relative to the local 
vertical through the source is not constant, no information on the 
directivity can be obtained from a single observing site. Simultaneous 
observations should be made in at least two widely different directions 
(Steinberg and Caroubalos, 1970). 

We have seen that the suprathermal electrons which produce type 
III bursts are guided along open magnetic lines of force. These lines 
are carried av-ay by the solar wind into interplanetary space so that 
type III are observed from low in the corona to the earth orbit and 
beyond. If we were able to map the successive positions of the type III 
source, we could also draw 3-D maps of some solar magnetic lines of force. 

Stereoscopic observations of radio bursts are powerful tools to 
study the corona and solar wind. This may be illustrated by some recent 
results obtained with the STT:RE0-1 experiment carried out in 1971-1972 
at 169 MHz; in cooperation between France and the Soviet Union. 

(Caroubalos and Steinberg, 1974; Caroubalos, Poquerusse and Steinberg, 1974 
Steinberg, Caroubalos and Bougeret, 1974). At 169 >31z, radio bursts of 
types I and III occur at altitudes in the range 0,3 “ 0,5 solar radius. 

STEREO RADIOASTRONOMY OF TYPE I BURSTS 

Let H be the stereo angle between the tv;o observing directions. 

When 6 increases, the corrolaticn between the two intensity- VS tine- 
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records taken simultaneously (in the source time scale) decreases in gene 
ral. Even with 0 - 15® » this is clearly visible, but when 0 - 35®, the 
correlation coefficient is less than O.I. This means that the beamwidth 
of type I radiation is sometimes smaller than about 25®. However, on 
some consecutive days, the same intensity may be received at both 
observing sites and then the beam pattern looks nearly isotropic. 

Such an apparent contradiction can be resolved if ue note that the 
STEREO- 1 observations were carried out in the ecliptic ; so that we 
are actually analyzing only a plane section of a 3-D beam pattern 
and we do not know the configuration of the beam pattern out of that 
plane ; we cannot, for instance, know if the 3-D beam pattern is solid, 
multilobed or even hollow. It is easy to conceive beam shapes whose 
cross sections by different planes can be either narrow or broad. 

The ratio R of the burst intensity measured in Space Ig to that 

measured at the earth varies widely from event to event ; so that the 
beam has to be randomly oriented if its shape is assumed almost constant 

the rms deviation of the orientation is about 0.25 of the beamwidth and 

this is a rather clear indication that the source does not contain a 

large number of inhomogeneities. 

This is, in turn, connected to an old problem : type 1 burst 
intensity can vary by a large factor in 0.1 second ; but the observed 
source size is about 3 arc min or 0.3 light-second ; so that it was 
suggested long ago (Hogbom, 1960 ; Fokker, I960) that what we see 
is actually the scattered image of a deeper and smaller source. In a 
scattering corona, the assumed small scale inhomogenci r ies do produce 
a broad scattered image of a point source but at the same time they 
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broaden the angular distribution of the radiation from that point source. 
Both effects are intimately connected together via the scattering power 
distribution along the path. The total rms random angular deviation of 
the radio rays over, their trip from the source to us cannot be larger 
than half the observed beamwidth ; STEREO observations yield directly 
a measure of the beamwidth and»thns» an upper limit to the rms angular 
deviation along the path* This limit is too small for the existing 
models to account for more than a small part of the image size. Therefore 
either there is less scattering than generally assumed to account for 
mos t of the source apparent size or the inhotnogeneities built in the 
models are inadequate. 

In any case it has been demonstrated that very interesting 
information on the beam orientation and shape can be obtained from 
stereoscopic observations. To -learn more, it is necessary to go out 
of the ecliptic plane for the following reasoni; ; 

- to compare the beam orientation to that of the density gradient 
we must know the latter and therefore the 3-D electron density 
distribution in the source region. This can be obtained from 
coronagraphic measurements on the limb where the electron density 
distribution as a function of latitude will always be better knouTi 
than the longitudinal one. It is therefore much more effective to 
measure the beam orientation in a plane perpendicular to the ecliptic 
than in the ecliptic. 

- operating a stereoscopic cxperirient between nn out-oi “’the- 
ecliptic probe and the earth will also provide a larg(*r variety of 
cross sections of the beam pattern using solar rotation (fig, 1); in 
the ecliptic, solar rotation moves the same cross sect ion of the 
beam across our lines of sight. Using an out-of--the-ccliptic set* up 
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one should be able to get much closer to a conqplete descriptiun of the 
3-D beam pattern of the bursts. 

If the out-of-the-ecliptic probe is on a 1 AU orbit, in a plane 
tilted to the ecliptic, the stereo angle will also vary quite rapidly 
and this is again favourable to a detailed description of this beam 
pattern. 

STEREO RADIOASTRONOMY OF TYPE III BURSTS 

The spectrum of a type III in a frequency-time domain (dynamic 
spectrum) shows a band of noise drifting from high to low frequencies. 
This band is sometimes split in two components which, at a given time, 
are centered on harmonic frequencies. Some type Ill’s are therefore 
made of two components which ar;^ believed to be produced, one (the 
"fundamental") at the local plasma frequency f^, the other at twice that 
frequency. When observed with a single frequency receiver, the first 
con^>onent is recorded first and the second some seconds later, making 
up a "pair" of type Ill’s. 

This interpretation of pairs as fundamental-harmonic pairs has 
been questioned in recent times but not in a convincing way ; and to 
settle that question, directivity measurements are important : indeed, 
the conversion mechanism and the propagation conditions are different 
for the fundamental and the harmonic ; for instance, the fundamental 
is generated at about the local plasma frequency so that it should be 
beamed into a narrow cone ; and if the corona is assumed quasi spherical 
this cone should be about radially oriented so that few fundamental 
components should be seen in high longitude events ; this is not the case 
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fundamsntal components ere seen nearly all over the disc. 

Stereoscopic observations showed that the first (fundamental) 
component of a pair is systematically more directive than the second 
(harmonic) and this is a strong argument for the fundamantal-harmonic 
interpretation of pairs (fig.2). 

The time profile of a type III at a fixed frequency is also 
rich in information as it is the convolution of an exciter function 
by the transient response of the corona which includes the effects 
of multipath propagation. At 169 MHs» it was found independent from 
the direction of observation ; therefore propagation conditions do not 
play an inq>ort-:jtt role in the formation of the time profile. 

Vniile the time profile is independent of the direction of observation 
the intensity ratio R ■ Ig / Ig can take values very different from unity 
the rate of change of the intensity with the observing directions can 
reach i 10 dB and more over 30°. This proves, again, that coronal 
scattering is less effective than previously thought ; even less effective 
than necessary to account for some other observations ; for instance, 
scattering- has been invoked to explain that, at a given observing 
frequency f, the fundamental component (local plasma frequency fp > f) 
and the harmonic component (2fp " f) are observed at the same position 
although the first should take place at the f critical level and the 
second at the f/2 critical level, higher up in the corona. This 
observation is indeed explainable in a scattering corona (Riddle, 1972 ; 
Leblanc, 1973) but the scattering power las tc be larg.er than the one 
deduced from directivity measurements. 

Another result from STIiKEO observations is that the directivity 
ratios of various type Ill's can be very different. This can be 

pi ; ■ muCIBILITY OF THE 
PAGE IS POOR 


- 72 - 



interpreted 1£ the observed directivity is produced by coronal macro-* 
structures mostly in the form of streamers. On Nov.l4| 1971 » for instance 
type III have been observed through a streamer ; the overdense streamer 
material reflects » absorbs and scatters the type III radiation, away 
from its source and produces the observed directivity. If this 
interpretation is correct, the streamers oi "lames coronates" detected 
by Axisa et al (1971) do control the type III image sise and shape as 
a piece of ground glass or a light shade* Observations of these images 
with radioheliographs together with Stereoscopic observations can be used 
to study the streamer structure which is hard to resolve optically because 
of line of sight integration effects. We still do not know where aru the 
type III sources located as compared to streamers ; to settle that question 
2-D position measurements at radio frequencies are necessary but at the 
present time, they are no more accurate than ! arc min or so. Occultation 
effects are only detectable with Stereoscopic observations but they are 
very sensitive to the position of the source relative to the occulting 
structure. They open up new ways to localize the path of the type III 
electrons relative to streamers and the site these electrons are 

accelerated in the active region. Here again che stereoscopic observations 
dould be carried out on an out-*of-the-ecliptic probe because the macro- 
structure of the corona is better known from optical observations as a 
function of latitude than of longitude. 

One of the most useful properties of type Ill’s is that they arc 
produced over trajectories which span the interplanetary ricdium. At each 
altitude the 10-100 keV electrons induce plasma waves which are scattered 
into electromagnetic W£vcs at the local plasma frequency f. or at twice 
that the frequency. Therefore from the measurement of the position of 
type Ill’s at several frequencies, a map of the electron density along 

the trajectory and a map of that trajectory itself can be drawn. 

11 



As the energetic electrons travel most probably along magnetic lines 
of force, we have a way of plotting such sun*rooted lines of force up 
to the earth orbit and beyond even out of the ecliptic. 

Such an erperiisent should be carried out at frequencies lower 
than 10 MBz if we are interested in the coronal and solar wind structure 
hi^er than 10 ; this means that observations should be made from 

space as radiation of these frequencies do not reach the earth. As a 
matter of fact has just done that (Fainberg and Stone, 1974). 

IMP-6 was spin stabilized around an axis perpendicular to the ecliptic 
and carried a dipole perpendicular to that axis. Using the nulls in 
the receiving pattern of a short electric dipole it is quite possible 
to measure the direction of a source as projected on the ecliptic plane. 

An experiment jointly designed and built by Paris Observatory and 
Goddard Space Flight Center teams will measure die direction of the type 
III source at 24 frequencies on ISEE-C. Using a spin plane and a spin 
axis dipole, the experiment will measure a conqilete direction (two angles) 
at each frequency aqd will produce 3-D maps of some magnetic lines of 
force from 10 altitude to the earth orbit and beyond. It will, 
however, be necessary to assure that these lines of force rotate with 
the sun as a solid body. The use of a second remote satellite equipped 
in T.ach the sane way as ISEE-C could eliminate this restriction, (fig. 3) . 

There are some indications from the radioastronomy experimen 9 
on lMP-6 and other experiinents that few type Til have been detected 
far ojt from the ecliptic. This m^ght very veil be due to son^ 
dirccLlvity of the radiation, but this difficulty can be overcome by 
going out of the • cliptic. 





Oae o£ the main purposes of any out*-o£'*the-ecliptic mission 
will certaimly be to e^lore the 3-D topology of the interplanetary 
magnetic field and more specifically its latitude variation ; the role 
of solar active regions in the determination of this topology will be 
studied, ii^quipiaents designed to measure the local magnetic field vector 
will be flown to achieve this goal but it will be very hard to reconstruct 
the magnetic configuration in- the whole heliosphere from local measurements 
mly. Type 111 tracking at several frequencies can provide the overall 
description of this field topology which will be essential to the 
interpretation of OK>st local measurements made on the 0/E probe : for 
instance the modulation of cosmic rays by the interplanetary ragnetic 
field cannot be understood without a description of this magnetic field 
in the whole heliosphere. 

CmCLUSlOti 

Stereoscopic observations of solar radio bursts are not needed 
only to inprove our knowledge of the physics of these transient radio 
emissions. In the deci- to decameter-X range, they can be used to 
probe the macro and microstructure of the corona. In the hm to 
km-X range, type Ill's are natural tracers of sun-rooted magnetic lines 
of force ; tracking them as a function of frequency will give a 3-D map 
of some lines of force from low in the corona to the earth orbit and 
provide an overall picture of the interplanetary medium which is essential 
lo the interpretation of local measurements. 

The choice between the two wavelength ranges depends upon the 
scientific object *ves of the mission ; but it has been shown that, in 
h< *s, observaticr -• should be carried out from an out-of-ecliptic 

probe. 
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Steinberg - Caroubalos 


3-D Solar radioaatronony 


Figure Captions. 

Fig. I Exploration by stereoscopy In a plane perpendicular to the 
ecliptic of the radiation pattern of a radio burst using solar 
rotation. The plane sections I to 3 are analyzed at different 
t'lnes. If the exploration was carried out in the ecliptic » 
only one plane section would be studied. 

Fig. 2 A typical pair of type III bursts at 169 (fi^z as recorded 

irom the earth (top) and from the Mars-3 Soviet s'lce probe 

on Nov. 14, 1971. The intensity ratio I_ /I„ , 

’ ^ Space Earth 

first component (ftmdamental) is always greater than that of 

the second coinponent (harmonic) . 

Fig. 3 3-D mapping of a solar magnetic line of force using a type 

III burst as a tracer. 
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Abstract 


Interplanetary scintillation (IPS) observations from 1971-1975 show 
that the average solar wind speed increases away from the solar equator, 
with a mean gradient of 2.1 km/s per degree. These results are compared 
with spacecraft observations over the t 7° attainable in the ecliptic and 
with those deduced from comet tails. The role of temporal variations, 
especially those caused by latitude dependent solar wind streams, is 
emphasized. This points to the need for extensive ecliptic and ground- 
based observations during an out-of-the-ecliptic spacecraft mission. 



Introduction 

The solar wind both In and out of the ecliptic can be studied from 
the earth by the method of Interplanetary Scintillations (IPS). The 
method was pioneered by Hewlsh and his colleagues, who deduced the solar 
wind speed from multiple station observations. In 1966 Dennison and 
Hewlsh (1967) found an Increased speed out of the ecliptic, while In 1967 
Hewlsh and Symonds (1969) found no such Increase. Their multiple station 
observations then ceased. In this paper we report solar wind speeds 
deduced from 74 MHz IPS observations made at DC San Diego from May 1971 
through to April 1975. The observing system was described by Armstrong 
and Coles (1972) and by Coles et al (1974). Results from 1972 were reported 
by Coles and Maagoe (1972). We will also mention briefly the relevance 
of IPS observations simultaneous with an out-of-the-ecl iptic spaceprobe 
mission. 

IPS Method 

The scintillation signal is the sum of waves scattered along the line 
of sight from a given radio source. Most of the scattering occurs where 
the line of sight is closest to the sun, because of the steep decrease with 
solar distance in the strength of the electron density microstructure which 
causes IPS. For a spherically symmetric solar wind a weighting function can 
be defined and the IPS "mid-point" speed (Coles and Maagoe 1972) can be 
shown to be a spatial average of the solar wind speed centered on the point 
of closest approach. However, in the presence of solar wind streams spherical 
symmetry does not apply. 

The effect of the spatial average through such streams iias been inves*igated 
by comparing the IPS observation with those expected by mapping point observa- 
tions made on the IHP-7 spacecraft out along the line of sight in question. 
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Harmon (1975) and Coles et a1. (1975) demonstrated a close agreement 
between IPS "mid-point" speed and the IMP-7 data mapped to the point where 
the line of sight is closest to the sun. More detailed comparisons are in 
progress, investigating the precise form of the spatial weighting caused 
by streams. However, for the present purposes the comparison demonstrates 
that each IPS observation is representative of the solar wind speed at the 
point of closest approach. 

Results 

This effective observing point changes in solar latitude, longitude 
and radial distance as the sun rotates and the earth orbits the sun. Thus 
for radiosources not in the ecliptic, about two months of high solar latitude 
data can be obtained each year. The geometry, however, is such that the 
high latitudes occur together with small solar distances. The distance 
dependence can be separated out by studying the four ecliptic radio-sources, 
for which the latitude remains within 10° of the equator; Figure 1 shows 
the solar wind speeds averaged into intervals o^" 0.1 AU in radial distance 
during 1971-1975. We conclude that there is no significant variation of 
average solar wind speed with radial distance between 0.4 AU. and 1.1 AU. 
Figure 2 shows a similar plot for all sources versus latit ' and, because 
there is no radial distance dependence, it can be interpreted as showing 
the solar wind speed as a function of latitude. The vertical bars are 
t 2 standard deviations in the average solar wind speed over latitude inter- 
vals indicated by the horizontal bars. The r.m.s. variation in a single 
speed observation is remarkably constant at about 120 km/s, showing no 
significant change with latitude. The vertical error bars are larger at 
higli latitudes because there are fewer data points at high latitudes. (In a 
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typical year* 300 observations from 0-10°N decreasing to 25 observations 
from 50°-70\) The major conclusion Is that there Is a systematic Increase 
of solar wind seed with latitudes both north and south of the solar equator. 

This Is evident In each year from 1971 through 1975 as well as In the grand 
average of all data. The average gradient Is close to 2.1 km/s per degree 
of latitude. However, the curves are not quite symcetrlcal but centered near 
10°N, giving an apparently steeper gradient in the south than In the north. 

’’'his asymmetry Is only marginally significant and we are still checking for 
second order systematic errors which could cause this. 

Discussion 

Our observed latitude gradient must be compared with other data. As 
already mentioned, the Cambridge IPS observations detected a latitude gradient 
In 1966 but not In 1967. Whereas this could conceivably be Influenced by 
solar cycle effects. It Is more likely that the small number of observations 
at high latitudes is responsible. Their measurements included only about 
30 days each year which corresponded to latitudes above 20°; the long term 
average behavior could well be masked by the day-to-day and month-to-month 
variability found in the solar wind speed. 

Spacecraft observations have been analyzed to look for effects due to 
the + 7° latitude range available in the ecliptic (e.g. Hundhausen et al. 1971). 
Smith and Rhodes (1974) and Rhodes and Smith (1975) deduced large apparent 
gradients (10-15 km/s per degree) by comparing solar wind speeds observed 
near Earth (Explorer 33,34,35) and at Mariner 5. They analyzed data from 
nearly six solar rotations over latitude differences from 0° to 6°. We 
suggest that a possible explanation of their large gradient over a few degrees 
of latitude comes from solar wind streams. From our IPS data it is clear that 
solar wind streams often exist for several solar rotations with steep latitude 
gradients near the equator (see for example the wide southerly stream in 
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Figure 13 of Coles et a1. 1974). It Is likely that such features contribute 
strongly to a six month average. Figure 2 shows that In 1972 the 
apparent gradient between t 5** was ^ 4 km/s per degree, while In 1973 the 
apparent gradient reversed to - 3 km/s per degree. The Influence of 
steep latitude gradients from specific recurrent streams was probably the 
cause for such large values. 

More difficult to reconcile are the results from comet tail observa- 
tions. Brandt et a1. (1975) have analyzed 678 comet observations spread 
over 75 years and conclude that the latitude gradient is - 0.9 + 0.7 km/s 
per degree north or south; that is not significantly different from zero. 

Their observations are concentrated in the range 0° to 50°N (as are the 
IPS observations) and are scattered fairly well through the phases of the 
solar cycle. It would be of interest to see if the high latitude data 
were uniformly distributed over the phases of the solar cycle. Except for 
solar cycle effects we cannot suggest any simple explanation for the 
discrepancy; though long term systematic or random temporal variations could 
be responsible. 

Conclusions 

We have presented strong evidence that during 1971 through 1975 the 
average solar wind speed increased out of the equatorial plane giving an 
average gradient of 2.1 km/s per degree of latitude either north or south. 

Our observations show that stream and also slower variations can obscure the 
average latitude behaviour in the solar wind (as they do also for average 
properties in the ecliptic). In planning out-of-the-ecliptic spacecraft 
missions such changes must be expected. Jupiter "swing-by" missions would 
give 12-18 i.ionths at more than 30° from the equator. During this time it will 
be important to maintain regular observations in the ecliptic in order to 
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disentangle temporal and latitude effects. A continued program of IPS 
observations throughout the period would cover a range of longitudes and 
latitudes and further help build a picture of the spatial and temporal 
structure in the solar wind. In addition to the average latitude behaviour 
such joint observations would allow the latitude structure of individual 
streams to be explored. 
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Figure 1. Solar wind speeds from 1971-1975 measured within 10° 
of the solar equator, averaged into 0.1 AU intervals 
of radial distance. The vertical error bars are ± 
twice the standard deviation in the mean. 

Figure 2. Solar wind speeds from 1971 through 1975 averaged 

into latitude intervals shown by the horizontal bars. 
Vertical bars are ± twice the standard deviation 
in the mean. The lower right graph is the overall 
average from the other five graphs. 
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Abstract - Analysis of the orientations of ionic comet tails 
gives no support for the suggestion that the radial solar wind speed 
is higher near the solar poles than near the equator. These results 
refer to a long-term, global flow pattern and do not refer to short- 
term variations. 
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Evidence from comets concerning the latitudinal variation of 
solar wind parameters has been discussed previously by Pflug (1966), 
Brandt (1967), and by Bertaux, Blamont and Festou (1973). In this 
short report, I summarize the evidence based on Ionic comet- tail 
orientations as recently analyzed by Brandt, Harrington, and Roosen 
(1975), and show the distribution of the sample In latitude, time, 
and phase of the solar cycle. 

The basic observation Is the position angle of the tall axis on 
the plane of the sky (Belton and Brandt 1966). The position angle 0 
is interpreted In terms of a taj.l vector T^ whose direction in space is 
determined by dynamical aberration, viz., 

T = w - V (1) 

where is the solar wind velocity vector and V is the vector velocity 

of the comet. The astrometric technique developed by Brandt, Roosen, 

and Harrington (1972) does not assume that the comet tail lies in the 

plane of the comet's orbit. Each observation determines essentially 

a half plane in velocity space. A preferred solar wind velocity vector 

(w , w., w^) is determined as the one which minimizes the sum of the 
r* G 0 

squares of the residuals between the computed and observed position angle 

At present, a sample of 678 observations are available and these 

are spread over approximately 75 years in time and between roughly 0,5 

to 1.5 a.u. in heliocentric distance. Ihe basic results arc a radial 

velocity, < w^ > 400 km/ sec, an azimuthal velocity < > « 6-7 km/ sec 

2 3^5 

(varying with solar latitude b and distance r as cos ’ (b|/r), and a 
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RMS dispersion of 3?7. An additional result from comet tall 
orientations, but not from the astrometric technique, is that 

> 225 km/ sec (Brandt and Helse 1970). These values are in good 
agreement with results from spacecraft and provide confirmation of 
the basic approach. 

This technique has been previously used to search for a 

meridional flow pattern in the solar wind. A value w «« 2.5 km/ sec 

m 

(at 0 = 45^, varying as sin 20) has been found in the sense of a flow 

diverging from the plane of the solar equator by Brandt, Harrington, 

and Roosen (1973). This result implies a radial variation in the 

-2 013 

equatorial density of the solar wind of N « r * . If this law 

held from the sun to earth, the density would be VL smaller than on 
spherically symmetric models. 

The basic technique can be used to search for a latitudinal 
variation of the radial solar wind speed by assuming that it varies as 

dw 

= ''o ^ df] I*’ I <2) 

where w^ and dw^/d|b| are constants to be determined. The results 
are given in Table 1. 
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Table 1. 


Solar-Wind Speeds With and Without a 
Latitudinal Variation in Radial Speed 

w or w 
r 0 

(km sec"^) 

dw /d |b 1 

-1 -1 

(km sec deg ) 

w 

® -1 
(km sec ) 

-1 

(km sec ) 

RMS Dispersion 
in (8 - ©g) 

402.5±11.9 


+2.6±1.2 

7. Oil. 8 

3?749 

418,5±27.3 

-(0.9*0. 7) 

+2.9±1.3 

5.3±2.2 

3?750 


The latitudinal variation found in our sample, if any, is in the 
sense of decreasing radial speed with increasing latitude. However, 
the error in fdw^/djbl] is almost as large as the value of -0*9 km sec ^ deg ^ 
found, and there is clearly no trend. In addition, the best solution as 
judged by RMS dispersion is still the solution with [dw^/dlb|l ^ 0. 

When an additicnal significant parameter is included in the model, the 
dispersion must decrease even if only marginally. The lack of a decrease 
is a definite flag that the additional parameter has no significance. 

The slight increase in RMS dispersion is simply due to round-off error. 

The errors in the components of the solar wind speed increase with a 
latitudinal variation included because then all components are functions 
of heliographic latitude and can be correlated. Possible correlations 
between components are calculated and are used to assign the probable 
errors. This is the explanation for the increase in errors while the 
RMS residuals remained essentially constant. 
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The negative result for slgnlflcrnt latitudinal variation In 
refers only to a long-term, global situation. It does not rule out shorter 
term results such as the one presented by Rlckctt (1975; preceeding paper). 
It does appear to imply that such short-term vr.rlatlons average out 
over the long term. 

The sample is concentrated in the range 0 ^ |b { ^ 50^ as shown 
in Figure la where the solid line represents the present sample of 
678 observations and the dashed line represents the sam» number 
distributed at random; only 60 observations lie in the range 
50° s |b| i 90°. If we plot the distribution against latitude 
instead of the absolute value of latitude as shown in Figure lb, we 
find a strong concentration of observations in the northern hemisphere. 
There is no obvious reason to expect an adverse effect from this 
observational bias. The sample by year of observation is shown in 
Figure Ic which reflects the irregular nature of cometary apparitions 
and the reduction of observations. Figure Id shows the distribution 
of the sample with phase in the solar cycle; the observations a’-e 
concentrated toward solar maximum. 

Because the astrometric technique can be applied only to fairly 
large groups of observations, results on short-term variations in the 
solar wind speed cannot bo obtained directly. Work on an indirect 
technique is currently in progress. 

Nevertheless, there is ample direct evidence for large, short- 
term variations in solar wind properties. The time required to 
establish a meaningful measurement of an average property at a 
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particular Ictitude Is prubably at least one solar rotation. This 
vas found to be the case for spacecraft observations of as reported 
by Lazarus and Golwsteln (1971). Hence, direct out>of-the-ecllptic 
observations of the solar wind s'tculd utilise an orbit with a slowly 
changing latitude. Several passes througli all solar latitudes and 
possibly several spacecraft will be required to map out the basic 
structure of the solar wind in three dimensions. 
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FIGURE CAPTION 


Figure 1. The distribution o£ the sample of comet observations 
(a) in absolute value of solar latitude; (b) in solar latitude; 

(c) in date of observation; and (d) in solar cycle phase. See text 
for discussion. 
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Los Alamos, Nev Mexico 875^5 

ABSTRACT 

A comparison of polar solar wind proton flux upper limits derived 
using Saito’s coronal density model, with Ly a meastirements of the 
length of the neutral H tail of comet Bennet at high latitudes, shows 
that either extended heating beyond 2 is necessary some of the time 
or that Saito's polar densities are too low. Whichever possibility is 
the case, the fact that the solar wind particle flux does not appear to 
decrease with increasing latitude, indicates that the heavy element 
content of the high latitude id.nd may be similar to that observed in 
the ecliptic. It is then shown that solar wind heavy ion observations 
at high latitudes allow a determination of the electron temperature at 
heights which bracket the nominal location of the coronal temperature 
maximum thus providing information concerning the magnitude and extent of 
mechanical dissipation in the intermediate corona. 



1, Introduction 

There is, at present, a lack of information on the physical 
conditions in the polar regions of the solar corona and solar wind. This 
lack results in a corresponding uncertainty in the global characteristics 
and esctent of that plasma of solar origin vhich fills interplanetary 
space and thereby controls the near solar environment. For example, only 
little is knovn about the variation vith heliographic latitude of so 
fundamental a flow parameter as the solar vlnd mass flux. Similarly, 
hardly anyv ing is knom about latitude variations of the solar wind 
energy and m oattun fluxes. Yet, these parameters may be very important 
in determining .he pliysical state of the polar corona and the size of the 
solar dominated cavity over the poles vhich separates the sun from the 
local interstellar meditun. In addition, the existence of heavy elements 
in the polar solar wind may depend (Alloucherie , 196?, 1970; Geiss et 
al. , 1970 ) on whether or not the proton particle flux exceeds a 
temperature dependent lower limit. 

It is therefore useful to consider hypothetical variations of solar 
wind particle and energy fluxes with heliographic latitude. This task is 
approached in section 2 of this paper by calctilating upper limit values 
of the polar solar wind particle flux implied by the most comprehensive 
coronal density model developed to date (Saito, 1970). This model was 
determined from an average K corona brightness distribution constroicted 
using 15 solar eclipse observations as well as K-coronameter measurements 
all made at the minimum phases of the solar activity cycle. As a 
necessary result of the method employed, the model densities (and hence 
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the upper limit values of the solar wind particle flux derived below) 
determined for the polar regions are uncertain because it is not rtc:«ible 
to uniquely invert the convolution integral which relates the coronal 
brightness distribution to the average line of sight electron density. 
Nevertheless » it is shown that if the polar coronal densities are as low 
as calculated \ising Salto's model, then without extended heating, the 
emitted polar particle fltuc should be substantially less than that 
observed in the ecliptic plane at 1 AU and less than that necessary to 
drag coronal heavy elements away from the sun. However, limited evidence 
based on Ly a measurements of the neutral hydrogen tall of comet Bennet 
(Bertaux et al., 1973; Keller, 1973), is consistent with a polar solar 
wind flux at leaJt as large as that observed in the ecliptic at 1 AU. 
These observations therefore require either an extended coronal heat 
soiirce distinct from electron heat conduction or that Saito's polar 
densities are too low. In any event since the 

particle flux in the polar wind may be comparable to that observed in the 
equatorial wind, it is possible that coronal heavy ions at polar 
latitudes do indeed expand with the protons into interplanetary space. 

Since it is reasonable to expect that heavy elements will be 
observable in the polar solar wind, the range of ionization state 
"freezing in" distances is estimated in section 3 for selected heavy ion 
species at polar latitudes. It is found that the polar coronal density 
may be sufficiently low that the ionization states "freeze in" below the 
nominal location of the temperature maximum. Hence high latitude heavy 
ion observations may allow a determination of the thermal state of the 
intermediate and low corona and provide an estimate of the magnitude and 



extent of mechanical dissipation. Section b summarizes the main 
conclusions. 

2. Latitude Variations of the Solar Wind Particle Flux 

It is currently thought that the solar wind ev'^lves from open field 
regions in the coroi.3i (see Hundhausen, 1972 for a review). Such regions 

are generally distinct from regions of activity and are generally characterized 
by low density. For th regions, the electron density, N, as a function 
of solar distance, r, and hellographlc latitude, 6, has been modeled by Salto 
(1970) with the relation 



where R ■■ r/R and R Is the solar radius. 

9 9 

Upper limit values for the polar solar wind particle fltuc can be 
derived using relation 1 if various subsets of several reasonable 
assumptions concerning the state of the intermediate corona are adopted. 
These assumptions are l) the coronal gas consists of H, Ke and electrons 
only, 2) there is no extended heating other than that due to 

electron heat conduction much beyond the coronal temperature maximum, 

3) the energy equation may be closed id.th the standard relation 
Q * (Chapman, 195^; Spitzer, 1956) which assumes that binary 

coulomb interactions limit the mean scattering length of a thermal 
electron, h) coronal electron and proton velocity distributions are very 
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nearly Maxwellian, 5/ wave-particle Interactions and macroscopic wave 
pressure effects are negligible above the heating region, and 6 ) the 
magnetic field is open but not necessarily radial. 

The purpose of this section is to show that if .he coronal*density 

over the pole drops off as quickly as implied by Salto's analysis then 

some of the above assumptions may be tested by in situ soleir wind 

observations. We begin with a standard single fluid formulation of t>e 

coronal expansion using equation 1 in place of an energy equation and 

derive upper limits for the particle flux at 1 AU. A separate treatment 

based on various possible forms of the energy equation is considered next 

to provide Independent estinptes of the 1 AU flux upper limit. The 

results of this analysis ai-e in agreement with those obtained by IHirney 

and Hundhausen (197^). As will be shown in section 3 these upper limits 

are substantially lower than that observed in the ecliptic at 1 AU and 

are sufficiently low, that if all of the above assumptions are correct, 

■H* 

He and many of the heavier ions should not expand with the protons away 
from the sim at polar latitudes. 


(i) Mass Flux. Momentum Flux and Density Equations 

The mass and momentum conservation equations are respectively; 


HVACR) = F 


NGM zn M 

^ p 



( 2 ) 

(3) 


2 

Here A(R) is the area of a fltix tube which varies as R If the expansion Is radial, 
G is the gravitational constant, is the mass of the sun, is the proton 
mass, M is the mean molecular weight * (1 -f* 4a)/ (2 ^ l»a) where a Is the He 
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abundance by number , k is Boltzmann's constant, B is the proton density, V is 
the bulk convection speed and T is the one fluid temperature. Concentrating on the 


region in the intermediate corona between R » r/R » 2 and 4, equation 
1 for 8 » 90^ can be simplified to the form: 


K(R) as 


7.9 X lO'" ^ -3 
/ cm 


If it is assumed that A(R) varies as w then equations 2, 3» and U can be 


integrated analytically to obtain T(r) 



Here the subscript o refers to parameters evaluated at the base radius 

In the following R in' chosen equal to 2- 
o 

Eauation 5 can be rewritten in slnplified form as follows? 

T(R) = (|-)^ {T^ - C^[l- )’^1- - ll) (6) 

o o 

Here and are constants which are readily evaluated iqr conparing 
equations 6 and 5. Inspection of equation 6 shows that T(R) depends parametrically on 
two variables, and (NV)^, Following the analysis of Brandt et al. ( 1965 ) 
it is possible to show that two physically reasonable assumptions imply 
stringent constraints on the range of realizable values of and These 

two assumptions are: l) the derived temperature, T(R), must remain positive through- 

out the range of validity of equation 1; according to Saito (l970), R < U, 2) there 
is not sufficient external heating >■ /ond R = 2 to produce a second peak In T(R). 
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It is seen from the third tern on the right hand side of equation 5 
that for a constant T^, increasing (NV)^ eventually drives T(R) negative. 
The radius at which this happens can be increased beyond R = ^ by in- 
creasing T^. However if is too large, the (R/I^)^ term in front on the 
right hand side produces a second peak in T(R) beyond R = 2. Therefore 

acceptable ranges of T and (NV) can be determined as follows* The 

o o 

minimum value of T^, is calculated for (NV)^ = 0 under the assumption 
that T(r) < 0 for R ^ R^ where R^ is the limiting distance of validity of 
equation 1. This gives 



(7) 


Given R = 2, M = 0.5**7 (corresponding to a hfj I{e abundance by number) and 
o 

assuming R^ = 3 and 3.5 then T. = 0*85 x lO^K and 0*89 x lO^K respectively* 

A ii 

Upper limits for aiid determined from equation 5 by finding 

the largest value of and (NV)^ such ^nat T{R) > 0 and dT/dR < 0 for R in the 

range R < R < R^* Thus for each R the following two relations must be satisfied 
® o X 




Inspection of eqxiations 8 and 9 shows that for (HV)^ = 0 both conditions 

can be satisfied simultaneously. However for each H both conditions cannot 

be satisfied if (KV) is leurger than some maximum value. This maximum is 

o 

obtained by equating the ri^t hand sides of equations 8 and 9* 


(HV). 


/GM f, /R /rV 


1/2 


s \ [7 


6 - / 


]\6 - 2 s) [\ 

»J V 

> 3 j[ 

U J; 


1/2 


2<S<6 (10) 


2 

Setting S = 2 (radial flow) and (MV) = (NV) (R /R^) (the subscipt e refers to 

e o o e 

parameters evaluated at the orbit of the earth) , equation 10 is plotted 
in Figure 1 for R^ = 3 and 3.5. The minimum value of the right 


I. 


r> 




- 
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hand side of equation 10 for R in the range R^ < R < is the maxiDuo flux 

at 1 AU consistent with the assumptions dT/dR < 0 for R > 2 and T > 0 for 

R < Rjp Thus for 3.0, (NV)^ <0.U x 10® cm"^sec"^ with 0.8| x 10®K < 

T(R = 2) < 1.1 X 10®K and for R = 3.5, (NV) < 0.2 x 10® cm“^sec"^ with 

e 

0.89 3C 10 K ^ T(2} ^ 0*98 X 10 K. The curves for T(R}^ corresponding to 
values of (UV)^ and T(2) determined from equations 9 and 10 evaluated near 
the minimtim of the curves in Figure 1, ere drawn in Figure 2. Drawn also 

A 

for conqparison are the polar scale height teaperature, ® (GM m M)/(kr dlnlT/Jr) 

and the curve T R ' . 

•2 

Since it is likely that the flow is more divergent than R inside of 

some radius, Rjj» it is. necessary to consider how this possibility affects. 

the upper limit of (NV)^* This may he accomplished by assuming the area 

S 2 

of a flux tube increases as R out to and then as R from there to 1 AU, 

R . Using this model, (NV),, = (NV) (R /R = (HV) (R /R and hence 
e D oop 

(NV)^ can be determined from equation 9 using the relation 


(NV)^ = 



R <R <R < R 
O X D 


( 11 ) 


Investigations of equations 10 and 11 for S in the rsmge 2 < S < U, 

= 2, (R,/R„) = 1.75 and R^ ^ R < R^ show that the maximum flux at 
1 AU consistent with a single temperature maximum below R is not significantly 
changed from its value for S = 2 ((NV]^< 0.2 x 10 cm* sec ■**). However, if 
S is siifficiently large and/or (Rjj^R^) is sxifficiently small, this upper 
limit is raised. For example choosing S = 5 with (Rj/R^) * 2 (which is 
equivalent to expansion from a polar region defined by 60^ < 0 < 90® at 
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6 

to the full hemlephere at I^), OlV)^ <0.28 x 10 cnT sec“ . It should 
be noted though that for all cases of nonradial expansion (NVl^ is significantly- 
raised o-ver that obtained for S = 2. 

ii) Energy Flux Supply to the Polar Wind 
An €aternative approach to the polar particle flux problem is possible by 
considering the energy equation. Here, a limit on (NV)^ may be established if the 
velocity at 1 AU is known and if the usual assuc^jtions about the state of the 
Intermediate corona are made. Using a one-fluid, steady-state, spherically 
symmetric model, the energy equation 


may be conbined with equations 2 and 3 and integrated to yield 


A(B) Q + F 




GM m MT 

0 p 

R R 

© J 


e = constant, 
o 


Cl3) 


If a supersonic solar wind exists at 1 AU but not at E then the dcninant tera at 1 

GM m M 

AU is F [4 m MV while at R , the dominant terms are A(R )Q + F kT ° ^ ]. 

2pe o o d c ^ w 

o o 

Therefore 


4 m MV 
2 P 


e 


2 



+ 



GM m M 
R R 


o © 


(iM 


Further progress is not ^possible without an additional closure 
relation which gives in terms of the lower velocity moments. Usually 
the Spitzer conductivity is assumed valid so that 

— 5/2 — 

Q = - K VT 
^ o 


- 117 -- 


( 15 ) 



vlth “ T.7 X 10 “'^ erg oaT^ sec”^ (Chapman, 195^, Spitzer, 1996). 

However, it is also possible that the density is sufficiently low over 
the poles that equation 19 is not obeyed. In particular, it is possible 
that the polar density is so low that the dimensionless third moment, 

1/2 

q » Q/[l.9 MKC(kT/m^) ^ ], becomes iaqjosslbly large at a low altitude 

(Parker, 196 U). For exanq>le if H^ * 1.23 x 10^ cm”^ (Saito, 1970), 

«= 0.98 X 10^ (see section (i) above) and T « then q » 

(0.15)(R/R^)^^^^^ or q > 1 when R > 1.9 It is therefore probable 

that below this altitude instabilities develop (Porslund, 1970) which 

will limit Q to a value less than the Spitzer upper limit. In other 

words, the heat flux will be limited within 1.9 R^ thus effectively 

producing an isothermal region at lower altitudes and a region of steeper 
-2/7 

than R temperature decrease at higher altitudes. 

It is thus not clear how to estimate the value of Q in equation 

o 

lU. For the saJce of concreteness tvo alternate approaches are adopted 

-2/7 

helow. The first asstunes equation 15 to he valid with T R and the 
second adopts an exospheric approach. In hoth cases the solar wind He 
abundance is asstJined to be a free parameter since its value is observed 
to be highly variable in the ecliptic at 1 AU and is not known at high 
polar latitudes. Such an assxmfiption is necessary since, in contrast to 
the analysis presented in section (i) where upper limits for (NV)^ were 
independent of M (see e.g. equation lO), the magnitude of the He 
abundance may be significant here. This fact results because most of the 
energy needed to drive the solar wind expansion goes into gravitational 
potential and kinetic energy which are both mass dependent. However this 
effect is more than compensated for by the fact that maximum values of 
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derived trcm the analysis in section (i) scale linearly with M (see e.g. 
equation 9)< 

Assuming first that equation 15 is valid, T « R**^^*^ and 

=(0.98 X 10^)(M/0.5*^7) K (see Figures 1 and 2) then upper limit values 

for (NV)^ can he calculated from equation ih for chosen values of M and 

V^, The results are summarized in Table 1 under the label (JJV)^ (Spitzer) 

for * 320, U50 and 750 km/sec and '! values corresponding to a He 

abundance of 0, O.OU and O.08 by number. 

In the appendix, an analysis is presented which shows that it is not 

clear whether or not electrons are ccllisionless below R = R^^. If indeed 

coronal electrons are collisionless near to but outside of 2 R then Q 

© 

must be calculated u.3ing exospheric theory (Jockers, 1970; Lemaire and 


ility of the 
ATE 13 POOR 
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Scherer, 1971a ,hj “Schulz emd Eviatar, 1972 i Hollweg, 197^^; Evlatar and 
Schulz, 1975). In this approach, only those electrons above the electric 
potential barrier, |eA<|>| vith velocities directed away from the sun can 
carry heat. The three fluid energy equations may be ccmbined and integrated 
to yield 

|eA4»| = M A 

where the A symbol signifies a difference between any two radial distances 

and e is the electronic charge. Choosing R and R « 1 AU as the two 

o e 

referenee-dfstascca Chan 


I e A(J) I — M 


1 ,, 2 
m V 
2 I) e 


+ 


0{ m 1 

® P I 

R R I 
o ®J 


(17) 


If both and the shape of the electron distribution at R^, f(V), are known 
then is readily evaluated using the relation 


^ ir/2 tt/2 

%‘J ] I 

Vj -ir/2 -ir/2 


m COS0 cos(J»)f (V)V^dV cos6d9d(t> (l8) 


where ^ m^Vg^ = |eA^|. Assuming a Maxwellian shajje for f(V) then 



( 39 ) 
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Equations lU, 17i and 19 can be combined to give a self-consistency 

condition for 3 « bence an upper limit for (NV)* if 

the bulk convection speed at 1 AU is to be greater than or equal to V , 

6 


Ciwlo 



-CS+2.5) 
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[3 + 7 + 


13.25 1 


where 


1 

2 kT 


04 Mm 
® P- 

R R kT 
o @ o 


1 

2 


and V >0 
e 


( 20 ) 


( 21 ) 


Using equations 20 and 21 along with the assumptions that = 2, = 1.23 

X 10^ cm”^ (Saito, 1970) and = 0.°8 x 10^(14/0,5^7) K (see e.g. Figures 1 

and 2), upper limit values for (NV)^ have teen calculated for various, values 

of M and V and are also listed in Table 1. 
e 

A comparison of the exospheric upper limits with the Spitzer upper 
limits for (NV)^ shows that if both the wave-particle collision frequency and 
the coronal electron density are low enough over the solar pole so that an 
exospheric formalism is appropriate, very severe upper limits can be placed on 
the solar wind flux at 1 AU whether or not the He particles expand with the 
plasma. These upper limits fall well below that calculated from the mass and 
momentum equations alone. However, since it is likely that the corona is 
sufficiently turbulent that an exospheric formalism is not appropriate, the 
true upper limits for (w)^ may be less than but closer to that calculated 
using the Spitzer conductivity, 

A comparison of the upper limit values for (NV)^ derived using the form 
of the energy equation which incorporates the Spitzer conductivity, with that 



derived using the mass and momentum equations* requires a knowledge of M and V^. 
Reasonable choices for values of these quantities are made as follows. First* 
inspection of Table 1 shows that by choosing the solar wind He abundance to be 
k% by number, the error made in estimating (KV)^ is probably less than 20%, 
Therefore, for the purposes of this comparison, M is chosen to be 0.5^7. 'Jon- 
cerning the speed of the polar solar wind at 1 AU several pieces of evidence 
have recently indicated that over the pole is higher than that observed in 
the ecliptic (Cole, 197^) and may be close to 750 km/sec (Gosling et al., 1976; 

Feldman et al., 1976). If this is the case then from Tabl' ', energy considera- 

/ . , 8 -1 * 

tions require that (NV) be less than approximately O.U x 10 cm sec . This 

6 

value compares favorably vith that derived using the n?.ss and momentum equations 

Op •! 

((NV) ^ (0.2 to O.lt) X 10 cm”* sec” ). It is therefore concluded that if heat 

conduction is the dominant node of energy transport at about 2 and if 

Saitohs polar density model is correct then the particle flux of the polar 

3 2 1 

solar wind should be less than about 0.5 x 10 cm'" sec*" , This upper limit is 

about a factor of 7 times less than the solar wind particle flux observed in the 

ecliptic at 1 AU (Feldman et al., 1976). 


3* Latitude Variations of Heavy Ion ”Freezing In” Distances 

Alloucherie (1967) derived an approximate criterion necessary for a 
heavy ion of Am^ and charge Z to diffuse upward in an expanding corona. 


It should be noted that this value is an upper limit. If the heat flux, Q, 

is regulated below 1.5 R as suggested earlier, then the region below 1,5 R 

o o 

becomes more nearly isothermal thereby reducing VT, Q, and hence the upper 

limit value derived for (lIV) . 

e 
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Kis result is: 


(NV)^> 



B/n e**lnX(Z®/A) 


( 22 ) 


where the sjnnbols are as previously defined and InX is the coulomb 
logarithm. Choosing ® 1.23 x 10^ cm""^, « O .98 x 10 K and expressing 

22 in terms of (WV)^ we get: 


(NV)^ > 
e 


2.0 X 10^ 
(Z^/A) 




(23) 


Since for a radial expansion (S = 2) this limit is approximately a 
factor of 1* to 10 times greater than the upper limit for (NV)^ derived above, 
it is reasonable to conclude that if Saito's model is correct and if the polar 
corona is not externally heated above r = 2R^, He may not expand with the 
solar wind. This conclusion remains -^-alid for S < U and (Rj^/R^) = 2 as veil. 

However, observation of the length of the neutral hydrogen tall of Comet 

Bennet (Bertaux et al., 1973; Keller, 1973) as a function of heliographic 

latitude indicates that the polar solar wind flux is at least as large as 
8 —2 “1 

2 X 10 cm sec at 1 AU. This value is in disagreement \/ith the upper 
limits deduced in section 2. It is therefore concluded that at lecst one 
of the assumption., made in the above analysis is not correct and that it 
should indeed be possible to observe solar wind heavy ions at polar latitudes 
at 1 AU. If true then measurements of the populati m densities of individual 
heavy ion ionization states will yield information concerning the temperature 
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structure of that region in the polar solai* coi’ona where the various ionization 
states "freeze in.” 

It is possible to determine the "freezing in" distances of the various 
heavy ion species as a function of heliographic latitude using equation 1 if 
the following assumptions are made; l) the flow is radial; 2) the velocity 
distribution is Maxwellian; and 3) the electron temperature, T, depends on the 
radius, r, as T = T (r/S Following previous work 

Q © 

(Hundhausen et 1968a, b; Bais> al., 1974) these distances are 

defined as those for which the expansion rate, « (VdlnN/dr) becomes 

larger than the ionization state changing rate, » N(R^ + C^). 

Hei: is the electron temperature at the base of the corona, V is the 

solar wind speed, is the rate of recombination from state i to state 
and is the rate of collisional ionization from state i to state i + 1. 

Changes in the "freezing in" distances with la^iitude of a sample of the 
most abundant ions are shown schematically in Figure 3 superimposed on scale 
height temperatures calculated using equation 1 for 0*0® and 90®. For 

purposes of illustration an isothermal corona with T « 1.0 x lO^K and a 1 

8 "2 ^1 •*! ^1 
i^U particle flux of 2.5 x 10 cm sec were assumed for evaluating + ^ci* 

The scale height temperature for a static corona is given by T « 

n 

2 

(CM. m M)/(kr dlnN/dr). Inspection of Figure 3 shows that the region in the 

o p 

corona for which temperature values can be determined from solar wind heavy 
ion data moves inward from above to below the temperature maximum as 6 
varies between 0® and 90®. Thus at some intermediate latitude, coronal 
temperatures bracketing the maximum can be sampled allowing the magnitude 
and extent of mechanical dissipation in the intermediate corona to be 
estimated (see e.g. the analysis of Brandt et al., 1965). 

^Collisional ionization and radiative recombination (including dielectronic 
recombination) coefficients for 0, Si, and Fe were kindly supplied by Dr. A. Dupree. 



l». Simmanr and Conclusions 

In this paper two related aspects of the physical state of the inter- 
plemetary plasma at high solar latitudes were explored. In the first part 
upper limits for the polar solar vind particle flux were derived using a set 
of reasonable assumptions concerning the base coronal conditions along vrith 
Saito's (1970) coronal density model. In the second part, it was determined 
whether this flux was sufficient to drag the heavier ions away from the «<un 
into interplanetary space. 

From the analysis in the first part it was concluded that if Saito's 
model is correct, the polar electron density is sufficiently low that in the 
absence of extended heating the solar wind flux at high latitudes should be 
at least a factor of from to 10 times less than that observed in the 
ecliptic at 1 AU. Such a low particle flux was shown in the second part 
to be small enough that most heavy ions would not be expected to expand 
with the protons into interplanetary space. 

However, indirect and limited evidence available at present is con- 
sistent with a polar solar wind that has at least as large a velocity (Coles 
et al., 197^; Brandt et al., 197^^) and as large a particle flux (Bertaux 
et al., 1973; Keller, 1973) as that observed in the ecliptic at 1 AU. From 
the analysis presented in section ??, these observations then require either 
that extended heating distinct from that provided by electron heat conduction 

is necessary some of the time above 2R or that Saito's polar densities 

o 

are too low. Whichever is the case, the fact that the 

solar wind particle flux does not appear to decrease with increasing heliographic 
latitude (Bertaux et al., 1973; Keller, 1973) indicates that coronal heavy 
ions may be expected to expand with the protons away from the sun. If true 
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then measurements of the population densities of individual heavy element 
ionization states in the polar wind vill provide information at 1 AU con- 
cerning the thermal state of that region in the intermediate corona where 
the respective ionization states freeze in. It turns out that the latitude 
variation of these freezing in distances calculated using Salto's model 
Is such that the region in the corona for which temperature values can be 
determined moves inward frca above to below the nominal location of the 
temperature maximum as 0 varies between 0® and 90°. Therefore, measuicments 
of heavy ions at high solar latitudes may provide valuable information con- 
cerning the magnitude and extent of mechanical dissipation in the intermediate 
polar corona. 
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Table 1 

Upper Linit Values of 
Consistent With the Ener?:y Equation 


km/sec 

$ 

-2 -1 
cm sec 

-2 -1 
cm sec 

V 

e 

He/K 

{lIV)^(Spit2er) 

(rrv) Exospheric 
e 

320 

0 

1.02 X 10® 

1.2U X 10^ 

1*50 

0 

0.70 X 10 

h.96 X 10** 

750 

0 

0.33 X 10® 

0.1*1 

320 

U 

D 

1.28 X 10 

1.3 X 10^ 

1*50 

1* 

0.88 X 10® 

u 

5.2 X 10 

750 

1* 

0.1*1 X 10® 

0.1* 

320 

8 

O 

1.5l* X 10 

1.1*0 X 10^ 

1*50 

8 

1.06 X 10® 

5.59 X 10** 

750 

8 

O 

0.1*9 X 10 

0.1*6 
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Appendix 

Comparison of Expected Electron-Electron Collision Lengths 
with Scale Lengths in the Polar Corona 
The magnitude of the electron conductivity in the polar corona 
depends critically on the electron-electron collision length > If 

is small enough then the Spitzer conductivity is applicable but if it is 
too large, then an exospheric approach is needed to evaluate the polar 
electron heat flux* It txims out that, according to Saito, is 
sufficiently low over the pole that it is not clear vhether or not 

electrons are collisionless above R = 2R . exanrie the self scattering time for a 

© 

thermal electron at 2R^ is = (l.l x = 87 sec (Spitzer, 1956) 

whereas at that distance the expansion time (assuming a radial magnetic 

field) is T = [(k? /m )^'^^dlnN/dRj“^ = 6o .5 sec. Furthermore the coulomb 
e 0 e 

scattering length, defined ^y 



(Al) 


may be either larger than or smaller than the temperature scale length, 

= I-dlnT/dR] ^ = 3-5 R^, depending on the value of the maximum altitude 
at which Saitohs density model is valid, R^. This is readily shown by 
assuming a density model consistent with Saito’s results (19T0): 




with Rjj ^ U* Combining and A2 and using = 2 with from equation 1 
evaluated at 0 = 90® it is found that “ 1*33, 2.U8 and 6.69 for = 

3.0, 3.5 and U.O respectively. 

Since the actual value of I may be either less than or greater 

c » 

than depending on the value of R^, it is not known whether the Spitzer 
conductivity or a conductivity calculated using exospheri^c theory is most 


valid in the polar solar corona. However, the fact that is of the same 
order of magnitude as sioggests that neither of the above is correct and 
that to obtain an accurate determination of the true conductivity a 


kinetic approach may be necessary. 
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Figure 1. 


Figure 2. 


Figure 3. 


Figure Captions 


Plots of equation 9 for = 3 and 3.5* Radial flow is assumeu 

2 

and hence (NV) = (HV) (R /R ) . The lainimum value of each 
e o o e 

curve corresponds to the maximuia flux at 1 AU consistent irtth 
the assmptions dT/dR < 0 for R ^ 2 and T > 0 for R ^ R^^, 

Coronal temperatures, T(r) for values of (NV)^ and T(2) determined 
frcMi equations 9 and 10 evaluated at the minima of the curves in 
figure 2. Drawn also for comparison are the polar scale height 
taaperature, and the curve T « 


Variations with- latitude of the ^freezing in" distances of a 
S 83 iq>lc of the most ah\mdant heavy ions. Scale height tempera- 
tures are calculated using equation 1 for 0=0° and 90 ° and a 
constant "freezing in" temperature of = lx lO^K is assumed. 
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AJSSTRACT 

The large-scale, three dimensional magnetic field in the interplanetary 
medium is expected to show the classical spiral pattern to zeroth order. 
However, systematic and random deviations can be expected, although their 
nature and magnitude cannot be predicted. The sector structure should 
be evident at high latitudes, but the actual extent is unknown and the 
shape of the sector boundaries is controversial. Interplanetary 
streams will probably determine the patterns of magnetic field intensity 
but the actual patterns cannot be calculated at present because of our 
limited knowledge of speed profiles and the source conditions. The 
large-scale spiral field can induce a meridional flow which might alter 
the field geometry somewhat. The non-uniformities caused by streams will 
probably significantly influence the motion of solar and galactic particles. 
Unambiguous and detailed knowledge of the 3-dimensional field and its 
dynamical effects can only be obtained by In situ measurements by a 
probe which goes over the s.;n*s poles. 
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I. IK'n;ODlT.TTC\’ 


One cannot be c* ;ain of what will be observed on an out-of-thc-ccliptic 
nission. It is basically c.'q>loL*atory. One CuXI try to predict what will 
be seen, using current theories and the available interplanetary observations, 
and tliis paper attempts to do so for the interplanetary suignccic field. 
However, extrapolations to as little as 10*^ above the ecliptic a'-c highly 
uncertain. Only in situ measurcacuts con provide us Che unaribiguous and 
detailed knowledge that we seek. 

Many of the properties of the magnetic field observed in the ecliptic 
plane follow from a simple relation which is valid when the magnetic stresses 
arc not so large as to appreciably alter the motion, viz. 

Hr) - :?(r)/p^l (1) 

where B(r) is the field in a radially moving voliime element with constant 
speed, p is Che density, ^ and are the field and density at some surface 
near the sun, and V X is the gradient of the displacenent vector which is 
detemined if the speed is Icnown on the source surface* Thus, if 
end are known at scnac inner boundary (say 0.1 or 0*2 AU), and if p(r) 
is known, then to good approxiitation one can project or map the field any- 
where within i. AU by (1), if the magnetic stresses can be neglected* Ihis 
approach to the interplanetary magnetic field is discussed in Schatten 
(1972), Nolte and Roelof (1973), and more generally in Burlaga and Earouch 
(1975) and Barouch and Burlaga (1975a)* It is valid for both time-dependont 
flows and steady flows* 

Very little is known about the latitudinal variations of V near the sun 
(V (0)) so one must be content to explore several reasonable alternatives 
in order to study the effects of V X on the three dimensional field. Of 
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course, the sisqilcsC ossuention i;: V-constant, whicii :;Iveu the spiral field, 
as discussed below. 

Hie values of Bj^(r^j,o,t) ca:i be estimated by projcctir .3 photospheric 
taagnetic field neasurca^ts upward tlirou^i the solar envelope to the Alfven 
point and beyond. Several techniques for doins this are available in the 
literature, althou^ none is coc^letely satisfactory. Tne problcn is 
illustrated in Figure 1 from Schatten (1971), which shows a sketch of an 
eclipse in which the lines presumably r^rcsent magnetic field lines. One 
sees a variety of structures. Near the sun, the field is complc:: with many 
closed loops visible at low latitudes. Farther from the sun, the field 
lines generally tend to diverge and to become nearly radial. This is 
represented formally by assuming potential fields near the sun and supposing 
that only the lowest order hamonics contribute to the interplanetary field . 
The transition to radial fields is generally made artifically at some 
distance, and structures such as helmet streamers can be modeled by post- 
ulating that currents are present only in thin sheets. Such methods are 
sufficient for us to make estimates of the field out of the ecliptic, but 
it must be ec^asized that these are only approximations based on models 
rather than firm theoretical predictions. There is no substitute for 
in situ measurements of B out of the ecliptic. 

II. MA.GNSTIC FIELD DIRECTION 

Parker (1958) presented a model for the zeroth order configuration of 
the magnetic field lines, assuming constant V^, B^, andp^, and steady 
corotation. In this case, (1) gives the well-known Archemedlan spiral. A 
good illustration of this pattern may be found in Hirose et al. (1970). 
Measurements have shown that Parker's modal given an acceptable zeroth 
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order approximation for the field in the ecliptic plane between 0.3 AU and 
S AU (see the review by Behonnon, 1975, and the paper by Marion! et al, 
(1975)). There is appreciable scatter of the observed points about the 
theoretical curve, whicli miglit be due in part to variations in ^ and 
(Burlaga and Barouch, 1975), and one can expect to observe similar scatter 
away from the ecliptic. One might also see systematic effects in the 
direction of B due to systematic variations in B^ associated with structures 
such as helmet streamers and polar plumes. 

Stenflo (1971) attosqited to confute realistic interplanetary magnetic 
field configurations by introducing a recisonable model for V(r) near the 
s\in, projecting measured photospheric fields to a source surface at 2.6r^ 
using the potential field mapping technique, and mapping fields from the 
source surface to 1 AU by a technique equivalent to Eq. (1) . An ic^ortant 
result is that althou^ there are coi<ralex 1 nop- configui^^t Ions close to the 
sim, farther from the s\m Che field becomes more radial and only the lowest 
harmonics are significant for the solar wind flow. In particular, the 
large-scale (m 1 AU) ,3-dimensional field which he c<mq>uted for the period 
Feb. 17 to March 10, 1970 was found to have the spiral form predicted by 
Parker's model (see Figure 2). 

The pattern that we have been discussing is altered by the presence 
of streams. The magnetic field lines will be more radial when the speed 
is hi^ than when it is low. This effect is small, however, being Just 
a few degrees in the ecliptic plane (e.g. , see Burlaga and Barouch, 1975) 
and probably even smaller at higher latitudes. 
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Ill, HXCNETIC FIKLD SCTSE (S1T.T0RS ) 

Wilcox and Ness (1965) showed that the interplanetary magnetic 
field is structured such that it tends to point away from the sun for 
several consecutive days, then abruptly changes direction by 180*^ and 
points toward the s\m for several days, etc.* In other words, the 'Mense” 
of the field is sectored on the mcsoscale. In 1964, four sectors were 
observed by IMP 1 (Wilcox and Ness, 1967). Several papers, beginning 
with that of Wilcox and Ness (1965), show that the interplanetary sector 
structure is related to the polarity of the photosphcric magnetic field. 

In a recent study of this sort, Kolte (1974) (see Roclof, 1974) computed 
the cross-correlation (calculated over nine solar rotations) between the 
interplanetary polarities which were mapped to their connection longitudes 
on the sun using the solar wind speed and an empirical technique an'* the 
chromospheric polarities. He found that the cross-correlation peaks 
at latitudes +30^ and -20^, suggesting that the base of the field lines 
is generally not in the solar equator; however, the correlation was low, 
;^«3« \4hen he computed the cross-correlation coefficient for magnetic 
field data corresponding to speeds >400Ion/sec, he found a peak of ^*5 at 
a latitude very near to the solar equator, suggesting that the streams 
which are observed originate near the solar equator. 

Ness and Wilcox (1967) showed that the sector structure changes with 
time. Two sectors were observed in 1962, four in 1964, and the structure 
was complex in 1965. Similar plots for 1970-1972 are given in Fairfield 
and Ness (1974) together with references to other papers on time variations. 
Generally, the pattern of 2 or 4 sectors is the dominant one. 
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Wilcox and Svaalsard (197') co.v idi.-rcd a;: ratcrc.-iCin^ narccm in 
1969, when 2 sectors were prodov.inant for several solar rotations. 
Comparison with the solar fields, as determined by the “hairy-ball** 
(potential field napping) model, showed that the sector boundary, projected 
from I AU to the sun, corresponded to an "arcade" of closed loops running 
approximately K-S, Taus, one expects that the sector boundary extends 
to rather high latitudes in this ease. 



- 148 - 


pri’pnnUCIBn.rrY OF THS 

page is rcK>a 



dipolo whoso axis is in the equatorial plane, is t’no doninant pattern, 

Tltcy also found that four sectors occur frequently and that a :>S dipole 
predominates only occasionally. Hius, one expects that generally there 
will be two sectors whose boundaries extend N-S, soraetines four sectors 
will be present, and occasionally the polar regions will tend to he 
unipolar with a "sector bouncar/* in the solar equatorial plane. Iho 
time variations can bo very irrjortant for a S/C nission lastin’ rare than 
a year. It will be important to correlate measurersents of the solar 
field, out-of-ecliptic ncasurenents , and measurements made near the 
ecliptic, to separate the space^tinc variations 

Svalgaard et al. (1974) proposed a model to describe the more complex 
configurations in which both a N-S dipole and a dipole whose axis is in 
the equatorial plane contribute to the interplanetary magnetic field. 

This is illustrated in Figure 3 (bottom), which shows that the polar 
fields and the equatorial fields night combine to give sector boundaries 
that are tilted with respect to the solar equator, the direction of tilt 
depending on the polarity of the fields, Ihe figure also illustrates the 
magr.etic arcades and the associated helmet streamers that arc presumed 
to be associated with sector boundaries. Svalgaard et al. ( 1974 ) 
presented evidence that strongly supports this picture for the period 
that they studies. An alternative model (see Figure 3) was presented by 
Hansen et al, (1972) based on data obtained in 1972 , when 4 sectors were 
observed at 1 AU and coronal streamers were observed by OSO-1. hansen 
et al. postulate a coronal bridge corresponding to the low-latitude 
arcades of Svalgaard et al., but they differ in postulating independent 
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atrcamcrs at hi^ latitudes. The relative merits of these two models 
has not been established. Out-of-ecliptic measurements could be decisive. 

Rosenberg and Coleman (1969) showed that the number of days with 
negative polarity varied sinusoidally with a period of 1 year between 
1964 and 1967, randomly in 1968-1969, and sinusoidally with a 180° phase 
shift in 1970-1973. This pattern showed the expected change in 1974-1975 
(Fairfield and Ness, 1974). Following a suggestion of Rosenberg and 
Coleman (1969), Schuls (1973) suggested that this is a consequence of 
"warping" of the equatorial plane (minimum B surface) of a dominantly 
N-S solar dipole with significant quadrupole contributions. This conflicts 
with the results of Altschuler et al. (1974) which show that the N-S dipole 
is not frequently observed at the sun. It also implies that the sector 
boundaries should trace a "sinusoidal" line with small amplitude about 
the solar equator, in contradiction with the results of Svalgaard et al. 
which suggested that sector boundaries are not tilted so much. This 
difference has not been resolved, but could be settled by an out-of-ecliptic 
mission. 

IV. MAGNETIC FIELD INTENSITY 

Parker's model, based on the assiin^tion of constant and uniform V 
and B^, provides the zeroth order approximation of the interplanetary 
magnetic field intensity. It predicts the measured value of fs 5/ in the 
ecliptic plane at 1 AU for reasonable values of B^ near the sun, and it 
predicts somewhat smaller fields near the poles. A somewhat more 
complicated model, assuming constant V and a N-S dipole giving B , was 
discussed by Parker (1958) and considered in more detail by Stem (1964) . 
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In view of the results of Altschuler ct al. (1974), which showed that a 
K-S dipole rarely dominates, this model is generally not appropriate. 

Superimposed on chc larsc- scale variations in Tna<pactic field intensity 
are non-uniformities due to streams. These are the result of 7 X i 0 in 
Eq. (1). Faster plasma ovcrtoltcs slower plasma, causins a compression of 
the plasma and (because the field is "frozen" to the plasma) an cnliancement 
in B. Shears in V can also cause a chansc in the magnetic field intensity, 
and this too is in^licit in (1) (see Burlaga and Barouch, 1975). Enhance- 
ments in B are generally observed at the leading edge of streams at 1 AU, 

often as large as four times the ambient value. Similar enhancements 

night be observed out of the ecliptic, depending on the velocity profiles. 

They night be the most important magnetic field intensity variations 
at high latitudes. Illustrative spatial configurations of the magnetic 
field intensity on a spherical shell v;ith radius 1 AU, relative to the 
unperturbed equatorial field at 1 AU, are shown in Figure 4 from garouch 

and Burlaga (1975a), At the top, is the result for V * (1-A cos 4,5) 

- . - -2 

cos 0; at the bottom is the result for V = (1+A cos 44/)e:o*-(5-=^) . 

In the first case, one expects that the field intensity might be approx- 
imately twice the unperturbed intensity in some regions out of the ecliptic, 
if the streams extend to nigh Ivtituces. In the second case, in v.hich 
streams are confined near the ecliptic, B(8) is strongly perturbed only 
near the ecliptic, 

Billings and Roberts (1055) sugp.estcc dial: streams come from r.^-ions 
of open and diverging magnetic field lines near the sun and that slow 
plasma is associated with closed loops. Tnis is consistent with ti.e 
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obiSi'rvation that the solar wind speed is ijencrally snail near r,o;;l:or 

boundaries, which according to Svalgaard ct al. (1974) are associated 

with "arcades" of closed loops near the sun. It is also consistent with 

mappings of streams back toward the sun (Roelof, 1974). Pneunon (1973) 

and Pneuman and Kopp (1970, 1971) modeled this situation with a dipol<'‘ in 

the solar equatorial plane. The basic idea is that when the field lines 

diverge, heat is readily conducted to the critical point where it con 

effectively accelerate the solar wind, whereas when the field lines are 

closed heat cannot be conducted radially because of the low por*pcndicular 

conductivity and energy is not available for acceleration. Of course, 

other models are also possible. TIic point that we wish to emphasize is 

that stream profiles might be related to the magne’ field near the sun. 

Calculations which explore the consequences of variable B^(9,({) arc given 

in Barouch and Burlaga (1975a). Tlie result is that there might be a 

latitudinal variation of the stream induced perturbations in B which 

results from the latitudinal dependence of B . 

* o 

V.'e conclude that streams induce ^ijr.if icoict distortions in the ;ic 

field intensity which must be considered in r.easuring a "zerotl." order 
field. They are also interesting in tiiem.selves and r.ay help to und. r^tend 
the source of streams. Finally, they have important effects on sol-r 
particles, galactic particles, end possibly s. ’ ar wind flow perturbation.. , 
as discussed in the next section. 

V. ImFECTS 0? 3(9, (^) Gb COSl^lC :f.VS PLASM. 

Several authors (e.g. see the reviev.’ of Montgomery, 1971) seggesu.^ 
that particles should have easier acce.c.o and .shorter pctr.s if L .y 
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'Ti'iCi presence of strear.-in ' :coc ^reeients in :: can also significantly 
affect cosnic rays* Barouch and Buriaga (1975b) showed that Forbush 
decreases and similar s^l^ctie cosr.ic ray intensity variations are strongly 
correlated with magnetic field ernanccncnts associated with streams. 

Tliey proposed that these are the results of perpendicular gradient drifts, 
and Barouch and Burlaga (1975a) shov;ed that the drift speeds are 
appreciably higher than the speeds in streanis, as required* If strear.'.s 
extend to high latitudes, one e:cpects to observe Forbush decreases there. 

If streans do not extend to nigh latitudes, there will be small drifts 
due to the spiral field configuration (VJinge and Coleman, 1968), but there 
may be no observable effects. 

Streams cause field lines to diverge less rapidly and eventually 
converge in the conrpression regions in front of streams, faus, solar 
particles moving in such nesoscalc conf igurations v;ill be collimated less 
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stronjly thai iii the larso-sci'.lc 3-.)i:*al field, :i:ul r.o::.i. v,.Mi r.i/ror 
(Barouch and Burlaga 1975a), giving jurean related i-.-.Cen:;i:y and a;jir;otropy 
profiles. One ocpects systcr.atic differences in tlieae profiles wirl. in- 
creasing latitude, depending on the variations of the streans. 

Low energy (thermal) particles arc also influenced by tlie magnetic 

field, although not as strongly as the hi^h energy particles vd esc .a'rgy 

2 

density is much less than B /(Sn). In fact, Parker's model, i.vici. . unes 
constant V and gives a spiral field, is not exactly self-consistent for 
this reason (Gussenhoven and Carovillano, 1973; Aleksoycr ct al, 1071). 

In oarticular, the spiral field gives a Jx_B “ C7>:b)x3 force that causes 
a meridional flow away from the ecliptic, Tliis flow has been studied by 
Kinge and Coleman (1974) and by Sucss (1974), Its magnitude might bo as 
much as » 1 kia/sec in the ecliptic at 1 AU. Rosenberg and Coleman (1973) 
invoked this flow and the frozen field condition to explain his observations 
that the magnetic field direction diverges av/ay from the ecliptic plane. 
Tnesc effects are small and have not been confinned. Tney vary with 
latitude in the spiral field configuration, and they might be strongly 
modified by streams. 

VI. SUIXARY 

One expects the large-scale, chrce-dimensional magnetic field lines 
of the solar wind to have the form of spirals wrapped on cones, as 
described by the solutions of Solar wind streams and solar 

magnetic field configurations probably will not alter this very much, 
although small, systematic effects due to the variation of the orientation 


uf B near the sun might be observable. 
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Tlic sector pattern possibly extends to latitudes and can change 

appreciably during a year. Tlic pattern and extent of sector boundaries 
is a matter of controversy. Extrapolations of the solar field and mapping 
of the interplanetary field to the sun suggest that the boundaries extend 
nearly north'SouCh, although tilted somewhat depending on the polar and 
sector field directions. On the other hand, the "Rosenbcrg-Colcmcn 
dominant polarity effect" and some calculations suggest that sector 
boundary surfaces are confined closer to the equator. Two distinctly 
different models of sector boundaries have been proposed by Svalgaard et 
al. and Hansen et al., but one cannot chose one or the other at the moment. 

The interplanetary magnetic field intensity will vary with latitude 
depending on the photospheric field configuration. For a solar monopolc, 
the polar field near 1 AU is a/2 smaller than the equatorial field. The 
intensity might vary more than this on a smaller scale, < 1 AU, due to 
the presence of streams. The actual configuration depends on both the 
stream profiles and the magnetic field intensity profiles near the sun, 
but these are not lenown at present. 

The gradients in magnetic field intensity produced by strea-ms cause 
energetic particles to drift away from the ecliptic, and they might be 
responsible for Forbush decreases. If so, these decreases siiould disa7,pear 
at high latitudes if the speeds arc confined near to the ecliptic. Stream 
induced magnetic field cT.h.ance;.'.;;'. t.^ r.igi;t also mirror soisr p./tic., ^ 
a..J Ti ..ai i' l!u:ir col 1 i ar. , c..’,:.; i :v’, t ...M- rcl .it c d i.. intan.ity 

..nd a.'.isotropy. The effect varies v.'it'r. latitude in a systemawic way, oO 
t..at cr.w can test t*.csc .....^as ’.’1^,'. ..a out"0.*t.'.C'"ecl mtec 
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Tlic ! •‘rgc-scalc spiral fiolJ causes a small neridlorjal flov; w, a 
consequence of the force, ilic v.;a;^aitu(lc of this no’..' be 

altered by streams and vary nil',, latitude for this I'cason. 

In conclusion, an out-of-tho-ccliptic mission will allow us to test 
present models of the Interplanetary magnetic field, resolve some 
controversies, provide iiiforr.;atioa needed to understand encr”,ttic 
particle and plasma motions, and it will probably give new results that 
we cannot anticipate. 
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Figure I 


Figure 2 


Figure 3 


Figure 4 


Sketch of a solar eclipse on 30 i:ay 1005 . The contours 
arc believed to indicate the direction of the magnetic 
field 

Interplanetary magnetic field lines on a scale of 1 AU, 
seen by an observed in the ecliptic plane. They were 
computed by Stenflo using photosphcric magnetic field 
measurements. 

Sector boundaries. This illustrates two conceptual models 
of sector boundaries and their relation to coronal streamers. 
Magnetic field intensity contours relative to the un> 
perturbed intensity in the ecliptic plane at 1 AU, on a 
surface with radius 1 AU. The top figure shows the 
pattern caused by a scream which varies with latitude as 
cos 9, and the bottom figure describes the result of a 
scream which is confined near the ecliptic. 
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Abstract 


Initially spherical blast waves are systematically distorted due to persistent 
latitudinal solar wind structure. This is to be distinguished from the non-systematie 
(random) distortion due to varying structure in longitude which introduces an ^ 30^ 
average deflection of shock normals from radial. The systematic latitudinal effect 
should be at least 25^ at mid-latitudes, and observable in the 30° noise level with 
14 or more shocks for statistic. The observed occurrence rote of shocks during 
solar maximum is sufficient to detect the effect. Corotating shocks should become 
deter:table between 1 and 5 AU. Identification could be a problem because of the 
30° noise level becoming greater beyond 1 AU. However, the three-dimensional 
geometry of corotating shocks show a strong latitudinal structure which can be used 
in an out-of-the-ecliptic mission for a statistical identification based on shock 
occurrence rates. 
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IntroducHon 


Most of the interplanetary shock waves observed with 1 All of the sun 
originate from some short lived solar event, such as a solar flare, and then 
propagate out as a more«or*less spherical shock wove until they leave the solar 
system. Beyond 1 All another class of interplanetary shock wave becomes commonr- 
the cortating shock pair formed by the interaction of long lived solar wind streams. 
We discuss here the three dimensional geometry of these two classes of interplanetary 
shocks and how these geometries can be studied with an out^of-the-ecliptic mission. 
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Oot-of“the-EclipHc Distortion of Solar Blost Waves 


Lack of spherical symmetry in the solar wind distorts the surface of a shock 
wave as it propogctes away from its solar origin into interplanetary space. This 
phenomenon is observed near the solar equatorial plane at 1 AU, where the inhomo- 
geneities associated with solar wind streams produce typically a 30^ deflection of 
the shock wave normal away from the radial direction ( Heinemonn and Siscoe , 1974; 
Hirshberg et o l., 1974). The distortion results from differential advecrion of the 
shock front due to solar wind speed variations and from differential propagation 
speed of the shock in the solar wind due to density variations— the shock propagates 
more slowly in gh density regions. In spite of the large distortion of individual 
shocks/ the averaged shock normal direction near the equatorial plane ot 1 AU is 
radial from the sun ( Chao ond Lepping . 1974; H einemonn and Shcoe , 1974). 

A systematic variation in latitude oi fhe solar wind speed and density pro- 
duces a systematic distortion in latitude of the surfaces of solar produced shock waves. 
That is, the averaged shock normal direction (n ) in general will not be radial from 
the sun. The angle between n^ and the radial direction will depend on latitude in 
a manner which reflects the average latitudinal dependence of the solar wind speed 
and density. 

A lower lim t on the deviation of n^ from radial is shown in Figure 1, Here 
a shock wave that was spherical at 20 solar radii becomes distorted into a quasi- ellipse 
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at 1 AU by the action of a differential advection of 400 km sec ct the equator 
increasing smoothly to 600 km sec ^ at the poles. This latitudinal gradient of solar 
wind speed, approximately 2 km sec ^dsg \ is at the lower end of the range of 
gradients suggested in the literature (see reviews by Gosling , 1975, Hundhausen, 1975, 
Dobrowolny ond Moreno , 1975). The figure probably also underestimates the dis- 
tortion for the assumed gradient since it neglects the possible latitudinal gradient in 
density-decreasing toward the poles (Huncfttausen et a l. , 1971)— causing the shock to 
propagate faster in the solar wind at higher latitudes. 

The maximum deviation of n from radial (AO ) occurs at mid-latitude a-d 

s max 

is about 25^ at 1 AU. In order to observed this effect, enough solar generated shock 

waves must be measured while a spacecraft is in the mid-latitude region to obtain a 

value for the polar angle (d^) of with o statisticol error considerdaly less than 

Ad • For the soke of having a numerical e,:ample, we take the requirement that 
max 

the expected standard deviation of 0^ be less than or equal to 1/3 of the conservative 

value obtained above for A0 , i.e. S. D. (0 ) s8°. 

max s 

Near the equo^orial plane the distortions produced by solar wind streams 
cause approximately a 30° standard deviation in the ongle between individual shocks 
and their average, radial direction. The standard deviation of the average angle of 
N silocks is 30 //N , For N = 14, this is ~ 8 . Thus in a 30 backgrounc >oise 
level ir the angle of individual shocks due to stream structure^ it takes 14 shocks to 
obr^ir a value of 0 with an expected stondarc’ deviation o^ 8° from the true value. 

It is likely that the noise level at mid-latitudes is less than 30° since the effect of streams 
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in producing solar wind inhomogeneiHe$ as measured by enhanced radio scintiiioHons 
is apparently confined to within 40*^ of the equator (Houminer^ 1973). Thus the 
requirement of 14 shocks is probably more thon necessary. 

To estimate the time needed to observe 14 solar generated shodcs, we use 
the 100 years of SSC statistics compiled by Moyoud (1975), The validity of this 
procedure is based on the study of Chop ond Lepping (1974) showing thot at leost 
87^ of SSC's can be associated with solor activity such as solar flares and type 2 
and type 4 radio bursts. There ore on average about 10 SSC's per year during solar 
minimum and about 35 SSC's per yew during solar maximum. Thus about 17 months 
we required in the first case arxl 5 months in the second of mid-latitude observations 
to observe 14 shock woves originating from a solar surface source. We conclude that 
an out-of-the-eci iptic mission scheduled for solar maximum would have high probability 
of observing systematic shock wave distortion even if the worst-case example discussed 
oboveshould apply. 

The impwtance of measuring the systematic distwtion of the shock shape lies 
in its use in determinirtg the systematic latitudinal dependence of solar wind parameters. 
This method is independent of all other methods. It does not involve unravelling 
separate space and time vwiations. In the case of an out-of-rhe-ecUptic mission via 
Jupiter, the space dependence involves both radial distance and latitude angle which 
cal' give independent contributions to any variation : observed in in situ measurements. 
Knowledge of the three dimensional shape of blast waves is important also for determining 
the flow of flare energy into interplonetary space. A non-spherical shock shape implies 
that . crgy is not distributed uniformly but converges in some places — relative to 
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a purely radial flow— and diverges in others. In the case considered, energy would 
diverge away from the equotor ^otning the shodc strength to decrease faster in the 
equotoriol plane thun would be expected on the basis of sphoricol symmetry. 

Three Dimertsionol Structure of Corotating Shocks 

Near the equatorial plane the border between contiguous solar wind streams 
is a spiral (Sard>hoi , 1963; Dossier , 1967). If the trailir^ stream — in the sense 
determined by Hie direction of the solar rotation— is faster, a pair of shock waves 
will form at some distance from the sun (Hundiousen , 1973 o,b). Such shock wave 
pairs hove a parently been observed between 1 and 5 AU by the Pioneers 10 artd 11 
spacecraft (Hundhousen ond Goslir^ , 1975; Smith ond Wolfe , 1975). In a steady 
state situation, the streams, their spiral border, and the shock pair all corotate with 
the sun. In this section we estimate the heliocentric distance of shock formation in 
the equatorial plane os a function of the speed differential between the streams and 
give a qualitative description of the three dimensional shape of the shock surfaces. 
For the latter we assume that the stream border is perpendicular to the equatorial 
plane. This example illustrates the essential aspects of the geometry and possibly 
represents the typical case as indicated by radio scintillation observations (Houminer , 
1973) and the north-south alignment of coronal holes which might be the sources of 
fast streams (Krieger et al . , 1973; Noyes , 19751. 

Figure 2 shows the relevant geometry in the equatorial plane. The figure 
also illustrates one argument for expecting the existence of shock pairs which at the 
same time suggests a simple calculation for the approximate heliocentric distance to 
their point of formction. The spiral labeled stream interface is the border between 
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the streoms. Sample flow streamlines— labeled fast and slow~ore shown in the corotating 
reference frame in which all geometrical features are time independent. With the fast 
stream trailing, the pitches of the spiral streamlines in both streams are such that they 
would intersect the spiral interface unless prevented from doing so by forces that act to 
deflect the flows away. The build up of pressure at the interface resulting from the con- 
vergence of the flow there produces such a force (Siscoe, 1972). Relative to the flows 
in the two streams, the interfacw looks, like a curving wall and the flows are forced to 
follow the curve because of the increased pressure. The compressive deflection of a 
supersonic flow by o curving wall is known to produce a detached shock wave in the 
flow ( Landau ond Lifschitz , 1959, p. 429). Streamlines intersecting the shock waves 
are deflected parallel to the interface spiral, bringing to an end the compressive 
interaction between the streams. 

The opproximate location of the origin of the shock waves can be found by 
considering the characteristics of the flow eminating from a point. A, on the interface 
close enough to the sun that the streamlines are essentially radial but far enough from 
the sun to be in the supersonic region. The choract eristics ore generated by foil }wing 
the progress of a sound wave starting at A and subsequently expanding and being con- 
vected with the flow, as illustrated in the figure. However, before the shocks ore 
formed, the flow converges cn the interface; thus ^he sound speed must be greater than 
the speed characferizirtg the convergence in order for the sound wove to expand. As 
the wave moves out from point A, the sound speed decreases because the solar wind 
cools as it expands, and the speed of convergence increases because of the relative 

' 1 1- u y op jjj , , 

IS POOR 
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pitches of the spirals. A distance is reached when convergence exceeds the speed of 
sound and the wave begins to shrink. 

From the point of view of an observer moving with the flow along o streamline, 
and looking at the wall represented by the interface, he sees the wall approach him— 
that is, convergence in his frame of reference. If we think of the wall os a piston 
moving into the flow, at the point where the sound wave begins to shrink, tl.« piston 
is moving faster than the local speed of sound, and a shock wave will form upstream 
from the piston. Thus, the origins of the shock waves will be approximately at points B 
in the slow stream and C in the fast stream marking where the sound wave stops 
expanding awoy from the interface. 

To find the opproximote locations of these points and tl. ir dependence on the 
speed o’lfference between the streams, we consider an idealized ccse in which the slow 
stream has constant speed V^, the fast stream has constant speed Vp the speed at the 
interface is = (Vj+V^)/2, and the Mach number, M, is constant throughout. Using 
the procedure given in Heinemonn and Siscoe (1974), we find the equations of the 
characteristics to be 


t?f = ^0 -0nrAf-f/f^) 


0 ) 


where rj is the azimuthal angle in the corotating reference frame, r is the heliocentric 
distance, and = V r^ = V^/Mn . Without the Bn terms, these are the equations 
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of Porker's solar wind spirals ( Porker ^ 1963, p. 138). The £n ferms represent the 
movement of the sound wove away from the spiral. The equation of the interface is 

with ~ V O . To estimate w4tere the diocks form we determine where the radial 

separation between the characteristics and the interface begins to decreose, i.e. set 

dAr/di) = 0. The result is shown in Table 1 for different speed differentials, Vf 

.1 

with V = 400 km sec and M = 4. 
o 

As expected, the bigger the differential, the nearer the sun the shocks form. 

Ill the biggest cose considered 160 km sec \ they form near the orbit of Mors. Any 
differential bigger than approximately 50 km sec ^ produces shocks inside the orbit 
of Jupiter. Although this example is idealized, it gives a fair test of the argument 
based on ch<vacteristics to predict qualitatively the essential geometrical aspects 
of the formation of corotating shocks. The prediction that typical solar wind streams, 
which have differentials bigger than 50 km sec \ should form shocks between the 
orbits of Earth and Jupiter is apparently confirmed by the Pioneer 10 and 11 observations. 
This justifies applying the argument to determine the out-of-the-ecliptic shape of 
cororcting shock waves. 

The application is straightforeward, and the result is immediate if we consider 
the situation at the poles. Here the border between the streams is a radial line — the 
polar axis. The pitches of the streams lines are essentially zero and, hence, so is the 
speed of convergence of the streams. A sound wave starting here will expand forever. 
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although it will slow down because of the radial decrease in sound speed. Ilius 
corotating shocks will not form over the poles. At intermediate latitudes, the pitchy 
of the streamlines and the border ore not zero, but they are less than at the equator. 

A sound wave must travel further before the convergence speed overtakes the sound 
speed and diodes form. Thta, the distance to the formation of the shodcs increases 
with latitude. 

The three dimensional geometry of the stream interface and the shock pair 
is sketched in Figure 3. The interface is genmtited by an expandir^ meridian circle 
that rotates about the polar axis os it expands so that its intersection with the equatorial 
plane moves along the border spiral. The leading points of the shocks are in the 
equoto’ial plane and the leading edges spiral outward as they move away from the 
equator, maintaining a proximity to the interface. A sketch on an expanded scale 
of a single corotating shock surface is shown in Figure 4. The motion of this surface 
is analogous to that of a tapered paper banner attached to a stidc Hiat is being twirled 
around the polo- axis. The length of the banner is determined by the lifetime of the 
solar wind stream. 

To study these structures an out-of-the-ecliptic mission is needed that covers 
the radiol range between Earth and Jupiter. Knowledge of the three dimensional 
nature of these shock waves is essential for the interpretation of cosmic ray data and 
for applications to other astrophysicai situations. The distortion of solar shock wav^ 
as described earlier and the highly structured geometry of corotating shocked waves 
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illustrate the complexity of the problems faced by galactic cosmic rays os they try to 
enter the inner solar system. Three dimensional probing of the interplanetary medium 
is required to obtain a complete picture of the extended stellar envelope of a 
representative from a major population of main sequence stars. The interaction of 
such stars with the interstellar medium in various galactic situations can they be 
treated with a fuller understanding of the stellar parameters. Comprehension of the 
three dimensional aspects of structures generated because a star rotates has application 
to contemporary ostrophysical problems such as the interaction of the Crds pulsar with 
the Crab nebula. 
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v.-v 

f s 

dAr/dt} = 0 

16 (kiVsec) 

12.5 (AU) 

32 

6.2 

48 

4.2 

64 

3.1 

80 

2.5 

96 

2.1 

112 

1.8 

128 

1.6 

144 

1.4 

160 

1.3 


Table 1. Approximate heolcentric distances in the equatorial plane to the 
formation of shock pairs dcj to interacting streams with various 
speed differentials. In this example the average speed is 400 km sec 
and the Mach number is 4 throughout. 
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Figure Captions 


Figure 1. A meridian plane cross section through a solar origin shock surface assumed 

to be circular at 20 solar radii and distorted into a quasi~ellipse at 200 solar 

radii by the action of equator to pole solar wind speed differential of 200 km 
-1 

sec . 

Figure 2. A sketch of geometrical features of the flow in the equatorial plane and in 
the corotating reference frame. The stream interface separates a slow stream 
(leading) and a fast stream (trailing). Sample streomlines ond corotating shocks 
in both streams are shown. The shocks deflect the streamlines parallel to the 
stream Interface— with no deflection they would intersect the interface. The 
circles are sequential snapshots of a sound wave starting at A and expanding 
while being convected with the flow (distortion due to the flow speed dif- 
ferential is neglected). Shocks form where the sound wave begins to shrink 
which happens when the speed of convergence of the streams toward the 
interface (due to the differences in the pitches of the various spirals) becomes 
bigger than the local sound speed. 

Figure 3. A sketch of the three dimensionol geometry in the northern hemisphere of the 
stream interface and associated shock pair. The interface is generated by an 
expanding meridian circle that rotates about the polar axis as it expands so 
that its intersection with the equatorial plane follows the interface spiral, like 
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Q needle in the groove of a record. The interface separates the shock po!r. 
Their leading point is in the equatorial plane and their leading edges spiral 
away from the sun as they move owoy from the equatorial plane, maintaining 
0 proximity to the interfoce. 

Figure 4. A sketch on on expanded scale of a single shock surface. The location of 
the equatorial plane is indicated by lines radiating from the sun. The motion 
of the surfoce is similar to that of a tappered paper banner attached at its 
point to a stick that is being twirled around the polar axis. The length of 
the banner is determined by the lifetime of the streams. 
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CHAPTER V 


SOLAR AND GALACTIC COSMIC RAYS 
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Abstract 

The reasons for studying cosmic-ray transport theory are 
summarized and the fundamentally three-dimensional nature of the 
process is pointed out. Observations in the ecliptic plane cannot 
unambiguously test transport theories since the solutions to the 
transport equations depend critically on boundary conditions and 
variation of parameters such as diffusion tensor out of the ecliptic. 
Sample calculations are shown which illustrate the problem. It is 
concluded that out-of-the-eci iptic observations are essential to 
further test transport theory. 
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I. Introduction 

The study of cosmic-ray transport theory has an old and vener- 
able history, for it is clear that in order to understand anything 
in cosmic-ray astrophysics we must first understand transport. The 
observed near-isotropy of cosmic rays was early recognized to imply 
that the orbits of the particles have been severely distorted in 
their motion through space. Since the work of Fermi (19**9) it has 
been understood that the motion must be treated statistically, since 
the plasmas and magnetic fields through which the cosmic rays move 
are irregular and turbulent. This basic fact leads to the view that 
the spatial motion of cosmic rays is to a very good first approxima- 
tion a random walk in three dimensions. 

The solar wind provides an excellent local laboratory for testing 
our ideas of cosmic-ray transport by comparing In situ observations 
with theory. A comprehensive theory has been developed which has 
had reasonable success in explaining the various cosmic-ray phenomena 
in the solar wind (see reviews by Jokipii, 1971 and Volk, 1975). 

Before going on to discuss this theory In detail, it is useful to 
emphasize that irrespective of the detailed theory, cosmic-ray trans- 
port theory is fundamentally three-dimensional and that this problem 
is much more severe than for some other aspects of the physics of 
the interplanetary medium. The general picture is illustrated in 
figure 1 for particles originating in the galaxy or at the sun. It 
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■» avtden; for example, there Is no^ viev/Ing direction in vvhich 

one sees cosmic rays which have sampled only the ecliptic piane. 

This conclusion depends only on the turbulent nature of the inter- 
planetary plasma and the consequent stochastic motion of the cosmic- 
ray particles. 

One can, of course, dream up experiments which could be used to 
study cosmic-ray transport locally. For example, a useful experi- 
ment would be to inject a small beam of "tagged" cosmic rays at one 
point and measure them a distance -0.1 a.u. away. But such experi- 
ments are clearly not possible in the near future. We must be 
content with v.’nat nature provides, and work with the full three- 
dimensional problem. 

The preceding discussion has intentionally been as general as 
possible in order to emphasize the fundamental nature 
0 " the conclusions. In the next section, the current detailed theory 
of cosmic-ray transport is discussed, then some quantitative calcu- 
lations are presented v;hich illustrate the effects of uncertainties 
in parameter outside of the ecliptic. 
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II. Current Transport Theory 


The current theory of cosm{c**ray transport has been comprehen* 
stvely discussed in an earlier review by the author (Jokipii, 1970 • 
The theory works at two levels, which could be termed the "macro- 
scopic" and the "microscopic". This is illustrated schematically 
in figure 2 which shows the particle being scattered randomly by the 
magnetic fluctuations. Hence the particle pitch angle 9 undergoes a 
random walk. The details of this scattering process and its relation 
to the detailed microstructure of the magnetic fluctuations are 
studied in the microscopic theory, whereas the resulting spatial 
diffusion is considered in macroscopic theory. Thus, for example, 
in the simplest form of quasi linear theory, where only planar magnetic 
fluctuations with wave vector parallel to the average field are con- 
sidered, the pitch-angle scattering is characterized by the Fokker- 
Planck coefficient 


<&Q^> „ p (I, 

1u|wB„2 ^ 



( 1 ) 


v/here ■ qB^/mc is the particle cyclotron frequency in the average 
field B-, y “ cos0, w •» particle speed, r ■ w/u,, and P, (k) is the 
spatial power spectrum of the magnetic fluctuations as a function of 
wavenumber k (see, e.g. , Jokipii 1966, 1971). Other more complex 
expressions result if other magnetic fluctuation configurations are 
assumed. Similar expressions result for the other Fokker-PIanck 
coefficients such as <Ax^>/At, etc. Hence, if x is a direction 
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normal to we have 


<Ax^> 

At 


-^Px(k-O) . 
B 2 


( 2 ) 


The above expressions are examples of microscopic transport theory. 

However, in most situations the cosmic**ray angular distribution 

is driven essentially isotropic by the scattering, so that the 

transport must be described by the diffusion approximation. If 

is the particle density averaged over pitch angle as a 

function of position r, time t ^nd energy T, then the diffusion 

equation reads In the rest frame (Parker 1965, Gieeson and Axford 1967, 
Jokipii and Parker 1970) 


^ ^ r JiU 

at “ 3x. *-'^ij 3Xj 


1 


3 / 
aT 


I u) 


(3a) 


with an associated flux of particles in the rest frame 


F 
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-K ^ + 

Ij 


V 

w. 


.U 

I 


V . 
-SUI 
3 



(3b) 


where Kjj is the cosmic-ray diffusion tensor, is the solar wind 
velocity and o(T) «• (2mc^ + T)/(mc^ + T). The associated anisotropy 
is |6| “ 3|f(/ivU. In a coordinate syster:? with the z-di recti on oriented 
along the average magnetic field k.j may be written 
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where the parallel diffusion coefficient k,, , the perpendicular 
diffusion coefficient Kj^and the antisymmetric diffusion coefficient 
can be written 



The diffusion equation (3) embodies the macroscopic theory which is 
connected to the microscopic theory by equations (5) ~ (7)> At 
present, comparison of transport theory with observation must be 
done through the use of macroscopic theory and then working back to 
the microscopic theory. We are unable to measure <A^/At, etc. 
di rectly. 

With regard to the problem of out-of-the-ecl iptic exploration, 
it is clear that the solution to equation (3) depends on a knowledge 
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of Kjj and throughout the modulating region and on the value of 
U on the boundary. There appear to be no purely locally measurable 
properties of the solution which can be used in checking transport 
theory. 
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III. IMustratlve CalcuUtions 


To Mlustratia the rather large uncertainties Introduced by lack 
of knowledge of parameters out of the ecliptic, I present in this 
section a summary of some analytical calculations published elsewhere 
(Owens and Joktpll, 1971) • This problem has also been considered by 
Sarabhat and Subramantan (1966) and LlettI and Quenby ( 1968 ). Con- 
sider th'i ll^year solar cycle modulation of galactic cosmic rays by 
a solar wind which is not spherically symmetric. The 11-year solar 
cycle variations are usual regarded as slow enough that the time 
derivative In equation (‘S) can be neglected in solving for U. In 
this calculation the various parameters In the transport equation 
are allowed to vary with heliographic latitude and the effects of 
the variation on the density U and the flux F in the solar equatorial 
plane are considered. 

The following forms for the various parameters were assumed: 

The various components of the diffusion tensor are independent of 
energy T and are proportional to heliocentric radius r out to some 
boundary r ■ D where ic.j -*■ «. The cosmic-ray density U(r,T) is 
assumed i:* take on a given interstellar value U„(T) - AT“2*5 at r D. 

Finally, we follow the us'.:^! practice of circumventing the 
problem of proper boundary conditions at the Sun by requiring the 
solution for U to be finite at the origin. This will be shown to 
lead to negligible error in the present case. See Jokipii (1971) 
for a more complete discussion of this problem. 
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Equation (3a) for U(r,0,^,T) becomes» in spherical polar 
coorJi nates 


J-l- fr2(v ti - kh ~)1 - !_iL 

.2 I ** *'■ j s*"® 


kinaf-iilL . -3 

P ®'r ae r sine 


— , 
aUv 


JL 

a ' • 


2V 


au ^ ^ au* . _J1 1. /Tu i 
r sine a^ 'r sine a^ r ae^ ~ 3r 3T ' 


0. 


(8) 


In what follows, we consider only mocels which are axisymnetric 
(independent of^). lf0»>Ois the axis of solar rotation, this 
corresponds to a latitude-dependent solar wind. Effects of magnetic 
sectors could be represented by taking 0=0 along a sector. The 
two terms containing cancel from the density equation (8) if 
3<^/a3 = 3ic^/3$ = 0, as in our axisymmetric model. Thus the terms 
in :<^ do not contribute, and equation (8) becomes 


1_ ^ 
r^ 3r 



1 

r^sinO 



2V a 
w 3 

3r 3T 


It (TUa) 


0 . 


(9) 


Since K|, and are taken to be independent of energy, equation (9) 
separates, and one writes 


U(r,Q,T) = S(r,0)5(T). 


( 10 ) 
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Upon neglecting the sinall energy dependence of a(T), one finds 
inmedlatety that S must satisfy 4- 


r2 





( 11 ) 


with 

«(T) - T(3q/2a)-l . 


( 12 ) 


The separation constant q must be chosen to agree with 
U^(T) « in which case we choose q » -1. 

A number of different forms for the variation of the parameters 
V^, Kjj and %*ere chosen. Define 

V^(0,r) “ V^[l+6(e)J, Kj|(0,r) - K^[l+e(0)]r 

and » pic^ . (13) 

The problem may be solved in terms of a series of Legendre 
polynomials and the general solution is given by Oweis and Jokipii 
(1971)-. An illustrative velocity variation is given by 
V *» V^[l + .30P2 (cosG)] where is the average velocity and P2 is 

the Legendre polynomial of order 2. This velocity variation is 
illustrated in figure 3- Some typical results of the calculations 
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are illustrated iu figures 4, 5, and o\ • It is clear that many of 
the parameters observed in the solar equatorial plane can be sub- 
stantially changed by varying the solar-wind parameters out of the 
solar ecliptic plane. Of particular note is the fact that the radial 
anisotropy observed in the equatorial plane is extremely sensitive 
to parameters outside the equatorial plane. As shown in figure 6, 
even the sign of the anisotropy can be changed by relatively small 
velocity variations. 
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IV. Conclusions 

may conclude from the above discussion that fundamental 
ambiguities In testing cosrolc-ray transport can be removed by 
carrying out measurements out of the ecliptic. The out-of^the-ecliptic 
a»asureroents most necessary to study cosmic-ray propagation are: 


a) Heasurements of the flux of cosmic rays as a function 
of solid angle and energy with as much resolution as 
possible. The anisotropy measurements are easier to 
carry out on a spinning spacecraft. 

b) Simultaneous measurements on board the same spacecraft 
of the plasma and magnetic field. These measurements 
should be spaced In time so that good time-series 
analyses (power spectra, etc.) can be obtained. 

c) The spacecraft should go as far out of the ecliptic 
as possible to insure that any variations with helio- 
graphic latitude will be seen. 

d) Good base-line measurements at Earth or in the eclip- 
tic plane should be obtained simultaneously. 


It is not crucial that these measurements be carried out at 
constant heliocentric radius, although this might aid interpretation. 
It appears that measuring protons and possibly electrons will be 
adequate and it is better to optimize measurements for one species 
rather than compromise these in order to study composition. 
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Fig. 1. Schematic illustration in a meridian plane of typical 
trajectories of solar cosmic-rays (dashed line) and 
galactic cosmic rays (solid line) which reach a detector 
at 1 a.u. The fundamentally 3-dimensional character of 
the motion is apparent. 

Fig. 2. Illustration of the ''scattering" of a cosmic-ray particle 
by a magnetic irregularity. If the region over which A0 
is correlated is small, a Fokker-Planck equation results. 

Fig. 3* Illustration of the velocity profile used in some of the 
calculations. 0 is the angle relative to the z-axis. 

Fig. 4. Comparison of linearized and exact solutions for the case 
V = 0*30 P^(cos 0)], P = 3* The density 

U(r,9) is given for both solutions as a function of angle 0 
r = 0.2 D in terms of = U(r = 0,0). Straight line shows 
the corresponding result for an isotropic wind V = V^. 

The linearized calculation keeps terms only to first order 

in 0 as defined in equation {k) . From Owens and Jokipii (1971). 

Fig. 5. Total cosmic-ray flux as a function of r,0, and with Fa 

9 

suppressed. The pattern is azimuthal ly symmetric and even 

about the equator, (a) The parameters V(o) = Vp[l + 0.30 p 2 (cos ©)], 

Kj. (r,3) = x^r and p = i were used. The position of Earth 

for 0=5 a.u. is indicated. Of particular interest is the 
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virtual source of particles in the equatorial plane at 
r * 0.5 0. (b) Same as (a) except that V(0) = + 0.46 P^Ccos ©)]. 

Note that the flux is much larger and virtual source has 
moved out to r - 0.9 0. From Owens and Jokipii (1971). 

Fig. 6. Anisotropy in the equatorial plane at r ■ 0.2 0 for 
v(e) = + V^Pj,(cos O) , as a function of V^/V^. 

Note that in our model the anisotropy is radial. From 
Owens and Jokipii (1971). 
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COSMIC RAY MOOUUTION IN THREE DIMENSIONS 


N76-241S1 

J. J. Ouenby 
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Abstract . A brief critique of spherically symmetric conventional 
modulation theory is supplied. Estimates are made of the cosmic 
ray intensity at high solar latitudes. Direct evidence for 
significant off-ecliptic cosmic ray gradients is reviewed in 
support of the requirement for an off-ecliptic spacecraft mission. 
The possibility of measi'ring the galactic spectrum is discussed. 


Spherically symmetric modulation theory and its problems . 

Cosmic ray modulation arises because motion of the energetic particles 
along the spiral interplanetary magnetic field lines is controlled by 
scattering due to magnetic irregularities which are being convected 
outwards by the solar wind. Inward diffusion is balanced both by outward 
convection and energy loss of the particles as they suffer adiabatic 
deceleration in the expanding wind. A Fokker-Planck equation expressing 
these effects in a spherically symmetric, steady-state situation is 
(Parker 1965, Gleeson and Axford 196?) 

7 ^ If 1 - I ; ^ <»™) (') 

Here U is the differential particle number density at position r and 

T + 2T 

kinetic energy T,a= with as rest mass energy, V is the solar 

I + I o 

wind velocity and K,., the effective radial diffusion coefficient, is given 
by Kf ■= Kji cos^^» + Kj^sin^iJ» for K ar.d K , respectively the parallel and 
perpendicular diffusion coefficients with jJ* ■= cos”^ / I®.! l£.l 

Particle streaming ^ is due to diffusion in the wind frame plus an 
additional term involving a Lorentz transformation to the rest frame via 
the Con^ton-Getting factor C such that 
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S - cyu - K.grad U 


( 2 ) 


with C ^ ^ (oTu) and K as the t^sor diffusion coefficient. 

At magnetic rigidities exceeding about 1 GV, radial streaming is 
usually negligible and under these circumstances the Fokker-Planck 
(1) can be integrated to give a relatiwi between the Galactic <tensity 
U(rg, Tg) and the observed density U(r,T), 

U(r. T) - U(rg, Tg) exp (- z/® dr) (3) 

for spherical symmetry with determined by the effective outer bound 
to the interplanetary scattering process. 

It is important for Astrophysics to know the energy spectra of 
various nuclear species of cosmic rays in the galaxy and a common method 
for achieving this "demodulation'* is as follows: 

(a) Estimate the Galactic electron spectrum from radio synchrotrcxt data. 

(b) Compare the near-earth electron spectrun with the Galactic spectrum 

to find the magnitude and rigidity dependence of (^/Ky.).dr. 

(c) Use the results of (b) to correct the near-earth proton and heavy 
nuclei spectra for modulation. 

Various objections can be raised to this scheme. First, the average 
Galactic electron spectrum derived from radio data may not represent the 
local spectrum and in any case one must be sure of the model enabling the 
effects of local, cold, absorbing interstellar clouds (Goldstein et al . 
1970a) to be taken into account. Second, there is no adequate theoretical 
explanation as yet for the magnitude, rigidity dependence and radial 
dependence of K^, as witnessed by investigations of solar proton diffusion 
'e.g. Webb et al. 1973 ) and the measurement of cosmic ray radial gradients 
(<I0%/AU) which are much less than those expected on the basis of 
theoretical values (e.g. Webber and Lezniak, 1973)- It is important to 
know the r dependence of the integrand in (3) since adiabatic deceleration 
is inversely proportional to r. This importance becomes apparent in our 
third point, based on the work of Goldstein et al. (1970b) and Gleeson and 



Urch (1970 «4K>t for certain choices of r dependence and reasonable 
magnitudes of Kr» find that adiabatic deceleration may completely 
exclude all <100 HeV/nucleon particles from the inner solar system. 

Thus no certain knowledge of the lowest energy Galactic primaries is 
possible. Fourth, the cause of time variations in the modulation is 
unknown, with not enough solar cycle variation in the solar wind 
paraaeters being observed near earth to account for the known changes in 
the cosmic ray flux (e.g. Hedgecock et al. 1972). For example, the 
power spectrum of magnetic irregularities at 10'*^ Hz should change by a 
factor >2 in order to account for the observed modulation of 1 GV particles 
between 1965 and 1969 according' to rescxtant, wave/particle interaction 
theory while in practice the change is <10« (Hedgecock 1975)* Off- 
ecliptic control of modulation via the effects far beyond 1 AU of solar 
streams emerging from the zones of maxiimim solar activity may be the 
only way to explain the cosmic ray 11-year cycle. Hence careful study 
by off-ecliptic spacecraft is required. 


2. The off-ecliptic route to the boundary . 


It has been thought that a direct determination of the Galactic cosmic 
ray charge and energy spectrum can be achieved employing ecliptic plane 
spacecraft trajectories to the outer planets. In this way the problems 
of section 1 are all by-passed. However, the boundary to modulation 
could be as far as 100 AU and the low, measured cosmic ray density 
gradients seen by Pioneer 10 render this approach uncertain. Alternatively, 
we note that cosmic rays have an easier inward motion over the solar poles 
where the geometric path length along the interplanetary field lines is 
much shorter than in the tightly wound spiral regime of the equatorial 
plane (Lietti and Quenby 1968). An appropriate Fokker-Pl anck equation 
for the steady state which relinquishes the requirement of spherical 
symmetry and takes into account the spiral geometry is then 


1 A. 

A 6s 


( ak 


>1 6 s 


(r^evu) 


if “ 0 for path length s along a magnetic flux tube of area A. 
Allowing — K^r sin^9 where 0 is solar latitude yields 
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(5) 



CVr 


(Vrg) 




C\fr *sin2-<l 6 

* ‘ - in.r- ' 


and if q ■■ o and ^/Ug “ Ve at I GV, we obtain the following table for 
the percentage residual modulation at lAU: 


e “ 90® 0-30® 

'^/Ug (at 1 GV) 63% 22% 


0 “ 0 ® 


9% 


Thus on a sin 4 >le model for scattering, a spacecraft passing over the poles 
at 1 All may see ^90% of the unmodulated intensity and therefore get a better 
measurement of galactic conditions than is available at Jovian distances. 


3. Direct evidence for off-ecliptic gradients 


Observations confined to the ecliptic plane can only reveal the existence 
of off-ecitptsc effects by noting modifications to cosmic ray streaming 
from the expectations of the spherically syimnetric model case. Equation 
(2) in a more general form is 

where v is energetic particle velocity and m is cyclotron frequency. It 
has been assumed that direct slippage of particles across field lines makes 
only a small streaming contribution ''' t being time to 

travel one parallel mean free path). Furthermore, we assume the short- 
circuiting by scattering of that perpendicular gradient which would cancel 
out the anisotropy due to (E^x drift in the non -scattering limit when 
Liouvi lie's theorem applies to the particle intensity. At high rigidities, 
some few to a hundred GV, the radial streaming is negligible over long 
periods with the third term on the right of (b) cancelling out on average. 
Then the first two terms combine to give the streaming from the east or 
1800 hr LT anisotropy. Vfhen, however, the anisotropy is studied in 
practice as a function of sign of the interplanetary field sector structure, 
two effects of the third term become apparent. A north-south anisotropy 

s II 

arises due to ( -r— ) , x B, or the effect of the radial gradient and an 

ecliptic plane anisotropy arises due to ( x or the effect of off- 
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ecliptic gradient. Hashio and Bercovitch (1972) find Gj ~ 5.5 R~^*®%/AU 
directed north south In 1967/68 for the latter effect, possibly 
physically resulting from the exoess northern hemisphere solar activity 
suggested at that tima by coronal green line 5303A emission. 

The previous discussion refers to the first derivative of density, 
but studies of the second harmonic of the cosmic ray intensity can reveal 
the presence of a rising or falling, symmetric, off-ecliptic gradient via a 
dependence on the second derivative. Lietti and Uuenby (1968) essentiaily 
use a version of (5) for the rising gradient case to predict a second 
harmonic^wi th direction of maximum perpendicular to ^ and amplitude 
a2 " J ■J07 '^6*6®5P % at rigidity P, cyclotron radius p. This 
expression is in reasonable accord with o&servations although Nagashima 
et at. (1972) claim that a cylindrical pitch angle particle distribution 
about with a different physical cause better fits cosmic ray anisotropy 
data. 


With finite K , the symmetrical gradient will either feed particles 

Jm 

into the equatorial plane or draw them off to higher latitudes, thus 
setting up radial streaming. A correction to the Fokker-Planck (1) Is 
employed by adding ^ 

div ( • ^ ) 1 to the right hand side. Dyer 

et a1. ( 1974 ) in particular evoke this streaming for a falling gradient to 
explain a sunward flow at > 1 GV seen by a satellite detector which is 
too large to be explained by any energy toss effects in a spherical ly 
symmetric model. These authors require maximum modulation over the 
sunspot zones with mart di ana I flow patterns set up to draw particles down 
from 0 and up from 0 ==90®. Cecchini et al. (1975) have developed 
a computational model to confirm the above model. Chief features are: 


K - K^0P f(0) I 


''"o/expd) 

(Vro)^ 

C/ro) 


r < 1 AU 
r > 1 AU 

r < 1 AU 
r > 1 AU 


f(e) =• 1 + cos 8 (-3 + 5 cos^ 0 ) ; 




e “ 0.5» 1AU; V„ — const.; Pq 


“ 100 MV; Ug oi(T + T0)'2‘75 

Thus ™ 0.05 at 1 GV. The results of employing an alternating gradient 

technique to solve the corrected Fokker-Planck, with finite 1^, is to predict 
an inward streaming *^0.3% in amplitude between 2 and 10~^ GeV, a radial 
gradient 'v10%/AU at 1.1 GeV between 1 and 10 AU and a ratio . IL) 

at 1 AU varying from 2.5 at 0 “ 0 to 0.7 at 8 = 60° for 1,1 GeV protons? 

Hence it is possible to explain the radial streaming with reasonable gradients 
and Ki/„ ratios. 

^ ii 

4. Conclusions . 

We have shown that study of cosmic ray modulation by integrating the transport 
equation outwards in the ecliptic plane, assuming spherical symmetry, 
encounters various problems. The transport processes and boundary conditions 
are insufficiently well understood, modulation may be controlled by off- 
ecliptic gradients and assymetries can have noticeable effects on solar 
equatorial plane observations. Three-dimensional study of the solar 
cavity cosmic ray distribution is required to: 

(a) Measure off-ecliptic gradients and streaming. 

(b) Enhance understanding of the solar control of intensity time 
variations. 

(c) Gain better knowledge of boundary conditions, especially the 
possibility of measuring a near-Galact ic energy spectrum over the 
solar poles. 

Objectives (a) and (c) are satisfied by a Jovian swing-by mission 
but (d, requires a direct injection at 1 AU spacecraft for detailed time 
variation studies on solar wind and solar parameters. 

reproducibility op the 

OE: . AT. PAGE IS POOR!' 
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Abstract 


Based on a modified NKB analysis of the Interplanetary Irregularity spectra, 
a discussion of the radial dependence of the radial cosmic-ray diffusion coef- 
ficient at polar heliographic latitudes is presented. At 1 AU radial distance 
the parameters are tedcen to equal those observed in the ecliptic. In the sense 
of. a present best estimate it is argued that relativistic nuclei should have signi- 
ficantly easier access to 1 AU at the pole than in the eci..£»tic. The reverse 
may very well be true for the direct access of very lt>w rigidity particles. 





Introduction 


The access of galactic cosmic raya to the inner solar system is regulated 
by the cumulative effect of the irregular fields which scatter the par- 
ticles on their way in. To determine the degree of this access i . helio- 
graphic latitudes that are significantly different from the solar equa- 
torial plane# one has to rely* on theoretical estimates. These conventional- 
ly assume a diffusive propagation scheme. Of primary importance is then 
the spatial diffusion coefficient or, more generally# the diffusion ten- 
sor #and in particular its spatial variation in the solar system. 

In the learlier work of V61k et al « (1974) the irregular fields were assumed to 
be Alfvfen waves of solar origin# convected outwards by the solar wind. 

In a (on the average) stationary interplanetary medium with axial symme- 
try around the solar rotation axis# this led to a radial direction for 
almost all wave normals beyond about 1 AU. Assuming the effective average 
magnetic field to be in the ideal soiral field direction and the wave am- 
plitudes to vary according to the (WKB) approximation of geometrical 
optics# a radial dependence of the coefficient for diffusion along the 
average field was calculated. As a sensitive function of the angle be- 
tween wave normals and average field, the value of this diffusion coef- 
ficient, or equivalently, of the mean free path, varies from a minimum 
at O degrees to infinity at 90 degrees. Since in the solar equatorial 
plane the angle between the radial and the spiral field direction is 
large already at 1 AV, and increases with radial distance, deviations of 
the wave normals from the radial direction have been discussed by Richte r 
(1974) in the context of the solar wind stream structure; see also Hollwcg 


(1975) . 



do not Intend to evaluate here the effects of deviations from the 


simple picture of Vdlk et al . (1974) at low latitudes and refer to a 
futvire publication ( Morfill et al«, in preparation) . Ka shall rather 
concentrate on the situation at polar heliographic latitudes. There the 
assumption of all wave vectors )c being radial and parallel to the average 
field ^ leads to the minimum value of the radial diffusion coeffi- 
cient K> as far as its dependence on the angle between ]c and B is 
concerned. Thusf within the WKB approximation for the radial development 
of the scattering centers « this provides a lower limit for Ai . Such a 
lower limit is of interest^ if- it can be argued, or at least be specu- 
d in a reasonable way chat even in this case there is essentially un- 
impeded access to eibout 1 AU for a significantly larger part of the galac- 
tic energy spectrum than in the solar equatorial plane. 

We shall give a short discussion here, motivated by three considerations. 

The first one is that in two earlier publications this author was associated 
with ( V61k et al . , 1974; V61k , 1975), a rather different result was be- 
lieved to hold true. Another reason is given by the present discussion of 
an ex-ccliptical probe to explore the sun and the interplanetary medium. 

The final reason is that via easy access at the poles much larger regions 
of interplanetary space might possibly be populated by particles of galac- 
tic origin. 

In the next section we shall present che general behavior of K. in a modified 
WKB analysis, using simple approximations to two rather different measured 
interplanetary power spectra. In the last section the relation to the expected 
actual situation at the heliographic pole is discussed. 
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2. WKB Analysis 


Consider the singles t« axlsyometrlc, interplanetary neditUDf where all quanti- 
ties depend on heliocentric distance r and where only the ideal spiral field 
B with components 








V*o 


V ’ 


and the power spectrum P (f,r>9) may in addition depend on heliographic la- 
titude Q . Here r is a radial reference level# B is independent of Q ; 

o to 

^ denotes heliographic longitude,0^c;2.6S x 10 ^ cps is the emgular fre- 
quency of the svin's rotation# V is the solar wind speed (assumed to be radial 
and constant) # and f is tlie wave frequency seen by an observer at test re- 
lative to the sun's center. Then it is simple to show ( V61k et al . , 1974) 
that for Q — ► ^ the radial diffusion coefficient is given by 




1 



In equation (2) we have ® ^ ratio of particle velocity w to the 

velocity of light c; R denotes particle rigidity; v^ “ B • (4X^ )”^ is 
the (vectorial) Alfv&n - speed# where ^ is the average solar wind mass den- 
sity;^ • is the cosine of the particle's pitch angle# where w^ is 

the velocity component parallel to B; r^ ■ 1 AU, The ^^1 i f ic^t^t^n ^^ ^c^^ 

(jWGIaAL page IB ^OOR 
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for th« Alfv6n wave aaplitudes 


( 3 ) 




V + '^AV* (^) 


li^ (t) I 

W+V^^(v) / '*M 


t yfi^M (Jk)^ 

1 


is independent o£ wavevector and v7 • If both and r are larcje compared 

to the AlfvAnlc critical radius (t$ 20 in the equatorial plane ) , then 

the approximation on the far r.h.s. of equation (3} holds in this lowest 

order WKB approximation . Kov/ever> following equations (1) and (3) we have 

(at 0 ■ 90°); (r ,©)^ / B^(r,0) ^ r, where is the total 

2 2 

power in the fluctuations. Thus, the possibility arises that <Sb^>/ 
beyond some radial distance in which case we expect wave propagation to be 
rather drastically altered by nonlineai. ofCects that ultimately should lead 
to dissipation oi wave energy. To take this possible effect into account, we 
modify equation (3) so that for all r 

(3a) ^ ^ B (•<") 


A similar device has been used by Hollveg (1973) and appears as a simple if 
crude way to take the inadequacy of the linearised WKB approximation into 

account. 


It is clear physically that equations (1), (2), (3) and (3a) also hold :f, 

at given r #^ro' S different from their values at ar, long a- 

their dependence on © is weak. enough, such that for example i •.n' 1. 

V dG 

For the rest of this section, however, we will consider (r^), V, ^ (r^) 
and P(f,©,rj) as given by their values in the equatorial piano, 
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Observed poeer ^ectra at Qs&O and (e.g. Joktpii and Coleaan » 1968; 

Bercovttcht 197t)f generally approxinate rather well a power law for P(f«r|) 
at high frequencies* while flattening at low frequencies* Qualitatively* an 
analytical form 


( 4 ) 


P (h = (^f 




With parameters f^^and q <0|^s not an unreasonable representation. With 
equation (4) the general character of 1C (r) can be inferred. For small enough 
rigidity R, so that [b^| • V/R> 27rf^ we have lc r^^^^ whereas ft^r^ for 
|Bj* v/R^ 2TCf^. The exponent of r is additively increased by unity ^en- 

ever the suppl^entary restriction (3a) comes into force. For our present 

3 2 

discussion we choose q = - /2. Thus# /C co.ist (or ^r) for small R • r # 

3 4 2 

whereas K A#r (or ) for large R • r . At suf fi^i^*itly large r# K will 

4 4 

be ^r (or Air ) for any R. At fixed values of q and R, the transition of 
3 4 

IC to the r (or r ) dependence increases with decreasing f jCf . equation (2) • 
As an aside we should mention that for this value of q and for the present 
case of k B B the quasilinear expression for IC used in writing in equation 
. (2) is numerically quite satisfactory and the modifications in the region 
around » O (e.g. Jones et al .# 1973) therefore not essential. 


Uumerical results are shown in Figures 1 and 2. Figure 1 is calculated using 


the spectrum of Jokipii and Coleman (I960) where, in reference to equation 
(4)# RJ 16 X 10 ^ * ^^(Hz) ^ f^<^7 x 10 ^ and Alfven wave 

propagation was started at 20 Rq (solar radii) • Figure 2 uses the spec- 

trum published by Bercovitch (1971) which we roughly approximate by C x 10 

2 1 -6 3 ^ 

^'(Hz)^^”^; f Q 2 X 10 *Hz; qsJ • /2. Although the results are not very much 


-3 
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different for both cases near r « 1 m, the very different f^-values lead to 

3 4 

strong differences in ths cmset distance of tte r (or rather r ) law* We 
^lould mention here that we used here tte ooe^onent of tte magnetic spectral 
tensor perpendicular to the ecliptic to represent tte total power per fre* 
quency interval* For true axisyoaetry of the spectrum all 1C values in Figures 
1 and2 should be multiplied by a factor 1/2* This is an unavoidable uncertain* 
ty 

The Interesting aspect of these results is that they inq>Iy little modulation for 
relativistic nuclei, where adiabatic deceleration is small^ considering the 
value of f^ in Figure 2 as a rather extr^ae lower limit to the actual situation. 
For 1 GeV protons, for exaii^le, a 10-20 percent modulation is estimated in 
the diffusion convection approximation. 

3. Discussion 

The above results were obtained by fitting^ at l.AU^the average solar wind 
parameters as well as power spectra by the corresponding quantities observed 
at 1 AU in the ecliptic; the spatial dependence of the spectra assumed a modi- 
fied WKB approximation for the irregularities. In reality, the polar region of 
the corona may well be a large, stationary coronal hole, resulting in a some- 
what (perhaps fifty percent) larger flow speed and, possibly, somewhat more 
power in the frequency region f All this would lead to a moderately in- 

creased modulation. 

On the other hand, it may very well be that the power 

in frequencies is much smaller at ^ = lC/2 than near 0 = O. We have taken 

these fluctuations also to be Alfvenic which leads to the amplitude variation 
given in equation (3) .In reality, the part in the spectrum with may well 

be due to the solar wind stream structure ( Goldstein and Siscoe , 1972) . If 
the latter is assumed to bo absent at 0 - 90^, also the power at f would 

be absent with a corresponding decrease in modulation. In this light, also 



tile poesiblll^ of increased scattering at larger distances due to local pro- 


dection of eaves - a situation that is quite likely at > appears 

zatlier weak. Irregularities produced l>y enhanced (coa^ared to the ecliptic) 
tntinilence due to radially increasing departures from thermal equilibrium at 
tiie poles should be assumed to have small scales, irrelevant for cosmic rays 
even in tbe case of the Firehose instability. It should be added that in con- 
trast to a popular feeling this result for 1C and the consequent argument for 
modulation has little to do with the shorter geometrical path along B of a 
galactic particle to. say. 1 AO. but rather to the strong decrease with r of the 
magnetic field at the poles. 

yhos. barring unknown new effects, the present best estimate is that cosmic 
ray access at the poles should be significantly better than at rela- 

tivistic nuclei. For very low-rigidity particles on the other hand, the sharp 
increase of fC with r occvirs only at such a large radial distance that their 
direct access may be at least as strongly prohibited as in the equatorial plane. 
However, for this last kind of statanent. the present estimate is not well 
suited. 
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Figure It Ihe radial diffusion coefficient K, at heliographic latitude 

9 B 90** as a function of radial distance r in solar radii 

for various proton energies. The power spectrtim P(f)=C^ • f*^ 

with the values = 16 x 10 ^ f^= 7 x 10 ^ Hz 

(Joklpii and Coleman , 1968) . Wave normals )c are assumed to be 

radial. The calculation was started at r » 20 . 

o O 

Figure 2: The same as Figure 1 with values of and adapted to the spec- 

tra of Bercovitch (1971) . 
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Introduction 

Observation of solar cosmic ray fluxes at high heliocentric latitudes 
would provide several new dimensions in specifying important physical 
processes in the solar atmosphere and in ter plane taty space. Valuable 
new approaches would be available for such problems as the steafy state 
acceleration in the solar atmosphere* propagation of fast charged particles 
in the solar coronal magnetic fields and in interplanetary space, the true 
interplanetary particle spectra free of transient solar influence, and accel- 
eration in interplanetary space by shock waves. For this brief review 
of what might result from a program of solar cosmic ray observations on 
"out-of-the-ecliptic" spacecrait the following outline will be used: 

A. Th'' magnetic fields of the Snn at high latitudes 

B. Propagation of fast charged particles in the solar corona and 
in interplanetary snace at high latitudes 

C. Origin of interplanetary particle populations 

D. Other particle phenomena in interplanetary space(e.g. , accel- 
eration of shock waves) 

E. Effect of spacecraft mission characteristics on solar cosmic 
ray studies at high latitudes 
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A. Solar Coronal Magnetic Fields 

It is now well known that the magnetic fields close to the Sun 
control the intensities of charged particles that appear in interplanetary 
space. Th^ do so in a variety of ways and to vazying degrees. For 
example, particles from flares near the East lim ) are detected at central 
meridian and more westerly longitudes only several hours after the flare 
(Van Hollenbeke et al . , 1975). The measured low cross field diffusivily 
of the particles in interplanetary space requires that the diffusion occurs 
in* the corona. This is a reasonable requirement since hydromagnetic 
wave activi^ should be more intense in the solar corona and the wavelength 
of such waves can be of the size needed to effectively scatter solar flare 
particles. 

Energetic solar particles stream from the solar corona in high 
intensities for many days, even weeks, following large solar flares. In 
general terms, this basic observation means that there is a sequence in 
which the particles are first accelerated, then stored for some time in 
the solar corona, then released into interplanetaiy space. The quantitative 
behavior of each of these steps is not well known and in fact is at the center 
of much controversy. For example, the acceleration may be essentially 
steady state and the storage transient. On the other hand there is the 
possibility that the acceleration is impulsive and the storage long term. 
Such a situation requires very low collisional energy losses, something 
that might be achieved in a "cosmic ray plasma" where all particles arc 
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energetic, approaching a Maxwellian distribution (Anderson, 1972). In 
either case the magnetic fields of the sun play a major role in each step 
of the three step sequence. All our experience with these processes has 
been in the magnetic fields of the solar activity zones. It seems unlikely 
that the full sequence df acceleration, storage and release takes place 
in the high latitude coronal fields. However, there are two inversely 
related questions of basic importance: 

1. Do solar particles, accelerated in the activity zones, propagate 

into the high latitude solar regions? 

2. Can we use the solar particles, if th^ do indeed move into 
the high latitude regions, as probes of the solar magnetic 
fields there? 

With these two questions as the basic motivation for solar cosmic 
ray studies at high heliographic latitudes, we next summarize some of 
the information now available on high latitude solar magnetic fields. 

Figure 1 is a photograph of the solar corona that shows the polar plumes. 
This old photograph is of interest because it shows an intermediate corona, 
the kind expected in the 1983 to 1985 period. The polar plumes have 
generally been considered to outline the magnetic fields over the Sun's 
polar cap and their appearance suggests that the field lines in these 
high latitude regions are open . Notice that the polar plumes appear only 
above latitudes of 50 to 60° . Because these field lines are so different 
from the activity zone and transition zone field lines, a strong impetus 





is given to going to the highest possible solar latitudes in order to encounter 
a qualitatively different phenomenologjr. 

Magnetograph observations suggest that the S^:n's polar magnetic 
fields are not uniform bundles of lines. Figure 2 shows measurements 
by Stenflo (1971) that lead to the following conclusions about the high altitude 
fields: 

1. They are patcl^ with large regions (channels) of weak (less 
than 10 Gauss) fields. These channels extend from the poles 
down to latitudes less than 50^ . 

2. The strengths in the strong field patches range from 10 to 20 
Gauss. 

Polarization methods of magnetic field determination (Charvin. 1971) 
indicate that the polar fields may not consist entirely of open field lines. 
Figure 3 shows the polar field topology based on the optical technique. 

On the other hand the Stylab coronograph shows Llie polar regions as 
large coronal holes which are evidently stable at least over the duration 
of the Stylab mission (Newkirk, 1975). This observation implies that 
most of the field lines over the poles are open. 

B. Solar Cosmic Ray Propagation at High Solar Latitudes 

A recent statistical study of solar flare cosmic ray events by 
Van Hollenbeke ct al . (1975) has revealed that: 

1. The onset time of 20-80 MoV proton fluxes at Earth systemati- 
cally decreases with heliolongitudc away from a minimum 
point at 50° ± 30° W of central meridian (I'igurc 4). 



2. The intensity of solar cosmic rays at Earth is strongly depen- 
dent on the distance in longitude from the flare site. The 
preferred connection region is 20° W to 80° W longitude with 
some of the spread being due to the variable solar wind velocity. 

3. Flai^ sources on the back side can sometimes be identified 
by vezy large X-ray and radio emission regions which extend 
above the limb. Particles from flares located away from the 
preferred longitude region as much as 150° are sometimes 
detected at Earth in this way. 

4. The energy spectra of the protons from individual flare events 
loci ted in the pi*eferred connection region (20 to 80° W) are 
remarkably similar: most of these spectra are fit fay power 
law exponents between -2 and -3, 1 (Figure 5). 

5. The energy spectra of flare parcicles measured near Earth 
changes systematically with distance from the preferred- 
connection longitude. At 40° E the range of power law 
exponents is -3.7 to -5.0 (Figure 5). 

The authors have eliminated interplanetary diffusion and adiabatic 
cooling as causes for the above observations. They explain the results 
in terms of energy dependent diffusion of the flare particles from the 
flare site through the coronal magnetic fields out to the spiral inter plane- 
tary field lines that connect to the observing site. The progressive softening 
of the flare particle energy spectra is attributed to the more rapid loss 


- 235 - 



of the higher energy particles onto spiral field lines lea^ring in the vicini^ 
of the flare site. The lower energy particles are preferential^ retained 
so that as the diffusing particle population ages, its average energy and 
spectral slope decreases. 

Such observations will make it possible to specify certain physical 
conditions in the coronal regions, especialfy the particle diffusion coef- 
ficient. Experiments at lower energy can determine the total amoimt of 
matter traversed fay the particles and thus a physical model of the diffusion 
in the high latitude corona can be built up. 

The out-of-the ecliptic missions offer the possibility of doing an 
analogous experiment in which latitude is the variable, thus complementing 
the longitude studies. In the latitude case the situation is even more 
complex but also should be more revealing about the spatial structure 
and characteristics of the solar coronal magnetic fields. To reach high 
latitudes, the particles will have had to propagate through the transition 
zone as well as appreciably into the polar fields. There are bound to be 
great surprises in such observations too far outside our present experience 
to anticipate. Nonetheless. Figure 6 attempts to qualitatively indicate 
what might be the outcome of solar cosmic ray observations at high latitudes. 

C. Origin of the Interplanetary Particle Populations 

Figure 7 shows the so-called quiet time proton energy spectrum in 
interplanetary space. At the lowest energy end there are the solar wind 
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protons. At high energies (greater than 10 GeV) the protons almost 
certain^ belong to a galactic population of particles, and at the very 
highest energies, the particles msyr even be intergalactic in distribution. 

In the range 30 MeV to about 10 GeV the protons are probabty galactic 
bat their intensity' is strongty modulated by the Sun's magnetized plasma 
wind, at least near the ecliptic plane. The range from 2 to 100 MeV has 
recently been studied extensivety by toe University of Chicago and 
Goddard Space Fli^t Center Groups. A detailed stucty by J. Kins^ (1 970) 
is based on the hypothesis that the protons in the energy interval come 
from two distinct populations. Figure 8 shows some spectra in this 
energy regi<m. Those observations may be represented by 

^ = F e"® + F E® (1) 

dE s g 

dJ/dE is the number of protons per square centimeter-steradian-second- 

MeV. F^, F^, s and g are parameters that give best fits to the obser- 

vatiems. It is fo^ind that F and s are highty variable while F and g 

s g 

vaiy much less. The interpretation is that the first term on the right 
hand side of Equation (1 ) represents a proton component of solar origin 
while the other term is presumed to be galactic in origin. 

No doubt the proton fluxes in this energy range will prove to be 
quite different when observed at high heliographic latitudes. There 
detectors presumably would be far removed from the solar sourc'C regions 
and the reduction of intensity in the galactic component by modulation 
should be much diminished. 
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Figure 9 shows the energy spectrum of electrons at times of 
lowest solar flare activity often referred to as the "interplanetary electron 
spectrum". As in the case of the protons* the highest energy electrons 
are presumed to fill the galas^, while the lowest energy electrons are 
known to be of solar origin — they are the solar wind electrons. Between 
these two components the situation is complicated, perhaps even more 
so than in the case of protons: 

1. The solar wind electrons have a non- Maxwellian tail at energies 

above about 100 eV. This tail appears to extend to 1 or 2 keV. 

-5. 1 

It can be described by a power law of E * . 

2. The non-Maxwellian becomes submerged in a new particle 

population which extends to about 1 MeV. This spectnun 

_3 

approximately follows the power law E 

There appears to be no spectral flattening of the 2 keV 1 MeV 
component as there is in the case of the protons. 

The stucfy^ of the complex electron spectrum would be greatly 
advanced by observations at high solar latitudes. Much will be learned 
by investigating how the electron spectrum changes relative to the proton 
spectrum. 

Among the questions we can now pose about the interplanetary 
fluxes are: 

1. Will the medium energy component decrease significantly as 
the detectors move away from the active sunspot zones? 

11 ,;] Y OF TH« 
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. From obserWng such a change ^stematicalty with latitude 
can we learn more about the sources and the transport of the 
particles through the coronal magnetic fields at high latitudes? 

3. Can we identify a component that is generated by instabilities 
in the solar wind? Possibly the 190 eV to 1 keV non- Maxwel- 
lian tail arises in this way. As the solar wind flow changes 
as a function of latitude does the character of the non-Maxwel- 
lian tail change? Possibly even the 2 to 10 keV region 
originates, at least in part, from internal energy sources 

of the solar wind ly means of plasma microinstabilities. 

4. Is there an electron component, for example the 10 to 100 
keV region, which is supplied, at least in part, ly strong 
cosmic X-ray sources? 

The high latitude electron and proton measurements will be free 
of planetary source effects: the Earth's bow shock and magnetosphere 
and Jupiter's magnetosphere. This is another important motivation for 
out-of-the ecliptic missions. 

D. Shock Waves and Other Interplanetary Phenomena 

A variety of complex energetic particle phenomena take place in 
interplanetary space. Although they must involve basic plasma processes 
they are not yet well understood. The motivation for further study is 
strong, however, since these processes must occur throughout our galaxy 
in ^sterns of similar and larger scale sizes. 
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One such phenomenon is the acceleration of particles by shock 
waves of solar origin in inter plane taxy space. Figure 10 illustrates 
this effect in a particular case. The important observational features 
in this event are (R. E. McGuire* Ph.D. thesis* University of California. 
Berkel^, 1976): 

1. The shock velocity was 520 km/p and the shock normal was 
close to the Sun- Earth line. 

2. Preceding the shock passage the electron and proton flux 
increased by a factor of 10 due evidently to a flare. 

3. Twenty minutes before the shock passage the protons above 
200 keV increased fay a factor of 20. The flux maximum 
occurred about one gyroradius in front of the shock. A mo»'e 
gradual increase in the medium and high energy electron 
flux also occurred. 

4. Electrons 0. 5 to 14 keV increased 20-fold just behind the 
shock front. This flux increase may be associated with dissi- 
pation of energy by the shock. 

Although several detailed models for shock acceleration exist, 
none seem completely satisfactory. By observing shock effects on par- 
ticles at high solar longitudes where the characteristic of the interplane- 
tary field are presumably much different it should be possible to arrive 
at a satisfactory shock acceleration model. This is a problem of Jirst 
order importance to cosmic ray studies since shock acceleration occurs 
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in many cosmic systems. For example* shock waves probably play an 
important role in the production of energetic flare particles: shocks may 
accelerate the entire fast electron component. And it is known that 
shock waves carzy energy into interplanetazy space in amounts equal to a 
large fraction of the total flare energy (Hundhausen and Gentiy* 1969). 
There is still no satisfactozy theoretical solution to the problem of mag- 
netic field line mergiz^ as an energy source for flaz^ phenomezia and it 
thus becomes vitalfy important in the study of the flare process to know 
as much as possible about shock acceleration. 

At high solar latitudes the magnetic field will make a smaller 
azigle with respect to the solar wind flow direction* on the average* com- 
pared to lower latitudes. Also* one expects that the power spectinim of 
fluctuations in the interplanetazy field to be considerably redu«'ed at high 
latitudes as compared to the latitude range of the activi^ zones. These 
differences can then be used to explain differences in the shock accelera- 
tion of particles as a function of heliocentric latitude. 

There seem to be several solar-interplanetazy phenomena which 
emerge clearly during part of the solar cycle but then become less apparent 
or disappear later in the solar cycle. In his presentation at the sympo- 
sium Dr. Hundhausen (1976) mentioned such an effect in the solar wind. 
Another example is the abundance of long lived particle streams (recur- 
rence events) early in the solar cycle (McDonald and Desai* 1971). Yet 
another example is the electron- proton "splitting" effect reported by 
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Lin and Anderson (1967). This effect occurs in solar co-rotating streams 
that appear near Earth one or two days following solar flares. The elec- 
trons in these streams are displaced to the west of the protons and thus 
are observed to pass over the detectors before the protons do. This time 
displacement is a few hours up to 10 hours (a few degrees in longitude). 

An example of electron- proton splitting is given in Figure 11. 

None of the proposed mechanisms has been shown to quantitatively 
account for this effect, and no explanation has been offered for their 
apparent disappearance later in the solar cycle. Lin and Anderson (1967) 
thought that the effect could be due to the larger drift velocities of the 
protons compared to the electrons (the larger gyroradius of the protons 
means that the protons sample any gradients in the magnetic field to a 
greater extent). Jokipii (1969) proposed that the splitting can occur in 
interplanetary space due to combined gradient and curvature drift in 
the interplanetary magnetic field. 

One expects that these phenomena and perhaps others still unex- 
plained can be effectively studied by an out-of-the ecliptic mission. 

E. Jupiter Swing-by Out of the Ecliptic Missions 

For a two-spacecraft launch in late 1980 by a Titan-Centaur 
vehicle as described by the NASA -Ames Research Center group. Figure 
12 shows when the spacecraft would reach high solar latitudes with 
respect to the sunspot cycle. A 1980 launch is probably not ideal for solar 
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cosmic ray studies since the spacecraft arrive at high latitudes near the 
expected decline of solar cycle 21. However, it is worth noting that 
several veiy large flares occurred in this portion of solar cycle 20 
including the veiy great August, 1972 flares. The spacecraft remains 
above latitude 35° N and 35° S for somewhat more than 2 years. 

Figure 13 shows where the spacecraft are positioned in solar 
latitude on a Maunder' s butterfly diagram. This figure shows that there 
is some advantage in arriving at these high latitudes late in the solar 
cycle since the activity zones have retreated toward the equator. 
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Figure Captions 


Figure 1. The intermediate solar corona. This photograph was taken 

following solar maximum and before solar minimum. At this time 
the polar plumes are clear^ \nsible. Their appearance suggests 
that the polar zone magnetic field lines do not close in the vicinity 
of the S\m. (Yerkes Observatory photograph* ) 

Figure 2. Synoptic charts of the polar magnetic fields H|| /cos 6. Solid 
line&i represent N polarity, dashed lines S polarity. The isogauss 
levels are 10 and 20 G. The areas covered by the observations 
are enclosed by curves with shadings on the outside, (a) North 
pole. Observations from August 1-23 1968. (b) South pole. 

Observations from July 31 -August 23, >8. (From J. O. Stenflo, 

Observations of the polar magnetic fields, in Solar Magnetic Fields 
(ed. R. Howard), D. Reidel Publ. lAU Symposium No. 43). 

Figure 3. Schematic map of the polar coronal magnetic fields (July 26, 
1970). The prominences and filaments observed on spectroheli- 
ograms are shown on the map. (From P. Charvin, Experimental 
study of the orientation of magnetic fields in the corona, in Solar 
Magnetic Fields (ed. R. Howard), D. Reidel, I.A.U. Symp. No. 43). 

Figure 4. The difference AT of the time uotween onset of 20-80 MeV 

m 

proton and maximum particle intensity is plotted as a function of 
the heliolongitude. The solid line is a least square fit through 
the data. It shows a minimum at ~ 50° ± 30° W of the central 
meridian. 
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Figure 5. Variation of the spectral index y in the 20-80 MeV range as 

P 

a function of the heliolongitude The open circles are 'long 
rise time events' with a rise time longer than 24 hours. For 
these events, effects of interplanetary propagation may be sig- 
nificant. The dashed countour lines enclose 92 % of all 
the other events. The solid line is a least square fit obtained for 
them. represented approximately by “ 

2.7 [1 +a\/2]. 

Figure 6. This figure idealizes how solar flare particle fluxes might 

depend on latitude for various coronal magnetic field configurations 
which vary with heliographic latitude. 

Figure 7. The interplanetary quiet time proton spectrum is made up of 
several components: at the lowest energy the solar wind protons 
and at the highest energies, the galactic cosmic rays. In between 
is a spectrum which probably has solar origin. 

Figure 8. This figure shows the region of presumed overlap between 
interplanetary protons of galactic and of solar origin. 

Figure 9. The interplanetary quiet time electron spectrum is also made 
up of several components; the solar wind electrons at the lowest 
energies and the galactic component at the highest. Several com- 
ponents appear to exist at intermediate energy. Little is known 
about their origin. 



Figure 10. An interplanctazy shock (speed 540 km/sec at 1 All) that 
produced large effects on particles. Protons in the 0. 3 to 0. 5 
MeV range are swept up or accelerated by the shock. A blanket 
of hot electrons appears just behind the shock possibly the result 
of energy dissipation fay the shock. 

Figure 11. This figure shows an example of electron -proton splitting 
early in solar cycle 20. The particles arrive at Earth in a co- 
rotating stream but the electrons are displaced toward the West 
by a few degrees. The origin of this effect is not known and it 
appears to become less frequent later in the solar cycle. 

Figure 12. This figure shows that for certain Jupiter swing-by out-of- 
the-ecliptic missions launched in late 1980, the spacecraft will 
arrive at high (> 35®) latitudes in late 1983 and remain there for 
somewhat over two years. 

Figure 13. For the same missions the spacecraft is seen to rapidly 
rise past the sunspot zones which, late in the solar cycle, have 
moved close to the equator. 
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Abstract t 

Observational results on the East-vest effect are summarized and 
discussed in the context of existing models of coronal propagation. 

The variation of the number of events with solar longitude is 
soirprisingly similar for p^o^ticles covering a large interval of 
rigidities* Over large longitudinal distances, time delays to the 
event onset and maximum intensity are independent of energy and 
velocity. This has important implications and vill require probably 
a transport process vhich is determined by fundamental properties of 
solar magnetic fields, e.g. reconnection processes betveen open and 
closed field configurations. 

A fit of Reid*s model for diffusive propagation in the corona to 
the observed delay times gives a (tvo-dimensional) diffusion 
coefficient K corresponding to r //. hours (r a distance 

of the thin diffusing shell from the Center of the sun). Limitations 
of the diffusion model are given by the existence of a fast propagation 
region vhich nay extend up to 40... 50” from the flare site, by the 
possible existence of an energy independent drift process, and by the 
influence of solar sector boundaries. The relative role of open and 
closed field configurations is extensively discussed. Some evidence 
is presented that the acceleration of protons to higher ( ^10 MeV) 
energies is related with a shock vave traveling in the solar atmosphere. 

The importance of measurements performed from spacecraft out of the 
ecliptic plane is stressed, in particular vith respect to the 
fundamental problems of particle acceleration in the flare process 
and for tinder standing fundamental dynamical properties of large-scale 
solar magnetic fields. 


« 

Extended version of a talk presented at the "Symposium on the 
Study of the sun and Interplanetary Medium in Th’^ee Dimensions", 
Goddard Space Flight Center, USA, May 15/16, I97t>. 
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1 . Introduction 

The story of coronal propagation begins with the East-West-effect 
for solar cosmic ray events: with increasing longitude of the 
parent flare on the Eastern hemisphere of the Sun the number 
of events detected at the earth decreases considerably, ai*d, 
for those events which are detected, the delay between the 
flare and the arrival of particles at the earth increases 
(see Burlaga, 1967, for asummary of some earlier results). 

The reason for the East-West-effect is obviously the asymmetric 
nature of the interplanetary magnetic field with respect to 
the central meridian on the Sun. For average solar wind 

conditions, a bundle of interplanetary magnetic field (iir.f) lines 
observed near the earth connects back to a point on the Sun 
which is close to 60® W. However, the controlliiig nature of 
the regions close to the Sun for the azimuthal propagation 
of energetic particles became only clear, v/hen it was established 
that the propagation of energetic particles in space occurs 
preferentially along the imf . The arguments for negligible 
particle mo'‘'ron perpendicular to the imf have meanwhile been 
summarized by various authors. In addition to the arguments 
presented e.g. by Roelof (1974) we wish to point out that the 
variation of delay times with solar longitude is independent 
of energy , which is another stror argument against interplanetary 
perpendicular diffusion (Relnhard and V/ibberenz, 1974, Ma Sung 
et al, , 1975). 

These effects of ‘'coronal propagation" whi<'V; depend on the 
relative azimuthal distance between a parent flare on the Sur\ 
and the point in space where the energetic solar particles 
are observed, could be studied so far only nr. a function of 
solar longitude. It is the purpose of this x^^iper, (i ) to 
summarize the observational results and to i-rder them v/ith 
respect to existing models, (il) to point out which fundamentally 
new results we should expect by studying va. iations with solar 
latitude , l.e. by using observations in space out of the 
ecliptic piano. 
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We shall start in section 2 to susunarize the methods • by 
vdiich various transport processes can be separated. In section 3 
we discuss the statistical methods, where the variation of 
characteristic parameters of solar events with solar longitude 
is studied for a large number of events. Different models 
have been developed to describe the average longitudinal 
variations. An independent method described in section 4 
ccnsists in detailed studies of individual events, including 
multi-spacecraft observations at different heliocentric longi- 
tudes, simultaneous intensity and anisotropy measurements, and 
the relation to observations of solar surface structures. Some 
indications to the acceleration process itself are treated 
briefly in section 5. Finally we summarize the open questions 
in section 6 and try to relate them to studies of solar particle 
events off the ecliptic plane. 

2. Separation of various transport processes 

In the sequence of events between the first acceleration of 
particles on the Sun and their final observation in space, 
we can ask different questions: How is the solar atmosphere 
filled up with energetic particles follov/ing the original 
acceleration process? How do the particles escape into space? 

How can the interplanetary propagation be separt ted from the 
solar transport processes? Let us start with some terminology 
related to the different steps. 

(a) The acceleration process is visualir.ed in many models as 
a pulse-li'.ce process limited in spatial extent to the flare 
area Itself, approximated by a delta-function in space and 
time. In principle, the acceleration could also occur over 
an extended area in the solar atmosphere for long periods of 
time (see below). 

(b) The accelerated particles spend a certain time in the 
vicinity of the Sun. We wish to make a distinction between 
pi opagatlon when they move away from the acceleration region 
and finally occupy a large area on the Sun, and storarre when 

the particles remain confined to a certain region. The difference 
Is depicted sc' ^matically in Figure 1. In reality, we may have 
a mixture of roth processes. 
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(c) The number of solar particles observed in space is determined 
by the probability per unit time that a particle v;jll leave the 
solar atmosphere by finally reaching an open field line leading 
out into space (Reid, 1964 ; Nevkirk, 1973 ). This release mechanism 
has a very important influence on the azimuthal distribution 

of particles, because it also determines how many particles are 
left for further propagation along the solar surface. 

(d) Sufficiently far away from the solar surface v/e only find 
open field lines which lead out into the interplanetary :.cdiuir>. 
Along this ’’source surface” (Newkirk, 1973) the processt ■ 

(a) to (c) define an i nfection '^ini'tion K(7 , t). I.easuro- 

ments on a single spacecraft connect back to a certain soiar 
longitude / close to the solar equator arid for a riven 

solar v;ind velocity the connection longitude 7*(t) as a function 
of time may be determined (see e.g. Nolte and Hoelo: , 



A singltt tatellitc therefor* "sees” an iHulection profile 
l(t) ■ M(i^(t), 0* t), Roelof and Krittigia (1973) have described 
an effective method to separate real longitudinal changes 
d R/d0 from changes in the injection function dN/^. Since 
dl/dt « (dK/d^)d0/dt -I- dN/^» the relative contribution of 
the two terms to the observed dl/dt depends on the motion 
d0/dt of “^e connection longitude. For a negative gradient 
in the solar wind velocity a solar wind "dwell”, d0/dt 
ist very souill, so that dl/dt er^H/^t. For a positive gradient 
in V^, d0/dt is large, the second term can be neglected, and 
one may directly construct the coronal distribution N(0). 

A more direct method for determining N(0) is of course the 
use of multi-spacecraft obsezvations (McCracken et al. , 1971; 

Bukata et al, 1972), in particular if combined with actual 
solar wild measurements for determination of the connection 
longitude 0(t). Results of this method are summarized by 
Roelof (197^). 


(e) It is clear that any attempt to determine the (coronal) 
injection profile, . (t) or the related longitudinal distribution 
N(0) has to start from observations in space, and, therefoi*e, 
one first has to separate the effects of interplanetary 
propagation . Methods available to perform this separation are 

(1) statistical studies in which the properties of solar 
particle events (maximum particle flux, times of onset 
and of maximum flux, shape of energy spectrum etc.) 
are ordered with respect to the longitude of the parent 
flare; 

(2) multi-spacecraft observations, by which longitudinal and 
temporal changes can be separated; 

(3) simultaneous intensity and anisotropy measurements, v/hich 
allov.' to separate long-lasting solar injection processes 
from long-lasting interplanetary storage j 

(A) the "mapping" of observed interplanetary particle fluxes 
to the high coronal source longitude, by vasing the 
simultaneously measured solar wind velocity. 



We shal l sta rt with a summary of results obtained by method (1), 
TTrTclJTiriilah c si' of the largest amount of observational data and 
gives insight into the average behaviour, i'ethods (2) to (U) con 
then be used to check predictions of models which have been 
developed and to provide additional insight into the relation 
with certain features observed on the solar surface. 



3. Statistical studies of longitudinal effects 

Let us first discuss how the total numher of observed events 
varies with solar longitude. Figure 2 shows the longitudinal 
distribution of solar particle events for four different sets 
of observations. The dashed line is for non-relativistlc 
electrons (after Lin» 1974), the dotted line for relativistic 
electrons (after Simnett, lS7^i plotted for longitudinal bins 
of 10® or 30^ • respect ivel> The full line is the original 
curve of Van Hollebeke et al. (1975) for 20-80 MeV protons, 
and the hatched area indicates results for ground level neutron 
monitor data (GLEs), as taken from Poraerantz and Duggal (1974). 

All four curves show the largest number of everts observed 
when the parent flare is on the western hemisphere of the 
Sun, with a broad maximum somev/here between 30 and 90®V/. The 
number of events clearly decreases as one goes to the Eastern 
hemissphere and beyond the West limb. Note that because of 
the difficulties of flare identification no electron data have 
been plotted beyond. 90®W. 

The overall similarity of the curves is rather surprising; 
only for the non-relativistic electrons the decrease in the 
number of events seems to f tart for a more v/esterly longitude. 
The similarity in the other three curves suggests a common 
propagation characteristics for the different particle species, 
which cover a range of Larmor radii of at least three orders 
of magnitude. Apart, from the clear decrease of the distributions 
of Figure 2 east of about V.30 solar Ic-ngitur’e, the functional 
shape of the curves cannot be determined precisely. From a 
different set of data. Smart et al. (1975) have fitted a 
Gaussian distribution xo the longitu.tinal distribution of 154 
flares on the visible hemisphere of the Gun, which have 
produced proton events, Th^’find th'^ highest frequency of 
flares grouped betv/een W30 and W40, and a standard deviation 
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of the Gaussian of 55^. A Gaussian curve with these parameters 
supplies a rather good fit also to the three similar curves 
in Figure 2. Prom this fit, one can now extrapolate to the 
invisible hemisphere of the Sun, with the result, that 
rou^ly 20 of all events should not be associated with a 
flare on the visible hemisphere of the Sun (for details see 
Smart et al., 1975)* A similar result has been found by 
Van Hollebeke et al. (1979)> 


As we shall see below the decrease in the number of events is 
mainly determined by the escape rate of particles from the 
Sun, so that distributions of the type sho\m in Pigur' 2 can 
be used to determine the escape rate. Ve note then thot tno 
similarity in the distribution function for different particle 
species suggests escape rate s wh ich do not depend strongly 
on particle type and energy. j us turn now to the variation 
of characteristic times of solar events with longitude. If 
the parent flare of a sf«lar particle event is located on the 
Eastern hemisphere of the Sun, the arrival of energetic 
particles becomes more and more delayed v/ith increasing solar 
longitude (Burlaga, 1967; Englade, 1971; L'atlowt, 1971; Barouch et 
Simnett, 1972; KcKibben, 1972; Lanzerotti, 197?; Reinhard 
and Wibberenz, 1974; Ma Sung et al. , 1975). 


As discussed above, azimuthal propagation ox particles in 
the interplanetary mediu:.i cannot account f‘>ir tt.e observations. 
The first quantitative model for- particle rope, Ration in 
a surface layer around, "^he -un vms devclopi-vi by i.eic (r»b't) 
and extendod by A;^ford ( I9cb). Tn iieic's i..' del, the in.iecTicn 
function can be written a.s 



Here is the loss t.ime which descrihrs ti e cscrye of 
particles into interplanetary .space, K , j •• thr corom l 
diffusion coefficient, and r^ iho <U.stance of tr.- c; i i : 
layer from the center of the Sun. So r^‘/K -r^. • •• -i. or; i c 
time it takes the particles to diffuse by .cn angvilar distance 


2fi7- 



of Y 2 * ra<J ■81^ from th« flare origin. In vhat follows, we 
put the flare origin at ^ « 0 and measure the angular 
distance ^ from the flare to the root of a bundle of field 
lines leading out into space where an observer is located 
at longitude 0^^ and latitude %«0. Note that flares occur off 
the sola.* equator, say at longitude 0c and latitude O,?* so 

p p 2 * o * 

that ^ s ®p ^h#* average Gj^jJr20 ; in what follows, 

we shall disregard the depeiidence on solar latitude and put 

0p - 0p. Observationally it is not possible to find 

systematic differences over a 20^ angular interval. 

Both and Xq sre assumed independent of but may vary 
with particle energy. Reid (1964) originally obtained an 
estimate of Tq = 3.f hours and Tr^J»l.2 hr. Applying Axford’s 
(1965) version of the coronal diffusion model, Kirsch and 
MUuoh (1974) obtain values for the coronal diffusion coefficient 
of the same order of magnitude for the event of Nov 2, 1969. 
Lanzerotti (1973) used the Reid model to describe the 
variation of onset times for 0.6 to 25 MeV proton.; with 
longitude. He esximated 7 hr$ T_ S16 hr. These values are 
probably underestimated because V,, refers to the bulk of particles 
to diffuse and net to the "first" particles which define the 
event onset . It is remarkable that shows almost no energy 
dependence. This energy independence of the coronal transport 
times over large longitudinal distances is meanwhile firmly 
established (McKibben, 1972; Reinhard and Wibberenz, 1974; 

Ma Sung et al., 1975). 

The most careful study so far in applying the Reid/Axford model 
to the East-West effect ha.s been performed by Mg and Gleeson 
(1975). They have replaced the plane approximation by diffusion 
in a real spherical shell, and the coronal injection profile 
is then used as the boundary condition for interplanetary 
propagation, taking into account anisotropic diffusion along 
the splraj. shaped Imf, convection, adiabatic deceleration, 
and corotation of the flux tubes past the observer. With their 
two stage propagation model they reproduce many features of 
solar events. They have used in particular the results by 
McKibben (1972) on the variation of the time-to-maxi..'.uin with 
longitude. 
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Their best estinatea of the coronal parameters are 
SOI* t^£lOO hr and lO$tT^£l5 hr. For the interplanetary 
propagation they obtain a value of the (radial) mean free 
path, %fhich corresponds to 0.03 AU for 10 MeV protons. This 
value is probably too small (see Wibberenz, 1974), since in 
many cases even for Western hemisphere events part of the 
delay is due to coronal, and not to interplanetary propagation 
(Reinhard and Wibberenz, 1974). But for the discussion of the 
Bast-West effect this difference is not critical. 

In Figure 3 we compare the computations of Ng and Gleeson (1975) 
with tjjj-values for >10 MeV protons (Reinhard and Wibberenz, 1974). 
For both curves (a) and (c) there is a well defined minimum 
close to V60, and a systematic increase on both sides of this 
idea] connection longitude. The apparent linear relation 
between t^^^ and ^ results from the escape term, which largely 
influences the coronal in;jection profile (1) for t^r^. The 
minimum value of t^^^ is related to the interplanetary propagation 
and could be shifted dovnav'ard by a factor of 2 or more, without 
changing the general conclusions. We can see now the limitations 
of the (modified) Reid model. For curve (a), a small value for 
the coronal diffusion time has been taken, 13 hours. This 

gives the desired flat longitudinal dependence for western 
events, but does not explain the large values of t^^j for eastern 
events. In model (c), the larger value for the coronal diffusion 
time, = 100 hr, gives the larger increase of t^^^ on the 
eastern hemisphere, but this increase now starts right away 
at W60 and is in strtng disagreement v;ith the bulk of data 
between 0 an<i W90. 

It had already been pointed out by Reinhard and Wibberenz (197^0 
that observational evidence speaks against a well defined 
minimum in the propagation times somewhere between W^O and W6n, 
and that very fast propagation with small or negligible coronal 
propagation times can be found for events where the parent 
flares are located between about 0 and W1«'0. The horizontal 
lines in Figure 3 are meant to indicate tiie existence of a 
fast propagation region** (FPR) where minimal propagation 
times can be found. The extent of this FHi may vary from one 
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•vent to the other. The existence of a very efficient solar 
propagation for certain- longitude ranges had also been pointed 
out by Pan et al. (1968). The "open cone of propagation" for 
>40 keV electrons found by Anderson and Lin (1966) and 
Lin (1970)* vdiich also has an extent of about 100°, may be 
identical with the FPR. A "region of preferred connection 
longitudes" ranging from about W20 to V80 is defined by the 
work of Van Hollebeke et al. (1973). It may Indeed be more 
appropriate to talk about a preferred connection to the 
acceleration region than about a fast propa/yatlon from it. 

We shall return to this point in section S. 

Let us return to the slow coronal propagation outside the FPR. 
One will obviously get a better fit to the data (see Figure 3) 
if e.g. curve (c) describing the coronal propagation does not 
start around W60, but with the same slope on both sides of 
the FPR. Clues to the possible nature of this slow coronal 
propagation may be found by a study of the dependence of t^^^ 
on particle parameters. Reinhard and V/ibberenz (1974) have 
studies time histories for 56 events in the energy range 
10-60 MeV and found that the most probable travel distance 
vtjjj is a linear relation of velocity, or with other words, 
tjij can be written as ® c-j -«• C2/V. The velocity dependent 
term in t^^ is found to be independent of solar longitude and 
describes interplanetary propagation. On the average, AU 

(with variations between about 2 and 10 /U). On the other hand, 
increases steadily with solar longitude East of central meridian 
and is independent of proton energy. 

These results have been confirmed by Ka ?ung et al. (19" 5) and 

extended to higher velocities (near-relativistic elev;trons) and 

to the inclusion of the event onset times They stuv/ that 

for the onset times a similar relation namely 

+ A_/v, with two additional con -.tants and E, , 

on on on on on 

which have to be determined for each event in additional to 

<=1 (s =2 (3 "max’- *ot. *niax ''“'‘y 

systematically with solar longitude, so tliat they can be taken 

to describe interplanetary propagation. The numerical values 

( Oon^ “ ^ ^ confirm that interplanetary 
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propagation plays a relatively minor role up to the time of 
maximum particle intensity at 1 AU. and 3^^ - c<| deaerlhe 
the onset time and the time of maximum of the solar injection 
profile, respectively. Both parameters are relative small on 
the Western hemisphere (of the oz*der of 1 hour if averaged 
over many events), but increase systematically with increasing 
longitude on the Eastern hemisphere. Both parameters are in- 
dependent of energy and/or velocity of the studied particles 
(0.5 to 1.1 MeV electrons, 4-80 MeV protons). 

The energy independence of the coronal transport puts severe 
limits on the possible physical mechanisms. Any particle motion 
in a given static magnetic field configuration gives transport 
times with (velocity) as one factor. This excludes e.g. 
gradient or curvature drift as one basic mechanism as v/ell as 
the current sheet diffusion (Fisk and Schatten, 197?), since 
a diffusion mean free path X which varies inversely with 
velocity (to give K'^vA * const) independent of the particle 
type is physically tinrealistic. Ma- Sung et al. (1975) propose 
a "leaky box model" with a diffusion coefficient k^‘''(aL) /^r , 
where is the scale size of the boxes \;hich open randomly 
on a time scale dr. The leaky boxes could be idealized models 
of large coronal magnetic field loops, and the process of 
field line reconnection provides the r'lndom opening of the 
boxes. In this model, the direction and speed of transport 
processes is not governed by the particle parameters (like 
velocity or rigidity) , b\>t by fundamental properties of solar 
magnetic fields; it is the randomness of the reconnecticn 
process which would be i-c-sponsible for a dJffusion-lih.e 
behaviour, v/hereos the particle i;,otj.on is 'leterministic. liote 
that this concept night .also lead to non-i' ffusive processes. 

If there is only a small numbor of boxes i i’ chanr.els along 
v/hich the particles propagate. A more deli aiuinist ic pi-ocess, 
namely an energy- indepeiulc-nt drift, v/hicli :-,hould act in 
addition to coronal diffusion, has been proposed by Reinhard 
and Wibberenz (1974). One of the original :-.upports for this 
idea, namely the linear relationship betv/ecn t^,, and , can 
no longer bo maintained, since because of the escape processes 
a quasl-linear relation can also be simult.ted in a diffusion 
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model (Ng and Gleeson, 1975). However, the total injection 
time profile depends critically on the nature of the processes. 
Reinhard and Roelof (1975) have studied the relation between 
onset times and maximum times and confirmed the necessity to 
include solar drift processes. ‘Sheir (linear) drift-diffusion 
model contains a drift rate a diffusion time « r^ /K^, 
and a loss time VThen the corresponding injection profile 
is convoluted with interplanetary diffusion processes, one 
gets a good description of the measured time profiles of 
10-60 MeV protons. The parameters of the solar injection 
profile and the interplanetary scattering mean free path 
are independent of proton energy. From a fit to several Eastern 
events, Reinhard and Roelof (1975) determine average values 
for the coronal parameters as = 0.42 hr /grad, IJ, » 400 hi', 
“13 hr. 

This choice of takes care of the observed widening of the 
time profiles with increasing solar longitude, whereas it is 
essentially the drift which determines the increase of the 
absolute time delay with solar longitude. The corresponding 
dependence tjj^(0) for the drift-diffusion model is indicated 
by the dashed line li* Figure 3. Here an extension of the 
FPR of - 40® has been assumed, centered around W50 longitude, 
so the increase sets in East of W10 and V/est of W90®. The 
average contribution of the inte-'planetary propagation 
corresponds to a time dtlsy of 4.4 hr for 10 MeV protons. 

The loss time of 13 hours is within the range assumed to be 
realistic by Ng and Gleeson (1975). Let us see how the loss 
time can be determined observationally . 

There are in gjrnal tv.'o processes by which for a given point 
on the solar surface the intensity is dir.ii ‘lished. (1) the 
lateral spread of pcu tlcles ctiuser. a corresponding decirease 
in surface density, (2) the injection into the interplanetary 
medium causes a general loss proportional to the number of 
particles present. These two point.s also determine the maximum 
intensity in the Injection function. If we take the 

Reld/Axford model seriously, sec equ. (1), *,/e can der-ive the 
follov/lng predictions; 
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(a) For small distances 0 the propagation times are small 

to that the escape term exp(-t/rj^) has no influence on the 
time of the intensity maximum c_ of the coronal injection 
function. In this case Cjjj ■ (r and we get from equ.(1) 

•WV * "tr- ^ 

(independent of the coronal parameters). Thus, a strong dependence 
of in the coronal injection profile on angular distance 

from the flare is predicted by the Reid/Axford model. Its 
verification by observations in space is limited, (a) by the 
convolution with interplanetary propagation, (b) by the fact 
that most flares occur remote from the solar equator (see above, 

©P«f20®). 

In addition, the existence of the FPR will initially fill up an 

extended area close to the flare site, vdiich will preclude the 

sharp dependence of on ^ as suggested by (2). In any case, 

the relatively flat distribution of the number of events with 

longitude on the western hemisphere (see Figure 2) favors a 

model with a moderate variation of N over small distances 

— — i ■ max 

from the flare. 

(b) For large distances ^ , it is mainly the exponential terra 
e7cp{-t/Zj) in equ. (1 ) which deter^.anes the decrease of 

(It should be noted that in the raore realistic version of the 
Reid model treated by fJg and Gleeson (1975), at late times 
this is the only term in the temporal variation, because the 
whole solar surface is covered with particles, there is no 
1/t-f actor left). 

Let us describe the average deci*ease of the number cf particles 
at the naximum of the Injection function wi'h P(9'). P(^) is 
allowed to vary with energy, but shall be the same for evei’y 
event. Here we measure 0 along the solar equator and let P(0) 
be normalized to 1 for events close to It is this 

variation 


N (0) 
max"^' 




( 3 ) 


hii.ITY OF THB 
PAGE IS POOR 
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tdileh determinta how the number of events detected varies with 
solar longitude (Reinhard and Roelofi 1975). For normally an 
"event" is identified vdien the number of particles during the 
maximum phase exceeds a certain threshold » which is determined 
by detector background and counting rate statistics. Let t)ie 
sise distribution (the number W of events where the maximum 
intensity of particles exceeds a given value be a 

separable function of and longitude, i.e. the longitudinal 
distribution is independent of the size of the event. 

“('W' 

With ■ T « threshold for detection of The event type under 

study we get W(T.O) > total number of events above threshold 
close to the preferred connection longitude around ,V/40. . .1/60. 
and W(T,0) « f(T)g(0) is the longitudinal f?.are distribution 
(see Figure 2). 

On the other hand, we have from equ. (3) V/(i',{5) = w|t/P( 0). 
Equating the two expressions for V.', we finally obtain 

f(p^^) ® f(T)g(0) (5) 

Relation (4) has been verified by Reinhard c.nd Roelof (1975) 
for protons >10 MeV, > 30 MeV, and>60 MeV, and they showed 
that f(N__„) can be described by a power lew (N,„„„)*^ with 
0.36. Inserting this into (5) they get r(0) = 
as relation between the longitudinal distriVution g and the 
size variation P. ffote that this relation i: plies that g(0) 
is rather insensitive to the .size variation l'(0). This is 
even more so, if we take an indc-’pemlent del* nninotion of the 
size spectrtim from Van nollehoko ot al. (1<^ '‘)); tliey obtain 
for the difierential size spectruin dv;/di:, with 

tC = 1.15 - .05, v/liich i!:ipiios that tlu; intf^'i'a! spectrum W 
is rather flat and^us»f .15 instead of 0.j;L. i.s uViove. 

The longitudinal distribution has been detei.'incd by Reinhard 
and Roelof (1975) as g(0) ~ exp(in/) with m -o.ni (degree)”"^ 
Insertion into the drift-diffusion model witli the average 
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paraiMtttrs oitod above leads to a loss time ■ 13 hours* 
independent of energy. In contrast* Van Hollebeke et al. (1973) 
have oonoluded from a variation of the observed spectral shapes 
of proton spectra with longitude* that the escape rate should 
be energy dependent. They find an average loss time of 

« 1.83 hours for protons with a mean energy of kO MeV* and 
an increase in the loss rate of 33-A3^ from 20 HeV to 60 MeV. 

The smaller loss time for the higher energies would then be 
responsible for the observed steepening of the spectrum v/ith 
longitude. 

Reinhard and Roelof (1973) do not find a systematic variation 
of the spectral slope with longitude. The reason for the 
discrepancy is not clear, it should be partly related to the 
use of a different set of data (difference in the threshold 
T for event detection; different selection criteria for "solar 
events"). In any case* it appears that a loss time of about 2 
hours is too small to be compatible with observations; the 
corresponding decrease in the injection function by a factor 
of 10 every 4.5 hours would make events from the Eastern hemi- 
sphere of the Svm praktically undetectable. Because of the 
intensitivity of the size distribution on (see the discussion 
following equ. (5) above) more direct deterninations v/ill be 
necessary (see section 4 for some indications). 

Let us close this section on the variations of solar e^r&nt 
parameters with solar longitxide. Ue have di;^cussed the delay 
times (onset and time-to-r.axinium) , -he general shaije of the 
intensi ty-tinic profiles, and the distrihut j.( n of the number 
of events with longitude, i/e have tried to i elate tlie overage 
behaviour of a large number of events to specific corona] 
propagation models. Lot us '«• nmiorize tiio essential aspects 
of the various models, in particular v/i tli rtcpecl to predictions 
of latitudinal solar variations. 

In the Reid/Axfoi’d-model there is just one j undamental process 
acting thougliout tho solar surface, characterized by a diffusion 
coefficient K.^, We had pointed out the difficulty to describe 
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simultaneously, with a unique value for the increase of 

t^ on the Eastern disk and the small values of t^ over the 

Western disk (see Figure 3)* In principle, this difficulty 

might be overcome by assuming a large variability of and 

by ascribing the small tg^-values to occasional large values 

of Tests of the Rexd model from measurements in the 
c 

ecliptic plane are partly restricted because one never scans 
A spacecraft measuring solar equatorial latitudes in 
the range 20-30® should see the strong dependence with distances 
from the solar flare predicted by the variation of (see 

equ. (2)) or a well defined second maximum related with the 
corotation for events East of the connection longitude (see Ng 
and Gleeson, 1975). 

On the other hand, small latitudinal variations would be 
expected if the FPR (fast propagation region) exists. Reinhard 
and Wlbberenz (1974) have suggested that the extent of the FPR 
is related to solar sector boundaries which have already been 
found to play an important role in determining the efficiency 
of coronal propagation (Roelof and Krimigis, 1973). The azimuthal 
distribution of the number of particles shov.'s discontinuities 
at certain solar sector boundaries, so that 1s has become 
possible to assign "access probabilities" to individual unipolar 
solar cells (Gold and Roelof, 1974; Roelof, 1974). Studies on 
spacecraft ”>ff the ecliptic plane and the technique of mapping 
intensities back to the iiigh solar corona (soe Roelof, 1974) 
should allow to observe the very existence i f " FPl;, its 
latitudiiijl e'd;f.-nt and the influence of the polev;srdc vsides of 
the unipolar solar cells. 


(?) Propagation_oyer_li!r£e__dif>t.;;neysi_diiryii ;ji onal;^^ 
transport? 


A very schematic distinction bnt;/een verioc. tJiC’dol.s can be 
seen in Figure 4, where the broad time dc'V*.- 1 opinent over the 
solar surface is plotted for the pure diffuc.ion model (left) 
and the drift-diffusion model (right). For implicity , ve 
neglect the influence of the loss term and merely indicate 
how a cliarncterlstic nn{:ular extent and llio average jiarticle 
density vdthin are expected to vary with time. 



The essential point in the "drift* process is that the transport 
process in the corona is not totally statistical (like in a 
diffusion process. %diere the net streaming of particles is 
simply proportional to the density gradient), but that there 
is a preferential bulk motion of particles into one direction 
superimposed. The filling up of one large "box” vfhich was 
initially empty might be one such process (see the ab( /e 
discussion about the "leaky box concept"), A pure drift, where 
all particles move into the same direction, is depict*id schema- 
tically in the outer right part of Figure 4. /n gxji-drift of 
the required order of magnitude (corresponding to velocities 
of 7 km/sec if ir * 1 solar radius) is not very likely. 

However, a very direct proof would be the systematic depletion 
of particles from the region close to the original flare (see 
Fig. 4, right part). The shift of time intensity profiles 
between spaceprobes separated in heliocentric longitude hns 
been shown to be consistent with normal corotation for four 
individual solar events in 1968 (McKibben, 1973) and does not 
require an additional drift. Moreover, the azimuthal gradients 
are in general positive when one approaches the heliocentric 
longitude of the flare (see McCracken and Rao, 1970; McCracken 
et al. , 1971 ). 

So one should still regard the drift-diffusion model as 
hypothetical and a convenient mathematical description to 
describe coronal injection profiles at one Jongltude. The 
consequence of a depletion around 0^0 at Icte times l.as still 
to be confirmed. 

Let us turn nov; to the question of dimensii-.iality of transport. 
The distinction is sliown in the tv/o left sc;. ernes in Figure ?. 

The question Is still opc-n v/hethoi- the con i:al propagat J on is 
relateci to a fundamental solar process vdii,' . acts sit.iiarly 
all over fne solar surface, of if it is ro I :tcd to spt-'cific 
processes which are typic:_l for the activity belts say. If the 
propagation is somehow related to large scale magnetic field 
loops, one should expect a prefeienti al proi .agation along the 
East-Vi'est direction because of the preferominl orientation 
of the loops in thl.e way. This would favor /i one-dlmer’.sionnl 
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proptigati<m in a limited latitudinal range (see second sketch 
from the left in Fugure 4). In this case* observations beyoud 
about 40 or 50 ^ in latitude would hardly show an detectable 
amount of solar energetic particles* and the large coronal 
holes found sometimes at the solar poles would be totally 
free of solar flare particles. 

On the other hand* coronal propagation might be related to 
a process which is occurring all over the Sun* similar to 
the supergraxiulation* or to the numerous current sheets and 
minute dipoles with average strength of 500 - 1000 Gau(^ and 
12 hour lifetimes (see Newkirk* 1975* for discussion). In such 
a universal process* v/e should more expect a distribution 
sketched in the outer left part of Figxire 4. It is clear that 
studies of particle populations at large heliocentric latitudes 
offer a unique opportunity to distinguish between the two 
fundamentally different possibilities. 

4. Detailed studies of individual events 

It is not intended to give a detailed account of the numerous 
studies of longitudinal effects for single events or for 
selected periods of time. V/e simply want to describe various 
methods and to describe a few results which give the necessary 
and important supplements to the statistical studies discussed 
so far. 
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NcCracken et al. (1971) have studied data fron four Pioneer 
spacecraft separated by <^180^ in heliocentric longitude. At 
late tines (> 4 days) in the eTents they still found strong 
gradients in longitude, with e-folding angles for 10 HeT 
protons of the order ^ » 30^. and no tenporal change in the 
relative gradient. Phis corresponds to a factor of 10 decrease 
every 70^. Note that this value of has only been directly 
detemined for two events. Van Hollebeke et al. (1973) from 
the study of a ouch larger number of events conclude that on 
the average the event sise for protons around 40 HeV decreases 
by about two orders of magnitude every 60^ away from the 
preferred connection region. This value is also consistent 
with an exponential gradient in longitude corresponding to 
a change of two orders of magnitude over a longitudinal 
distance of 40®... 60® (Roelof et al., 1975) for MeV protons 
and alpha particles. A comparison of these gradients with 
specific coronal propagation models has not yet been per- 
formed . 

Persistent anisotropies along the imf from the general solar 
direction have been among the first indicators that the solar 
source has to be described by a long-lasting injection instead 
of a delta-function in time (Bartley et al., 1966; Pan et al., 
1968; Krimigis et al. , 1971). The interpretation of these 
persistent large anisotropies depends on whether or not inter- 
planetary propagation can be neglected. Schulze et al. (1974) 
gave an example where the simultaneous fit of intensity and 
anisotropy data for 22-60 KeV protons during the Nov 18, 1968, 
event allows to determine the approximate duration of the solar 
injection as well as the interplanetary mean free path. As 
pointed out later (Schulze et al. , 1975) the interplanetary 
data arc relatively insensitive to the form of the solar 
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injection profile. A change in the form (not the characteristic 
duration) of the solar injection profile can be canceled by 
a suitable change in the interplanetary scattering mean free 
path without changing essentially the intensities or anisotropies. 

The situation is different iidien interplanetary scattering 
within the inner solar system can be neglected. For the 
scatter- free proton event of April 20, 1971, Palmer et al.(1975) 
could directly determine the solar injection profile and 
obtained a solar decay time of about 7 hours for 7.6-55 MeV 
protons. Roelof and Krimigis (1975) have pointed out that 
for low energy protons (^1 KeV) scattering in the inner solar 
system is almost absent. Here the magnitude and the direction 
of the anisotropies are used to infer the small interplanetary 
scattering, and by use of the ”inapping" technique conclusions 
can be drawn on the coronal injection profiles. Various time 
periods have been studies in a series of papers (Roelof, 1973; 
1974; Gold et al. , 1974; Krieger et al. , 1975). Typical coronal 
profiles are ramp-like structures, which are relatively smooth 
as long as the observer is connected to the same unipolar cell 
on the Sun, and sharp changes in intensity when a neutral 
field line on the Sun is crossed. These results are seen with 
corresponding time delays at spacecraft widely different in 
heliocentric longitude, (Roelof and Krimigis, 1973) v/hich 
confirms the spatial rather then the temporal structure of 
the profiles. 

The ordering of solar energetic particle data by solar 
stnictures observable in -f iltergrarns becomes also clear 
in the large solar events of August 1972, where Roelof et al . 
(1974) have studied flux histories for protons >1:. 5 MeV 
from Pioneers 9 and 10 and IMP 5. The different access of 
particles to regions on both sides of a solar sector boundary 
is clearly established. The crossings of solar sector 
boundaries are therefore in many cases seen by abrupt changes 
in the intensities of solar energetic particles, but they 
can also lead to marked changes in the riese or decay times 
of the total profile (see Reinhard, 1975b). It is very 
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remarkable that the solar sector boundaries, which are 
inferred from the -f iltergrams (McIntosh, 1972) and which 

are found to play such an important role for the access 
probabilities of solar energetic particles, do not in each 
case coincide with interplanetary sector boundaries. A 
possible explanation by chromospheric neutral lines which are 
not continued into an interplanetary magnetic field sector 
boundary, has been given by Roelof (197A). 

The irregularities which are related with the solar cells of 
different polarity certainly have to be superimposed on the 
overall dependence on solar longitude which we discussed in 
section 3. Formally this could be described by a variation 
rf the coronal parameters with solar longitude. This ‘ heme 
might offer an explanation for the observations of "anomalous” 
injection profiles of solar particles. Keath et al. (1971) 
have shown that the favored path for cosmic ray propagation 
in the March 12, 1969 event was about 40° east of the nominal 
Archimedes spiral line of force from the flare location. 

Palmer and Smerd (1972) also found a deviation from the 
"classical" picture, where the best connection into space 
should be close to the flare site. They explain the appear- 
ance of a prompt low energy proton component far away from 
the original flare by the triggering action of a shock wave 
travelling in the solar atmosphere. Cherki et al. (1974) 
find by analyzing the March 29, 1970 event that particles of 
different rigidity are elected at different longitudes on 
the Sun. Barouch et al . (1971) studied the onset times of 
6-25 MeV protons for several flares from the same active 
re;rion, and conclude that the mae^netic fields close to the 
active region should be considerably distorted from the 
nominal Archimedean field. 

These examples show that for individual events the release 
mechanism from the Sun may become very complicated, and that 
the models and coronal parameters discussed in section 
only describe the average characteristics over many events. 

’TY of thb 
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In all caaftSy however, the unusual particle escape from the 
Sun is thought to be related to processes occuring in the 
vicinity of the Sun» in particular in the coronal magnetic 
fields. As discussed already in section 3, large coronal 
TOf loops probably play a prominent role. Over which longit- 
udinal range and how far up into the corona these loops extend 
is not yet clear. Simnett (1974b) in case of the August 11, 

1970 event has suggested the existence of a stable loop 
extending several solar radii above the solar surface and 
about 100° in longitude. Two release points should exist 
for solar protons on both aides of this loop, with quasi - 
stable trapping inside the loon. 

Observations from Srfvlab have cast some doubt on the existence 
of stable loops of this extent. Ciiase «t al . (1P?S) hove 
studied one hundred loons detectable in so^t X-rays and 
show that the number of Interconnections decreases steeply 
for longer distances; the longest interconnecting loop 
extends over an angular distance of 37°. 

The question how far the loonn extend and which portion of 
the solar surface is covered v/ith "closed" or "open" config- 
urations is crucial for the whole propagation problem. The 
energetic particles perform gvrations about the field lines, 
and the transport of particles from one field line to the 
neighbouring one does not denend on whether the field lines 
are closed (l.e. return to the solar atmosphere) or open 
(i.e. lead out into interplanetary/ space). However, the 
number of closed field lines determines the amount of trapping, 
and once particles have been transmitted to open field lines 
they will eacane into apace. This means tliat an efficient 
storage mechanism, and a transno^’t which finally allows to 
fill practically the whole solar atnosrhere (see HcCracIcen 
et al . , 1Q71; MoKibben, 1^7?) ourht to be nnlv possible if a 
large fraction of the solar surface is "'dosed". This is 
confirmed by observational evidence; there should be a relative 
amount of open field lines of the order lO-l^ ^ in ecuntcrinl 
regions, ?5-40 averaged over the whole Sun (Nowkirlc, private 
communication). These figures are based on the potential 
(current-free) coronal field calculated from the observed- 
llne-of-sight fields at the pliotospheric level (Altschuler 
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and Newkirk, and tnnv have to be modified by th< 

influence of the exnandln^^ aolar wind. 

Newkirk (1^7“^) concludes that the ambient field configuration 
around active regions also determines whether or not protons 
escape from a given flare. It is found that among all flares 
proton flares have significantly more open field lines 
emerging from the vicinity of the active region. Newkirk 
assumes an ”ln.jection surface" of x "^6° centered on the 
flare and finds that the spread in longitude of open field 
lines is characterized by a full-width-at-tiO percent of 10® 
to 20® which is insufficient to explain the observed longit- 
udinal distribution of energetic particle events. One possible 
explanation for the discrepancy is that perhap.e the original 
in.iection surface must be larger. This would occur if shocks 
are the principal sources of energetic protons in the corona. 
We shall discuss this point in somewhat more detail in the 
next section. 

5. Relation to the accel*’ vrition process in solar flares. 

Let us first summarize some of the properties of coronal 
propagation along the ideas of Siranett (l974) or McKibben 
(197'^). There is a "prompt component” or "phase 1" of solar 
particle events. This is due to particles which either have 
been directly accelerated on open field lines or have been 
iniected onto open field lines shortly after the flare. The 
longitudinal extent where these prompt particles ore found 
should correspond to the "fast propagation region" discussed 
above. These initially injected pi^rticles give rise to a 
relatively short decay time (McKibiien , IQ?'")). 

There is a "delayed component" or "phase ?" of solar particle 
events. It is very pT’obable tliat tliose delayed particl€?s have 
been accelerated on closes rioid lines. Tlii;y tlien propagate 
in the corona, from one clo.aed conf iguru i j on to tiie next, 
maybe according to the concept of the leaky boxes discussed 
above, and from thereon tiiere is only a gradual release of 
these particles into space. If tliis release wei’e instead vo'y 
fast and efficient, we would never observe that flare 
particles have finally occupied essentially the whole inner 
solar systerol 
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There are two different solar decay timea for the two types 
of particle populations. Wibberens and Reinhard (197S) and 
Reinhard (1975s) show that by oonvoluting the solar decay 
prooeasea with interplanetary propagation one can quite 
naturally explain the exponential or quasi>exponential nature 
of the interplanetary decay, and that no "free escape boundary" 
around 2... 3 AU is required, which would be difficult to be 
reconciled anyhow with the Pioneer 10 and 11 observations. 

We turn now to the phase-1 particles within the FPR. It has 
been shown that the injection profiles for these particles 
are not of the delta-f uncti on type in tine, but finite in 
width (Palmer et al., 1975; Reinhard, 1975a). It is certainly 
difficult to distinguish whether these finite injection 
profiles stem from continuous release or continuous accel- 
eration of particles. But if the decay in phase 1 is much 
steeper than in phase 2 (McKibben, 1973)* then the particles 
cannot be replenished by the neighbouring storage region 
through the same process as the phase-2 particles. There 
might be a small storage region close to the flare with a 
different release mechanism for the phase-1 particles (this 
mechanism might exist because of the high degree of disturbance 
in the solar atmosphere following the flare). Or we have to 
assume that the injection profile directly gives the number 
of particles as they are accelerated. 

The second possibility is very Interesting with respect to 
the two-stage acceleration process, which is discussed in 
detail e.g. by Lin (1974). In the first phase non-relottvistic 
electrons are accelerated. If a sufficiently large number of 
electrons is damped into the chromosphere and lov;er corona, 
explosive heating occurs and produces an ejection of material 
and a shock wave which accelerates electrons and protons to 
relavistic energies. This picture is confirmed bv Svpstka anil 
Fri tzova-Svestkova (1974). They present convincing evidence 
that proton acceleration to higher energies (J^IO MeV) la 
closely connected with type II bursts, i.e. shock waves 
travelling In the solar atmosphere. Our above interpretation 
that the finite injection profiles might resemble the finite 
duration of the acceleration process itself would favor the 
shock acceleration model for protons, and It also explains. 
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why protons may be found on a wide lon/»ltudinal range of 
open field lines. 

She small longitudinal variation of the slopes of the 
proton energy spectrum discussed In section ? Implies that 
the particles of phase 1 and phase 2 have the same energy 
spectrum, which means that the acceleration process should 
work in the same manner on open as well as on closed field 
lines. In addition, any model of the acceleration proces 
should explain why the spectral slones for high energy 
protons and relativistic electrons are roughly the same, 

(see Van.Hollebeke et al., 197S). It would also he 
Interesting to see if the energy dependence of the decay 
times. as suggested by Reinhard (1975a) can be reproduced 
by a shock model acceleration. 

We close with a remark which emohazises the role which the 
fast propagation region raav play for the study of solar 
acceleration processes. The He -rich events which are a novel 
feature of solar cosmic rays (Serlemitsos and Balasubrahmanyan, 
1975) appear to be observed only for sufficiently small 
events and are only found when the parent flare is on the 
Western hemisphere (McDonald, private communication). This 
might indicate a snecific configuration near the Sun, where 
the acceleration process supplies He nuclei only directly 
to open field lines. 

6. flonclnsions. 

There are two important features of the interplanetary 
nropagation which allow us to study coronal transport nhenomenn 
(1) the motion of energetic particles rernendi cul nr to the 
Iraf is small, (2) the scattering; mean free path along the 
imf is large, so that details of the solar in, lection profile.s 
can be regained from measurements in space. 

We have tried to order the material by models for accel- 
eration, Injection, and propagation nrocepses. None of these 
models has been proven to give the real physical picture, 
because the underlying processes could not yet he identified. 

In this brief summary we shall concentrate on those open points 
which could be further clarified by spacecraft raenaurements 
performed out of the ecllntic plane. 
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(1) What 18 the nature of the "fnat nrnT^O‘'at1on r#*"ior'' 

(probably identical with the "onen cone of rrooopnt.ion'* or 
the ’’repion of preferred connection lonpltudea")? Is there 
also a region of nre^er^ed connection latitudes ? If the 
latitudinal extent la limited, what la the nature of the 
boundary? 

(2) If there is a transport process Involved within the fast 
propagation region, what is the nature of this process'!* 

Could it be a diffusive process with a surficieptly la»’'^e 
diffusion coeffi cl ent? In this case one should be able to 
see the relat.ivelv at'oonn' denendencf* or the mfixlmim intensity 
on angular distance 0 for small 0 as discussed in section ?, 
because the connection point of an ex-eclintic snaceornft 
comes closer to the active regions. 

(■^) Is the acceleration process directly responsible for the 
longitudinal width of the FPH and for the fast access of 
particles to open field lines? If particles are accelerated 
on open field lines by a travelling shock and may then 
escape into space, the extent of this "piompt” region should 
be determined by the distance which the shock can travel in 
longitude and latitude. V/ill we see narticles arrive very 
fast over the poles? 

(4) Which role do the solar sector boundaries plav^ Is th«»re a 
similar change In the access probabilities if one leaves a 
unipolar cell on the nortViem or the southern boundary? 

If there is one large uninolar cell (e.g. a coronal hole) 
extending from the pole to the equator, 1 m there the a;une 
access probability all ove’^ this cell, or is there a gradual 
or drastic change with latitude? Would one detect He -rich 
events over the poles? 

(5) What is the nature of the ener/'y-inJenerident slow coronal 
propagation over large distances :n 1 oror; tuiie? Is tlicre a 
drift nrooess (possibly related to electric ficl.i!;) involved? 

Are the time delays and the intensity dec rfane;; merely u function 
of the absolute angular distance between flare and observer, 

or are the variations tyricallv different into the Knst-West 
and into the North-.'louth direction” If th'^ lorre scale mf 
loops in the corona nlay an important role for the prona<'ati on . 
one would expect ouch systematic di f ferences , and then the 
preferred direction of the propagation and the latitudinal 
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extent up to which particles are transported may vary 
with the solar cycle y hecause of the difference in the 
average orientation and location of the loops. If the 
interpretation in section 3 on the importance of magnetic 
reconnection processes is correct « the study of large scale 
coronal particle transport should give insight into a 
fundamental solar problem. 
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Figure 1: Scheinatlc representation of "propagation" and "storage" 
processes. The main difference is in the lateral 
distribu':.lon of particles, whereas an observer close 
to the original acceleration region may see the same 
injection function in both cases. 

Figure 2: Number of solar events as a function of solar longitude, 
for different particle types and en--rgies. Data have 
been collected from Lin (1974), Van Hollebeke et al. 
( 1975 ), Pomerantz and Duggal (1974), and Sinmett (1974). 

Figirre ?: The time of the intensity maximum (t^) as a function 
of solar longitude. The collection of experimental 
points taken from Reinhard and V.'ibberenz (1974) is 
compared with calculations. Curves (a) and (c) are 
based on the coronal diffusion model (Reid, 1964) in 
the extended version of Mg and Gleeson (1975), curve 
(DD) is based on the combination of a fast propagation 
region vrith the drift-diffusion model (Reinhard and 
Roelof, 1975). 

Figure 4: Angular spread over the solar surface for different 
coronal propagation models. Temporal development 
from top to bottom. The influence of the escape 
process is neglected. N is the average particle 
density within the cross-hatched area, ^ is a 
characteristic ruaxiriuin distance froi:. the ilaie site 
reached after time t. 
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ABSTRACT 

Measurements of interplanetary dust via zodiacal light 
observations and direct detection are discussed for an 
out-of-eciiptic space probe. Particle fluxes and 
zodiacal light brightnesses are predicted for three 
models of the dust distribution. These models predict 
that most 'f the information will be obtained at space 
probe distSi ^es less than 1 A. Li. from the ecliptic p?ane. 
Joint interpi '-.ation of the direct particle measurements 
and the zodiacal light data will yield the best knowledge 
of the three-dimensional particle dynamics, spatial dis- 
tribution, and physical characteristics of the interplane- 
tary dust. Such measurements are important for our 
understanding of the origin and role of the dust in rela- 
tion to meteoroids, asteroids, and comets, as well as the 
interaction of the dust with solar forces. 
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Introduction 


The microscopic dust particles are an important constituent of 
our solar system. A knowledge of the physical and dynamical 
properties of these dust particles in three dimensions should 
aid our understanding of the origin and evolution of the plane- 
tary system and furthermore of circumstellar dust clouds as de- 
tected by infrared astronomy. In this, paper we discuss the type 
of dust measurements suitable for an out-of-ecliptic mission and 
the information which could be expected from these measurements. 

The interplanetary dust can be explored from space probes by two 
complementary methods : 

1. Direct detection of individual particles intercepted by a 
sensor. Velocity and mass parameters can be derived and, 
depending on the specific experiment, information on the 
chemical composition. Recent experiments have extended the 
limiting sensitivity down to masses of about 10~^^ g at 
impact velocities of approximately 20 km/s (Hoffman et al. 
1975) . 

2. Zodiacal light observations. Measurements of the brightness, 
color, and polarization of sunlight scattered by the dust 
particles give information on the average scattering proper- 
ties and the spatial distribution of the dust along the line 
of sight (ref. Leinert 1975). 

Direct detection has the advantage of defining a complete set 
of parameters for individual particles but the drawback of 
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sampling only a small number of particles. The zodiacal light 
observations sample a large volume of space, but the information 
on the physical properties of the particles is indirect. 

Combined data from these two methods can be valuable especially 
if the direct detection instrument is sensitive enough to cover 
the entire mass range which may contribute to the zodiacal light. 
The combined data can in principle specify particle velocities, 
orbits, size, mass rangv . spatial distribution, and the general 
physical composition of the dust grains (structure, density 
etc.). These parameters are necessary for an understanding of 
the dynamical history of the dust particles, their relation to 
cometary and asteroidal material, and their interaction with it. 

Extensive ground-based photometry and polarimetry of the zodiacal 
light has been carried out (see for example V^einberg 1964, 

Dumont and Sanchez 1975). However, observations from the earth 
alone have the fundamental drawback that the spatial distribu- 
tion and the particle scattering function can not be uniquely 
separated, unless some assumptions are made about the decrease 
of particle number density with solar distance and about the 
independence of dust composition from the position in the solar 
system (Dumont and Sanchez 1975) . Information about the valid- 
ity of such assumptions can be obtained by zodiacal light pho- 
tometry and polarimetry from space probes still in the ecliptic 
plane but far from the earth's orbit (Giese and Dziembowski 
1969 , banner and Leinert 1972 ). 
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The zodiacal light experiments on Pioneer 10/11 and Helios 
provide the brightness and polarization of the zodiacal light 
a function of heliocentric distance in the ecliptic plane, 
(Hanner ct al . 1975, Link et al. 1975) from which the large- 
scale spatial distribution in the ecliptic between 0.3 and 
3.3 A.U. can be derived. 

Impact detectors with increasing sophistication have been used 
to measure particle fluxes and velocities by experiments on tlie 
space missions of Pioneer 8/9, Prospero, HEOS 2, Pioneer 10/11, 
and Helios (Berg et al. 1971, Bedford 1975, Dietzel et al. 1973, 
Hoffmann et al. 1975, Humes et al. 1974, Soberman et al. 1974). 
Pioneer 8/9 and submicronsized impact craters on lunar surface 
samples showed the existence of a component of small particles 
moving outward from the sun under the influence of non-gravita- 
tional forces (Berg and Grun 1973, Fechtig et al. 1974). 

The interplanetary dust forms a non-stable dust cloud in the 
solar system. The continuous sources are mainly c nets, aster- 
oids, and space erosion processes. It is, however, unknown 
which of these processes contribute to which extent. The main 
dust sink is, besides impacts on planets and their satellites, 
undoubtedly the sun: dust particles spiralling into the sun 

according to the Poynting-Robertson effect (Wyatt and Whipple 
1950) presumably lead to a vaporisation near the sun (Sekanina 
1975) and hence to a dust stream of submirronsizod remnants 
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leaving the solar system as a result of the radiation pressure 
force. These so-called beta meteoroids (Zook and Berg 1975; 
Zook, in press) are also produced by mutual collisions of 
meteoroids. Hemenway has suggested that a component of these 
beta meteoroids is being formed directly at and ejected by the 
sxin (Hemenway et al 1972) . In addition to these interplanetary 
particles, a component of interstellar origin can possibly be 
expected. 


Models of Out-of-Ecliptic Dust Distribution 

To predict the particle flux and the zodiacal light brightness 
which might be observed from cui out-of-ecliptic probe, we chose 
three models for the spatial distribution of the dust, as de- 
scribed by Giese (1975) . 
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n(r) 
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(Yz) ; 
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9.0 (EllLpsoid model) 


2.5 (Gauss model) 


3.0 (Fan model) 


Here n(r) is the particle number density as a function of dis- 
tance r (A.U.) from the sun, n^ is the number density at 1 A.U. 
in the ecliptic plane, z (A.U.) is the distance above or below 
the ecliptic plane, and is the heliocentric ecliptic latitude 
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(sinB = z/r) . We have set the parameter v = 1.0, consistent 
with the Pioneer 10/11 and Helios zodiacal light data. The 
flattening parameter y was adjusted for each model to give the 
ratio 0.32 for the brightness of the zodiacal light at the 
ecliptic pole to the brightness in the plane at elongation 90^, 
as observed by Dumont (1973) . 

Contours of equal dust density are illustrated in Figures la, 
lb, and Ic respectively for the three examples. All three 
models predict that the dust is considercd^ly concentrated 
toward the ecliptic plane. In comparison, the spatial distri- 
bution of radio meteors derived by Southworth and Sekanina (1973) 
shows a similar z-dependence at a x-distance of 1 A.U., but an 
increasing number density beyond 1 A.U. in the ecliptic plane. 

The 0.5 contour in our examples ranges from z = 0.2 to 0.4 A.U. 
at X = 1.0 A.U.. The spatial density at 1 A.U. away from the 
ecliptic is at most 0.1 n^. Thus, observations of the zodiacal 
light and direct detection of particles will take place mainly 
at space probe distances less than 1 A.U. from the ecliptic. 

Direct Detection 

Direct measurements from an out-of-ecliptic mission could help 
to identify the various sources and sinks of the non-stable dust 
cloud. On the basis of the results of Pioneers 8/9 and HEOS 2 
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missions in the ecliptic at 1 A.U. and on the basis of the 

three models referred to above, the number of events per orbit 

2 

were computed for a detector surface of 100 cm and a space 
probe on a circular orbit of 1 A.U. radius having an ecliptic 
inclination of 30° or 60° or 90° respectively. Two cases of 
vehicle stabilisations (spinner with axis perpendicular to the 
orbital plane and three-axis stabilisation) were considered. 

The results are presented in Table 1 excluding or including 
the case of particles coming directly from the sun (marked 
'no sun particles* or 'sun particles'). 

From these data it is evident that, while such a particle detec- 
tion experiment is not able to differentiate among the various 
models of the spatial distribution directly, it could provide 
information on the role of the sun as a possible dust sink or 
source. Furthermore, the event rate is sufficient to continue 
analysis of enough individual particles to look for differences 
in composition between particles with low- and high-inclination 
orbits . 

Zodiacal Light 

The brightness and polarization of the zodiacal light can be 
predicted theoretically as a function of observing direction 
and observer's position in the solar system for any spatial 
distribution. If one assumes, that the average scattering 
properties of the particles are independent of location in the 
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solar system, then the brightness varration with spacecraft 

position is directly related to the spatial distribution of 

the dust (Hanner and Leinert 1972) . Models for the zodiacal 

light distribution over the sky as seen from an out-of-ecliptic 

space probe were discussed by Giese (1971). By use of the 

Scune progrcum we have computed for our three dust models the 

variation of zodiacal brightness for an out-of-ecliptic space 

probe on a circular orbit of 1 A.U. radius as a function of the 

2 distance between the spaceprobe and the ecliptic plane. The 

maximum value z of z is related to orbital inclination i by 
m 

z^ = (1 A.U.) • sin i. For the scattering function a simple 
diffraction plus isotropic reflection form (albedo = 1) was 
chosen. The size distribution adopted was a 3-part power law 
n(a)da a with < = 2.7; 2; or 4.33 in the regions of particle 
radii from 0.008 to 0.16; 0.16 to 29; or 29 to 189 wm, and 
n^ = 3.7 • 10“^^; 4.3 • 10“^^; or 7.4 • lO"^® particles/cm^ 

in each regime, respectively. This is a fair approximation 
taking into account both the distribution of particle radii as 
derived by Griin (1975) from direct measurements (Fechtig et al. 
1974) , and the brightness of 200 stars of tenth magnitude per 
square degree (200 Sj^q) found by Dumont and Sanchez (1975) at 
€ = 90° elongation in the ecliptic. 

Figs. 2a through 2c present the decrease of zodiacal light 
brightness with increasing orbital altitude z of the probe for 
the three models of Fig. 1. The viewing direction from the 
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probe is parallel to the ecliptic plane and perpendicular to 
the sun (Fig. 2a); in the positive z direction (Fig. 2b); or 
parallel to the negative x direction (Fig. 2c) . In all cases 
the models predict brightness values easily observable ( > 10 S^^) 
at space probe z distances up to 0.4 A.U. (Model II) or up to 
approximately 1 A.U. (Models I and III) / and at the same time 
observable differences between the models, particularly in the 
case of Model II with its flattened outer contours. Even if 
the density does not follow the functional forms we have chosen 
as examples, the brightness variation with space probe z dis- 
tance at a constant observing direction will give a measure of 
the rate of dust decrease away from the ecliptic plane. To 
look for any systematic changes in the average scattering proper- 
ties of the particles (size, composition) with z distance, 
polarization measurements as a function of elongation are neces- 
sary, in addition to brightness observations (see Giese and 
Dziembowski 1969) . 


Interstellar Dus t 

The existence of an interstellar component in the solar system 
dust cloud has been proposed by Greenberg (1969) and others. 
Even if such small particles are excluded from the inner solar 
system by radiation pressure, they might be observable with a 
sensitive detector at large heliocentric distance, during the 
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transfer orbit of an out-of-ecliptic probe. If we take, for 

-13 3 

example, a number density of 3 *10 particles/cm of 

interstellar origin and a relative velocity of 30 km/sec 

2 

between the particles and a 100 cm detector, 8 such particles 
should be detectable each day. 
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Conclusion 


A joint dust experiment which combines direct detection of the 
velocity and mass of individual particles with measurements of 
the zodiacal light brightness and, if possible, polarization 
will best serve the purpose of an out-of-ecliptic mission. 

These two complementary methods together will give a picture 
of the three-dimensional dynamics and spatial distribution of 
the interplanetary dust. Such a picture will help to clarify 
the relation of dust particles to cometary and asteroidal ma- 
terial as well as the interaction of the dust with the solar 
radiation field. The influence of the sun on particle motions 
and physical characteristics has a significance for the role of 
the dust in the early evolution of the solar nebula and early 
phases of stellar evolution. 
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Piq» 1 Surfaces of equal number density of dust particles 
in the interplanetary cloud 
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X solar distance in the ecliptic plane (A.U.) 
z hight above the ecliptic plane (A.U.) 


a) 

Ellipsoid Model 
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Fan Model 
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III; 
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Fig. 2 Brightness of the zodiacal light as seen '>y a spaceprobe 
ascending to an altitude of z (A.U.) above the ecliptic 
plane on a circular orbit of 1 A.U. radius. 

S^^zStars of 10th magnitude per square degree 

dashed curve: Ellipsoid-model 

dots : Gauss-model 

solid curve : Fan- model 

a) direction parallel to the ecliptic plane at right 
angle to the solar direction (e = 90°) 

b) direction towards the ocliptic pole (+ z direction) 

c) direction parallel to the •'cl.i.ptic plane and 
parallel to the earth'sun diroctioi: (- x axis). 
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A MEANS OF IN SITO MEASUREMENTS OF NEUTRAL 
H AND HE ON AN OUT-OF-THE- ECLIPTIC MISSION 
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Tucson, Arizona 85721 


ABSTRACT 

On an out-of-the-ecliptic mission, in situ ineasuremcnts of 
densities and temperature of interstellar neutral H and He in the 
heliosphere should complement observations based on backscattered 
Lyman-alpha intensities. A means of performing the in situ measurements 
is briefly described. 
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The experiments performed by the groups at the University of Paris 
(Bertaux and Blamont* 1971} and at the University of Colorado (Thomas and 
Krassa, 1971) have provided us the first glimpse of the presence and the 
distributions of interstellar neutral H and He in the heliosphere. 

Like all optical observations, the observations of neutral H and He by 
backscattered Lyman-alpha photons are indirect and integral in nature. 
Therefore, their results are constrained by the assumptions one must make 
about the properties of the solar wind and those of the solar H and He 
Lyman-alpha emissions as well as the temperature of the interstellar H 
and He. To complement the optical measurements, in situ direct determination 
of the densities and temperature of the H and He will be necessary. 

An out-of-the-ecliptic mission would provide the unique opportunity for 
such an effort since the relative velocity of the sun to the interstellar 
medium points out of the ecliptic. 

The necessary ijri situ measurements can be performed with an 
instrument of low power (1 W), low mass (< 3 kg) and single-particle 
counting capability such as the field-ionization neutral detector (FIND) 
being developed at the University of Arizona (Curtis ^ al^. , 1975). Figure 
1 illustrates the principal parts of FIND. The ionization tips, the grid 
and the surface-barrier Si detector are encased in a chamber with an 
entrance aperture. Neutral H and He enter the chamber and as they 
reach the vicinity of the tips, they are field-ionized to become and 
He^, respectively. These +1 ions are immediately accelerated towards the 
detector and their electrical signals analyzed. 

Although all the ions arrive at the detector with essentially 
identical kinetic energies acquired in acceleration, He^, being more 
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massive, will interact with the crystal lattice of the detector more and 
thus provide a smaller electrical signal than an H of same energy. 

Figure 2 is a composite plot showing that the He peak is shifted by 2.8 KeV 
from the peak which appears at 26 keV, corresponding to the accelerating 
potential of 26 kV. We note that each pulse-height distribution is gaussian 
and has well defined peak position and FWHM. Therefore, using only one 
detector and one pulse-height analyzer the two species can be separated. 

Our laboratory results also indicate that with a bundle of 200 
needles at +26 kV facing a 1 cm* detector at a distance of 1 cm, sensitiv- 
ities of 3 X 10“’ counts sec"^ per unit flux (1 unit flux = 1 ^'m"* sec“M 
for H and 4 x 10“° for He in the same units can be attained. Assuming 
an H flux of 10** cm"^ sec“S e.g., n = 0.01 cm“® and v = 10 km sec“*» 
a detector background of 1.4 x 10“* E“^‘* sec"^ keV“^ (6. Gloeckler, 
private communication, 1974) and an FWHM of 3 keV, then H signals can be 
well separated from the background in one day's accumulation. (Actual 
background can be determined in flight by turning off the high voltage 
supply to the ionization tips.) Figure 3 is a computer simulation based 
on the above assumptions. In addition to H , He of three different relative 
abundances are also included. The varied n(He)/n(H) are due to the 
different values the parameter p might take (see review by Axford, 1972). 

A least-square fit of two gaussian distributions of known peak positions 
and FWHM's to any of the three curves shown in Figure 3 will yield the 
corresponding n(H), n(He) and n(He)/n(H). 

The above description of FIND leaves little doubt that in situ deter- 
mination of neutral H and He concentrations can be performed. In addition, 
if the aperture of FIND scans the part of the sky surrounding the direction 
of maximum flux, the angular distribution of the neutral flux would then 
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be a measure of the temperature of the neutral gas at the point of 
observation. With an instrument such as FIND complementing a Lyman-alpha 
spectrometer on a spacecraft that covers large heliocentric distances and 
latitudes, the local Interstellar medium and its interactions with the 
solar wind can be examined fully. 
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FIGURE CAPTIONS 


Figure 1. An over-simplified diagram of FIND. For actual space missions, 
the ionization tips and the detector assembly should be protected 
from low energy charged particles by the use of repelling grids 
at the entrance aperture and from y-rays and secondary 
particles by the use of an anti -coincidence guard counter 
immediately behind and surrounding the ion detector. 

Figure 2. Laboratory results showing the difference in the pulse-height 
distributions of and He^ signals by a 2.8 keV shift. 

Both ion species have an average kinetic energy of 26 keV. The 
noise of the detector is 5.5 keV (FWHM). 

Figure 3. Predicted pulse-height distributions of an H flux of 1 x 10** cm'^ 
sec"* and He fluxes of 10^ , 10“ and 10 * sec”* as seen 
in one day's accumulation by a KIND having an accelerating 
potential of 26 kV, a detector noise of 3 keV (FWHM) and a 
background following the distribution 1.4 x 10"^ ^- 1-2 50 c”* keV"*. 
The sensitivies to H and He used in this calculation are 
3 X 10"^ counts sec"* per unit flux and 4 x 10“® counts sec"* 
per unit flux, respectively. 
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